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Introduction 

Guidance activity generally relates to the perturbation of spacecraft 
trajectory or state by the application of translational control effectors. 

The Space-Shuttle mission Includes three relatively distinct guidance phases: 
Atomospheric Boost (characterized by an adaptive guidance law), Extra-At- 
mospheric Activities and Re-Entry Activities (where aerodynamic surfaces 
are the principal effectors). Guidance tasks include pre-maneuver targeting 
and powered flight guidance (where powered flight is defined to include the 
application of aerodynamic forces as well as thruster forces). Figure 3«1 
is a flow chart which follows guidance activities throughout the mission from 
the Pre-launch phase through touchdown. Table 3.1 lists the main guidance 
programs and subroutines used in each phase of a typical rendezvous mission. 

A brief description of each such program and routine follows. Detailed 
software design requirements are presented in Vol. III. 

1. Atmospheric Boost Guidance 

This program, when completed, will (at a minimum) provide a programmed 
pitch over and attitude hold until sometime after max-q, It may also include 
some minimal closed-loop guidance to limit dispersion in the presence of wind 
and gust disturbances. This program (not yet submitted) will fulfill funct- 
ional requirement 2G1 5 *'. 

2. Multi-S t age Boost Guidance 

These powered ascent guidance equations provide inertial steering commands 
during boost to insertion. The equations accommodate to engine throttle cap- 
ability and to discontinous boost, i.e., PSR shutdown and jettison with con- 
trolled thrust from the MPS. These equations fulfill requirement 2G2, 

3. Rendezvous Targeting 

This Targeting Routine has the capability of constructing inflight the 
rendezvous maneuver sequence which satisfies the requirements of the particular 
mission. The routine can handle sequences with any given number of maneuvers, 
each of which can have a variety of constraints. These equations fulfill re- 
quirements 3G2. 

4. Rendezvous Braki n g 

This Targeting Routine has the capability to bring the vehicle into the 
station-keeping zone by automatic line-of-sight corrections and braking cor- 
rections. These equations fulfill requirements 3G3. 

^Functional Requirement Module (See App. I). 
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r t si.qt ion Keepi ng Guido e ci. 

Tliis targeting Routine is for use during the station-keeping phase and 
is designed to maintain 1, orbiter in a small zone which may be arbitrarily 
located with respect to t' - Garget vehicle. The orbiter is maintained in the 
non by the periodl c . of small velocity corrections computed so as 

w two avoiigo v_~: i^iooi'e of propellant per orbit. These equations 



This targeting Routine is for the computation of an optional phasing 
maneuver to place the veht rip. in a phasing orbit prior to ^Landing, and an 
in-plane minimum fuel deorbiu maneuver satisfying entry-interface and landing- 
site constraints. The propram is designed to allow the crew to determine 
; ' . il": deorblu t . 4J . : and to select one desired. The program 

satisfies requirement 31 6, 

7. Powered Flight Guidance 

These Guidance Routines are for use in computing steering and engine 
cutoff commands during either a maneuver with a specified velocity change, 
or a Lambert aim point maneuver, or a deorbit maneuver. The concepts of a 
curr- ^ positir. of ft--. n .1 state-vector navigation are used for the 

Lrim: - , and deo it i V'.- Cross product steering is used for all man- 

euvers. These equations satisfy the requirements 3G1. 

8. Entry and Transition Guidance 

This routine will provide guidance commands from entry interface through 
the heat control phase, through transition, to 40*000 feet altitude, The 
requirements for this phase are 4G1. 

9. Approach (Terminal Area) Guidance 

This routine provides steering commands which guide the S/S from an 

It- ti’de of 40,000 ft to s^-oint on the final approach (glide) path at an 
• -.XL "v v.ude of approximately . .. 0 ft. These requirements are listed as 402. 

10. Final Approach Guidry ce 

This routine provides steering commands which maintain the shuttle on 
the two-flare approach path through touchdown. These equations satisfy re- 
quirements 4G3. v . 

11. Conic State Extr a no" 1 ; tion Subroutines 

These subroutines are for conic state vector extrapolation as a function 
of time (Kepler) or as a function of angle (Theta), and are required both for 
raid- -ice targeting and for "1 -igation. 
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1 2 . Pr ecision State and Fi Iter We i ghtin g Matrix E xtra polation 

This subroutine has an Encke integration scheme wJiich includes the 
capability for precision extrapolation of a vehicle state vector and the 
associated submatrix of the Navigation filter weighting matrix in the 
earth's Jp gravity field. Additional perturbing accelerations due to 
higher order gravity terms, lunar and solar gravity, and atmospheric 
forces have not been included since the requirements for them have not 
been established. This subroutine is presented in Navigation Volume II. 

13. Conic Required Velocity 

This subroutine is for the solution to the multi-revolution Lambert 
required velocity determination problem. 

14. Precision Required Velocity 

This subroutine is for use by a targeting routine to compute the 
parameters needed by the Powered Flight Guidance Routines to perform a 
Lambert aim point or de-orbit maneuver. 

15. Abort Guidance Targeting 

This program, available only In preliminary form, provides guidance 
targeting for the transition from booster failure to acquisition of the 
nominal entry trajectory with virtually empty OMS fuel tanks. This pro- 
gram fulfills some of the requirements of 2G3. 

The Multi-Stage Boost Guidance presentation is by R, F. Jaggers of 
the Boeing Co. There are other designs under consideration, but none differ 
significantly from the Linear Tangent Guidance which is the basis of Jagger's 
presentation. 

The On-Orbit guidance submittals are the work of G. S. Draper Laboratory. 
Rendezvous Targeting is documented by V. H. Templeman, Rendezvous Braking by 
P. M, Kachraar, Station Keeping Guidance by Gustafson and Kreigsman, Deorbit 
Targeting by Brand and Brennan, Powered Flight Guidance by Brand, Brown, 
Higgins and Pu, the Precision Required Velocity routine by T. J. Brand and 
the Conic State routine, the Precision State routine and the Conic Required 
Velocity routine by W. M. Robertson. An alternative Rendezvous design has 
been published by D. J. Jezewski of the NASA/MSC Mission Planning and Analysis 
Division, It is anticipated that the Linear Tangent Guidance, in a modified 
form, will replace the Powered Flight Guidance design documented herein. 
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Figure 3. 1 Guidance 
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SUBROUTINES 
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/ / / 

TERMINAL 
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* RDZ BRAKING 
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/ / / / 

/ / / 
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/ / / 

DEORBIT 
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/ / / 

ENTRY 

LANDING 

ENTRY GUIDANCE 
APPROACH GUIDANCE 
FINAL APPR. GUID. 



"targeting and powered flight guidance 

USED SEQUENTIALLY FOR EACH MANEUVER 


TABLE 3. 1 SHUTTLE GUIDANCE 






The Entr y through Landing guidance submittals are produced by Kriegsman (CSDL) 
r>r>i jj.arpold. (KPAD) for Entry, NASA/MSC Crew Procedures Division (documented by Trio) 

«- .i Elias (CSDL) for Approach, and D. Dyer of GCD for Final Approach. Other Entry 
guidance designs include Fast Time Integration by Sunkel (G&CD) . The Harpold design 
r-nd the Kriegsman design for Entry Guidance are both included in this volume because 
both have outstanding features and it is conjectured that the final design will 
represent a combination of the best features of each. Likewise both the Crew Pro- 
cedures and the Elias designs are presented in expectation that the final design 
v iH include features of each. An alternate Terminal Area Approach design by T. Moore 
(GCD) is more complicated, but is a strong contender for implementation because it 
includes capability for low-altitude redesignation. 

The preliminary Boost Abort Targeting submittal is by G, McSwain Oj. G&C Div., 
HASA/MSC. 

Except for Boost Abort Targeting all Guidance submittals are complete, according 
to the requirements of Appendix I. 




Submittal 53 Multi-Stage Boost Guidance 

Introduction 

The purpose of the powered ascent guidance equations is to provide inertial 
steering commands during the boost to orbit maneuver. Throttle setting can also 
be provided if this is a requirement of boost guidance* however* this feature is not 
incorporated in the equations presented here. 

The linear tangent guidance law presented here was developed to meet at 
least the following requirements: multi-stage capability, ability to handle flight 
perturbations and maintain orbital insertion accuracy, abort to alternate conditions, 
engine out capability, and throttle for constant acceleration. 
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C o i r« ] '• ; : 1 at i o n a I Flow of Multi-Stage Guidance E qua ti ons 

Guidance input and output parameters are listed and defined below. Essen- 
l: }] Vl guidance input is the present slate vector and orbital parameters defining 
l!(f , desired terminal state vector, and the output is steering commands and time- 

to-go. 


Guidance 


R 


D 


} 

D 


G 


21 


C 23 


ex. 

l 


r. 

i 


T 

1 Bi 


T Ci 


a Li 


n 


Presettings (PREFLIGTIT IN P UT) 

Magnitude of desired terminal radius vector 

Magnitude of desired terminal velocity vector 

Desired terminal flight path angle 

X component of unit vector normal to desired orbit plane" 

Z component of unit vector normal to desired orbit plane" 

Effective exhaust velocity of stage i 

Initial m to m ratio of stage i 

Nominal burn time of stage i 

Coast time between stage i and stage i + 1 

Acceleration limit of stage i 

Number of guided stages 


Navigational Quantitie s (INFLIGHT MEASURED INPUT) 
a Inertial platform measured acceleration vector 


R 


Inertial platform velocity vector 
Inertial platform radius vector 


Inertial platform gravity vector 
(Calculated function of R ) 


* Optional inputs are i^, desired inclination* and Ap, desired longitude of des- 
cending node. 
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G uidance Output 

0 Commanded inertial pitch angle 

Commanded inertial yaw angle 
T gq Time-to-go till orbital insertion 

f Throttle setting if applicable (to be determined) 

Guidance Precalculations 

(Coordinate transformation from platform system to desired orbit system) 

Input anc * ^23' ^ anc * ^ components of unit vector normal to desired orbit plane. 
Compute unit vector in desired orbit plane normal to launch vertical, and unit vec- 
tor defined by intersection of desired orbit plan and the plane containing launch 
vertical and vector normal to desired orbit. 


°22 - (1 ' C 21 2 - G 23 2 > 


G 11 " ^ G 22 + G 23 * 


G 32 = ' G 23^ G 11 
G 33 11 G 22 /?G 1 1 
G 12 = " G 21 G 33 
G 13 = G 21 G 32 


NOTE: For due East launch 
and no plane change 


[G] = 




• J *Vt* iU* k* 


ENTER GUIDANCE PRE-TURN ON 




/fr TwcCTMrer'TgrararsnLUi 


sinv^cosA^l r s 1 in:^ 


, s1n \ ” sini^cose^ 


cos^cosA^J 


r . (r ? + r 2\ H 

G n - (G 22 ^ 23 ) 


G 31 - 0 


C 32 ' " G 23 /G ll 


1 •. , ’ G 22 /G n 

G J2 ‘ ” G 21 G 33 
''...'■ G 21 G 32 


DESIRED ORBIT 

COORDINATE 

SYSTEM, 

G MATRIX 


f* T&wiv i<jv f-u dpaf’ar <>¥ ^ 


i = 1 
T G0 = 0 


T 11 * t bi 
u ri * G 2i 

T G0 “ T G0 + T Bi + T C1 


i = i + 1 <$ 



= Launch Azimuth 
= Geodetic Launch Latitude 

= Desired orbit inclination 

= Desired Longitude of descending 
node (measured from launch 
meridian) 


EXIT. TO GUIDANCE 
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ENTER GUIDANCE {MAJOR CYCLE LOOP) 



INITIAL RANGE ANGLE 


LOCAL GUIDANCE 
COORDINATE SYSTEM 


RADIUS AND VELOCITY IN 
LOCAL GUIDANCE FRAME 

MEASURED ACCELERATION 
MAGNITUDE 


AVERAGE RADIAL GRAVITY 
MAGNITUDE 
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Time till acceleration limit 

TIME FROM LIMIT TO BURN OUT 
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Submittal 3 7: Rendezvous Targeting 


1 . INTRODU CTION 

The rendezvous of the Orbiior (primary vehicle) with a target vehicle (e.g. 
the Space Station) is accomplished by maneuvering the Orbiter into a trajectory 
that intercepts the target vehicle orbit at a time that results in the rendezvous of 
the two vehicles. The function of rendezvous targeting is to determine the targeting 
parameters for the powered flight guidance for each of the maneuvers made by the 
Orbiter during the rendezvous sequence. 

In order to construct the multimaneuver rendezvous trajectory, sufficient 
constaints must be imposed to determine the desired trajectory. Constraints 
associated with the Orbiter mission will involve such considerations as fuel, light- 
ing, navigation, communication, time, and altitude. The function of premission 
analysis is to convert these — which are generally qualitative constraints-into a set 
of secondary quantitative constraints that can be used by the onboard targeting 
program. By judicious selection of the secondary constraints, it should be possible 
to determine off-nominal trajectories that come close to satisfying the primary 
constraints. 

The proposed onboard rendezvous targeting program consists primarily of a 
main program and a generalized multiple-option maneuver subroutine. The driving 
program automatically and sequentially calls the maneuver subroutine to construct 
the rendezvous configuration from a series of maneuver segments. /The main pro- 
gram is capable of handling rendezvous sequences involving any given number of 
maneuvers. Enough different types of maneuver constraints are incorporated into 
the subroutine to provide the flexibility required to select the best set of secondary 
constraints during premission planning. In addition, the astronaut has a large, 
well defined list of maneuver options if he chooses to modify the selected nominal 
rendezvous scheme. 

As the new approach represents, in essence, just one targeting program, 
there is considerable savings in computer-storage requirements compared to former 
approaches in which each maneuver used in the rendezvous scheme had a separate 
targeting program. The programming and verification processes of this unified 
approach will also result in implementation efficiencies. 

1. l Number of Independe n t Constraints Involved in a 

Rendezvous Sequence 

During the Gemini and Apollo flights and in the design of the Skylab rendezvous 
scheme various numbers of maneuvers were utilized in the rendezvous sequence. 

The range went from two (Apollo 14 and 15) to six (Skylab). 
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The number of independent constraints (i. <?• . the number of explicitly satisfied con- 
straints) in each rendezvous sequer * must equal the number of degrees of fieedom 
implicitly contained in the sequence' To establish this number, a rendezvous con- 
figuration can be constructed by.,”- .^0’" arbitrary, constraints until the configura- 
tion is uniquely defined. For r v- : >■ ^our maneuver coplanar sequence is 
shown in Figure 1. followed by tv . A - terminal point. Using the constraints 

v (velocity magnitude), r. easy to establish that the total number 

involved is 12, assuming the lime of the first maneuver has been established. Re- 
moving one maneuver will reduce the number of degrees of freedom by three. 
Hence, the number of independent constraints necessary to uniquely determine the 
maneuver sequences are 


Number of maneir * * 
in sequence 


Number of independent 
constraints required 


If the above rendezvous are not coplanar, one additional constraint has to be added 
to each sequence to allow for the oir vf-plane component. 

In some cases the number of * "'y constraints may be insufficient to 

uniquely determine a rendezvous tr ctory for the desired number of maneuvers. 
One way of overcoming this deficiency in constraints is by introducing sufficient 
variables to complete the determination of the rendezvous trajectory and then 
determining values for these variables by minimizing the fuel used. 

In order to take advantage of updated state vectors due to navigation or ground 
updates, the rendezvous targeting program is called prior to each maneuver to 
compute the upcoming maneuver. In general, each maneuver computation will in- 
volve a multimaneuver sequence as the nature of the targeting constraints do not 
allow the maneuvers to he independently computed. The relationship between the 
rendezvous sequence involving n ttV /ers and the maneuver sequences is shown 
below. Maneuver Segments 


Rendezvous Sequence 


Maneuver Sequences 


(maneuver points) 
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K.-ich maneuver s ■ ■ . is composed of a number of maneuver segments and is 

basically indepei om the other maneuver sequences. These sequences must 

have the same n : * v T - independent constraints as tabulated above. 

1.2 The Corr - : rue* mt? of a M a n e u v er S e gme nt 

Each n-i., . * ■ aezvous sequence can be divided into n-1 maneuver seg- 

ments. Each segm .rTa . ; s, basically, the addition of a maneuver to the primary 
vehicle’s velocity w. -ar mid an update of both vehicle's state vectors to the next 
maneuver point. 

A maneuver segment is herein generated in one of three ways: 

Forward’ generation A maneuver 4v is computed and added 

to the velocity vector in a specified di- 
- rection. The state vector of the primary 

vehicle is then updated through a speci- 
fied amount to arrive at the next maneu- 
ver position. 

Target generation The target vehicle is updated through a 

specified amount and then offset to establish 
a target vector. An option is available at 
this point to compute a coelliptic velocity 
vector and update through A t to establish 
. . a new target vector as shown below. 





The maneuver is then computed by uniquely 
specifying the nature of the traverse between 
the primary vehicle’s position and the target 
vector. 


1 

A 
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I nteg r ate d gene ration In this case, the maneuver segment 

is computed as an integral part of a 
maneuver sequence involving more 
than one maneuver segment. The 
nature of the constraints are such that 
the maneuver sequence cannot be sub- 
divided’ into uniquely defined maneuver 
segments. The maneuver segment will 
usually have one degree of freedom, 
which will generally be assumed to be 
the magnitude of the maneuver. 

Each of the above methods is defined by specifying trajectory constraints 
by setting certain switches and parameter values. Specifying a trajectory constraint 
is equivalent to specifying one or more independent constraints. On the other hand, 
specifying an independent constraint can also be equivalent to specifying one or more 
trajectory variables. (See Ref, 9) A trajectory constraint common to all three of 
the above methods is the state vector update switch s ■ The options associated 

with this switch are: 


1 

2 


s „ = < 

update 


3 

4 
b 


Update from time t to time tp 

Update through time interval At 

Update through n revolutions 

Update through 8 radians 

Update to be colinear with a 
specified position vector 


In the next three sections, the trajectory constraints associated with each of the 
above methods will be listed. 
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1.2. 1 Mane uver O pt ions in Forward Gen eration of 

Maneuver Segment 

The forward generation of a maneuver segment is accomplished in one of two 
ways. Either the maneuver magnitude is uniquely determined in terms of the state 
vector at the maneuver time or the maneuver is determined by an iterative search 
to satisfy a terminal constraint. 

The maneuver magnitude Av is either calculated or assumed depending on 
the maneuver switch s , and it is applied in a direction controlled by the di- 


rection switch s dlrect . 


, and it is applied in a direction controlled by the di- 
an 

The options associated with the maneuver switch are: 


1 A v is assumed specified 

2 Av is computed based on a post maneuver 
velocity vector being "coelliptic" with the 

s = J state vector of the target vehicle 

man A 

3 A v is computed from the conic circular 
velocity constraint 

4 Av is computed based on a Hohmann type 
transfer resulting in a Ah change in 
altitude 

The options associated with the maneuver direction switch are: 

-1 Apply Av is horizontal direction in plane 

of primary vehicle 

1 Apply Av in horizontal direction parallel 

J to orbital, plane of the target vehicle 

direct _2 Apply Av along velocity vector in plane 

of primary vehicle 

2 . Apply Av along velocity vector parallel 

^ to orbital plane* of the target vehicle 


The selection of the update switch s 


update 


determines the update of the 


primary vehicle f s trajectory following the maneuver to the position of the next 

maneuver. A terminal constraint can be imposed at this point by setting the ter- 

minal switch s * : 

term 
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j 1 Terminal constraint is a height constraint 

y-l Terminal constraint is a phasing constraint 



Following the computation of the height/ phasing error, the maneuver magnitude is 
varied in an iterative search to satisfy the height/phasing constraint, 

1. 2. 2 Maneuver Options in Target Generation 
of Maneu v er Sc g me* t it 

The target generation of a maneuver segment starts with the selection of 

the update switch for the target vehicle. If this switch equals four, 0 will be 

augmented by the central angle between the primary and target vehicles before 

being used. The position of the target vehicle is then offset through either (e^, Ah) 

or (A0, Ah), depending on whether is negative or positive, to obtain a target 

vector. The "TPI offsets" (ej , Ah) are discussed in Section 5 (see Figure 2 for 

definition of e r ). If | s | equals two, a coelliptic velocity vector is computed 
ij tar 

based on the target vector, and a new target vector is defined by updating the co- 
elliptic state vector through At . The options associated with s tar are: 



Offset target (e^, Ah). Compute coelliptic 
velocity and update through (negative) At . 

Offset target (e j , Ah) 

No target offset 
Offset target (A0, Ah) 

Offset target (A0,Ah). Compute coelliptic 
velocity and update through (negative) At 


The nature of the traverse between the primary 

vehicle's initial state vector and the target vector is 

controlled by the maneuver switch s*„ „ : 

J man 



Lambert - the trajectory is time constrained 

Horizontal - the' maneuver is constrained to 
be in the horizontal direction 

Tangential - the maneuver is constrained to 
•Be in the direction- of the velocity vector 

Apogee/Perigee - the trajectory has' an apogee/ 
perigee occurring at the target point. 


There is a minimum Av option associated with the above maneuvers which 
is controlled with the optimization switch s ^ : 
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| = Unit horizontal in forward direction for primary vehicle 

LOS = Line of Sight 

1. If the LOS projection on ji is positive-. 

a. When the LOS is above the horizontal plane. 0 < e|_ < 7 r/ 2 

b. When the LOS is below the horizontal plane, 3 tt/ 2 < ej_< 2rr 

2. If the LOS projection on Ms negative: 

a. When the LOS is above the horizontal plane, Til 2<ei<7r 

b. When the LOS is below the horizontal plane, tt < e[_ < 37 t/ 2 


Figure 2, Definition of the Elevation Angle e^ 
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“2 Minimize the sum of Ihe magnitude of 
the first and the next maneuver (based 
on a coelliptic parting velocity) by 
varying At , the time of update of the 
target vehicle. 

1 Minimize ihe magnitude of the first ma- 
neuver by varying £t, the time of update 
of the target vehicle. 

s = < 

1 Minimize the magnitude of the first ma- 
neuver by varying Ai . the time between 
ihe next maneuver and the initial offset 
position. (See sketch on page 1-4 ) 

'J. Minimize the sum of the magnitudes of 
the first and the next maneuver (based 
on a coelliptic parting velocity) by 
varying At , the time between the next 
maneuver and the initial offset position 
(see sketch on page 1-4 ). 

This minimization is accomplished by driving the slope ( Av / independent variable) 
to zero using a Newton Ruphson iteration scheme. 

^•2. d Maneuver Options i Integrated Generation 
of Maneuver Segn.em 


The integrated generation of a maneuver segment involves an iterative solu- 
tion to determine a maneuver sequence which cannot be sequentially solved for its 
maneuver segment components. The maneuver is computed by guessing its magni- 
tude, assigning a direction and plane through selection of the direction switch 

S dirert ' u Pdaling the primary vehicle* s state vector after selecting switch s 

c 1 6 update 

and then calling additional maneuver segments until reaching the point at which P 
the terminal constraint is to he attained. The maneuver is then iteratively de- 
termined by satisfying the terminal constraint. The number of additional maneu- 
ver segments and the nature of the terminal constraint are controlled by the ter- 
minal constraint switch s 
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2, -3, The terminal constraint is a phasing 

constraint and it occurs at the I s 

J term 

maneuver point from the start of the 
maneuver segment . 


s 

term 




2, 3, . The terminal constraint is a height 
. . <10 

constraint and it occurs at the s 

term 

maneuver point from the start of the 
maneuver segment. 


^ 10< S term Both a height and phasing constraint 

^ occur at the same maneuver point. 

The first digit n of s renre - 

1 term 1 

sents a phasing constraint that occurs 
at the n maneuver point from the 
start of the phasing maneuver segment. 
The last digit of s term represents 
a height constraint that occurs at the 
°2 maneuver point from the start of 
the height maneuver segment. 


1,2,4 Summary of the Maneuver Segment Con straints 

The maneuver and trajectory constraints that can be imposed on a maneuver 

segment can be divided into the following catagories (see Figure 3). 

Primary vehicle update constraints 
Target vehicle update constraints 
Initial velocity constraints 
Offset constraints 
Terminal constraints 
Traverse constraints 


Table 1 contains a detailed listing of the constraints. The three independent con- 
straints (four in the case of noncoplanar traverses) which govern a maneuver seg 
ment cannot be chosen arbitrarily from this list for two reasons: 

(1) There is not a one-to-one correspondence between the 
trajectory constraints and the independent constraints. 
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(2) Selecting some constraints negates the need for some 
otiiers (e.g. selecting a Lambert constraint negates the 
need for a maneuver direction constraint). 

In the case of a straight forward rendezvous profile, a basic understanding of the 
na ure of the constraints should allow the constructor of the rendezvous sequence 
o e oose a set of trajectory constraints which determine the required number of 
in epen ent constraints. For a complex rendezvous profile, such as Skylab a 
more formal approach such as presented in Reference 11 should be used. One of 
16 JUSUflcatlons for Presenting the three methods of generating a maneuver seg- 
ment was to aid the constructor of the rendezvous sequence in choosing compatible 
sets of constraints. ^ 
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TABLE 1 
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DETAILED LISTING OF CONSTRAINTS 
(Sheet 1 of 2) 


Primary and Target V ehicle Update Constraints 
Delta time 

Initial and final time 
Central angle 
Number of revolutions 
Terminal position vector 


Initial Velocity Constraints 
Plane 



Parallel to target orbit 
Parallel to primary orbit 
Direction 

Horizontal 

Along velocity vector 
Magnitude 

Circular 
Coelliptic 
Altitude change 
Specified 

Offset Constraints 

Angle 

Altitude 

- Elevation angle 

Terminal Constraints 
Height 
Phase 
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TABLE 1 


DETAILED LISTING OF CONSTRAINTS 
, (Sheet 2 of 2 ) 

Traverse Constraints 
Minimum Fuel 

One maneuver optimization 
Two maneuver optimization 
Apogee/Perigee designation 
Horizontal maneuver 
Tangential maneuver 
Lambert (time) 
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NOMENCLATURE 


i- ' axis of a conic 


c l’ c 2’ c 3 


Alarm code i 

Failure in fuel optimization loop 

Failure m height loop 

F , F i ; t ;-j phasing loop 

Fa re in obtaining Lambert solution in 
General Maneuver Routine 

Failure to find perigee/apogee in Search Routine 

Incompatible altitudes and elevation angle 

Fan are to find time corresponding to elevation 
amf, • . "earch Routine 

‘V,v- ' v o find desired position vector in 
Desired Position Routine 

Failure to update through 0 in Update Routine 

Iteration counter 

Height iteration counter 

Phase iteration counter 

R £ v diate variables 

Dciia altitude 


Delta projected position 




Delta time 


t 


4 


) 


At 
At 
Ay 
A v 
A v 

Av h 

Av 

P 

A V r 

A x 
A 6 


max 


LOS 


LV 


l N 


max 


m 


M 


Maximum time step allowed in Search Routine iteration 
Maneuver velocity 

Maneuver in line-of-sight coordinates 

Maneuver in local vertical coordinates 

Av used during height maneuver 

Av used during phasing maneuver 

Delta velocity used in fuel minimization loop 

Delta independent variable 

Delta central angle 

Error 

Eccentricity 
Height error 
Phasing error 

Elevation angle (defined in Figure 2 ) 

Unit vector 

Unit normal to the plane used in powered flight 
guidance 

Number of the maneuver 

Number of the last maneuver in the 
maneuver sequence 

Estimated vehicle mass 

Rotational matrix 

Vector normal to the orbital plane 

Number of revolutions 
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5 

\ 


| 

} 

i 

4 

* 





| 



p 


Partial used in Newton Raphson 
iteration 


Distance ratio 
Position vector 



s 

s 

astro 

s ^ 
coplan 

S direct 
s 


eng 


s 

s 

s 

s 

s 

s 

s 

s 


exit 

fail 

man 

opt 

outp 

pert 

phase 


proj 
S rdes 
S soln 

s search 

5 * 
tar 

s term 
S update 


Desired position vector 

Target vector used in powered flight guidance 

Switch used in Desired Position Routine 

Astronaut overwrite switch 

Coplanar switch 

Maneuver direction switch 

Engine select switch 

Program exit switch 

Failure switch 

Maneuver switch 

Maneuver optimizing switch 

Out-of-plane switch 

Perturbation switch 

Phase match switch 
Projection switch 
Desired position switch 
Solution switch 
Search switch 

Target offset switch 
Terminal constraint switch 
Update switch 


S37-17 











t^Slirttiiilig!i 'aid *&» W £&•» tkiinb tit I ..v B- - 1 


1 F 

: d time 

V 

' ity vector 

v v 

Vertiral component of velocity 

V 

c 

ir velocity 

X 

Independent variable in Iteration Routine 

y P > y P . 

Out-of-plane parameters (see Figure 6a) 



y t 


o 

p 


Radial component of velocity divided by v^, 
R - -ontal component of velocity divided 



e i 
e 2 
C 3 
f 4 
f 5 
f 6 
f 7 
f 8 


Y 

P 

6 


8 

P 


by v c 

Tolerance on time in fuel optimizing loop 

Tolerance on height in height loop 

Tolerance on central angle in phasing loop 

Tolerance on central angle in Search Routine 

To 1C ranee on elevation angle in Search Routine 

Altitude increment in Search Routine 

Tolerance on central angle in Desired Position Routine 

Tolerance on central angle in Update Routine 

Flight path angle 

Gravitational constant 
Central angle 
Perigee angle 
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Subscripts 



I 7 Final 

i Number of the maneuver 

I Initial 


LOS 

LV 

N 

0 

P 

S 

T 

TA 


Li no -of- sight 
Local vertical 
N e vv 
Old 

Primary 

Stored 

Target 

Target for primary vehicle 
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FUNCTIONAL FLOW DIAGRAMS 


1 

c? 


.s! 

i 



1 





i 


4 



iV 

'i 

4 


■1 


4 


-a 


The rendezvous targeting program consists of two major parts — a generalized 
maneuver subroutine which basically computes a maneuver and updates the state 
vectors of both vehicles to the time of the next maneuver and a main program which 
sequentially calls the subroutine to assemble a rendezvous sequence. These pro- 
grams call a number of subroutines which are briefly described below and in de- 
tail in Section 5. 


Search 


Phase 

Match 


Desired 

Position 


Update 


Coelliptic 

Mane uver 


Iteration 


To update the state vectors to either a 
specified apsidal crossing, a time, or 
an elevation angle. 

To phase match the target vehicle's 
state vector to the primary vehicle 1 s 
position vector. 

To compute an offset target vector or a desired 
position to be used in a phasing constraint. 

To update a state vector through a speci- 
fied interval. 

To compute a coelliptic velocity vector. 

To determine a new estimate of the in- 
dependent variable in a Newton Raphson 
iteration scheme. 


The functional flow diagram for the main program is shown in Figure 4. The 
main function of this program is to sequentially call the General Maneuver Routine 
to compute each maneuver segment for maneuvers numbered from i to i 
T ere are three major options that can be exercised prior to the calculation of the 
first maneuver segment: 


(1) A search for the time of the first maneuver. 
This time can be specified by: 

(a) An elevation angle, which is to be 
attained at the maneuver time. 
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(Ij) Whether the next maneuver should occur at 

the next apsidal crossing, the next perigee . 

crossing or the nth apsidal crossing. 

(2) A phase matching of the state vector. 

« 

(3) A rotation of the primary vehicle' s state vector 
into the plane of the target vehicle. 

There are three separate iterative loops built around the call to the general 
maneuver routine. One loop serves to minimize the fuel used during a maneuver 
segment with the options determined by the optimizing switch. 

The other two iterative loops involve maneuver segments which contain con- 
straints that do not allow the explicit calculation of the maneuver. These con- 
straints are height and phasing constraints imposed at the end of a maneuver seg- 
ment and controlled with the terminal switch. The iterative loop will involve 

several maneuver segments if sufficient constraints are not imposed to solve each 
segment uniquely. 

The functional flow diagram for the general maneuver routine is shown in 

Figure 5. This routine generates the departure velocity at the initial point in one 
of two ways: 

(1) As an explicit function of the initial state vectors. 

(2) By defining a target vector and then computing an 
intercept trajectory based on a specified constraint 
(as indicated by the setting of s man )« The target 
vector is determined by offsetting the updated posi- 
tion vector of the target vehicle. Depending on the 
setting of the switch s tar , a coelliptic velocity vec- 
tor is computed at the offset point and the coelliptic 
state vector is updated through At to obtain a target 
vector. 

Following an update of both vehicle's state vectors to the time of the next maneu- 
ver, the Av used or the terminal height/phase errors are calculated as required. 
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■'* • I NPUT A N D O UT P ( I T V A R Y A B LF S 

Ihe inputs to the orbiter rendezvous targeting program can be divided into 
five categories. 

Pro -Maneuver S w itches 

Upon selecting a maneuver from the rendezvous sequence, these switches 
(specified for each mane uver} serve in determining the state vectors at the maneu- 
ver point, the out-of-plane parameters and the calculation of a desired position 
vector. These inputs can also be used in determining the time of a specified apsidal 
crossing or the time at which a specified elevation angle is to be attained. 


’coplan 


*exit 


*outp 


*pert 


’phase 


■< 


Coplanar switc h 

0 Bypass 

1 Rotate primary state vector : nto plane 
of target vehicle's orbit 

Exit sw i t c h 

0 Bypass 

1 Exit from routine 
Out-of-plane switc h 

0 Bypass 

1 Compute out-of-plane parameters 
Compute out-of-plane parameters and 
modify maneuver by -y 

Perturbation switc h 
Do conic state vector updates 
Include oblateness based on J ^ 

Other perturbations as required 

Phase match switch 

0 Bypass 

1 Phase match state vectors (target leading 
primary) 

2 Phase match state vectors based on target 
leading primary by more than 360° 

1 Phase match state vectors (primary 
leading target) 

-2 Phase match state vectors based on primary 
L leading target by more than 360° 
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Desired position switch 


V 


j* 



■A 


1 


1 


i 

* 

1 

i 




i 

i 







S rdes 


s scarch 



By pass 

Compute desired position vector based on (e f , Ah) 
Compute desired position vector based on (A0. Ah) 


Search switch 

r 

-4 Compute elevation angle 
-3 Search for elevation angle 
J “2 Update to time t. 

-1 Search for next perigee crossing 
0 Bypass 

n Search for the nth apsidal crossing 
^ (n >0) 


Maneuver Switches 


These switches (specified for each maneuver) set the constraints employed 
in determining the maneuver segments. 


Direction switch ‘ 

•- — - 

Av in direction of primary's velocity 
vector, parallel to primary's orbital 
plane 

-1 Av in horizontal direction, parallel 

^direct " \ to Primary's orbital plane 

0 Bypass 

1 Av in horizontal direction, parallel to 
target's orbital plane 

2 Av in direction of primary's velocity 

w vector, parallel to target's orbital plane 

Maneuver switch 


s 

man 


1 Av is specified 

2 Av is based on coelliptic velocity 

3 Av is based on circular velocity 

4 Av is based on altitude change 

5 Lambert maneuver to offset target vector 

6 Horizontal maneuver to offset target vector 
^ Tangential maneuver to offset target vector 

L 8 p ^rigee/ apogee insertion at offset target vector 
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ti 



.ijneuvcr optimizing s witch 

1 : ■'ass 

* dnimize Av[ 

2 Minimize Av. + Av. , . 

i i + 1 



s 


soln 



'. . i . ; s . ; ■ evoluti on so l ution switch 

Solution with smallest initial flight path 
angle (measured from local vertical) 
Solution with largest initial flight path angle 


s 


tar 



Offset target (e L> Ah). Compute coelliptic 
velocity and update through (negative) At 
r « target (e ^ , Ah ) 
i Lrget offset 
Offset target ( A0 , Ah) 

Offset target ( A0 , Ah). Compute coelliptic 
velocity and update through (negative) At 


s term 


Tc rm i nal c o nstraint switch 

n (<0) Compute phasing error and back up 
fn + 1} maneuvers for start of phase loop 

0 ^ 1 1 '■ ss 

n . ' u^iO) Compute height error and back up 

n - 1 maneuvers for start of height loop 

n (10 < n < 1 00 ) Phase and height loop terminate 
on same maneuver. For phase loop back up 
x - 1 (where x is first digit of n) maneuvers 
for start of phase loop. For height loop back 
up y - 1 (where y is last digit of n) maneuvers 
for start of height loop 


s 


update 


Update switch 
0-,'! -pass 

1 Update through t p - t 
’ 2 Update through At 

3 Update through n r 

4 Update through 0 

^5 Update to be colinear with 





Parameter Values 

The parameter values (specified for each maneuver) are values for the con 
strained parameters. 
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A h Della alt it tide 

A h p Delta altitude, final 

A 0 Della central angle 

At Della time 

A v Maneuver magnitude 

n Number of revolutions 

t p Final time 

e^ Elevation angle 

Post- Maneuver Switch 

This switch determines the options available following the calculation of the 
maneuver. 


Astronaut overwrite switch 


s 


astro 


0 Bypass 

1 Overwrite maneuver in local vertical 

coordinates 

2 Overwrite maneuver in line-of-sight 

coordinates 


Maneuver Call Variables 

The maneuver call variables have to be specified for each call to the maneu- 
ver sequence. 

r p, Vp State vector of the primary vehicle 

r p, Vj State vector of the target vehicle 

i Maneuver number 


t Current time 

_ , . .th 

t . Time of the i maneuver 

l 

s Engine Select switch 

eng 

m Estimated vehicle mass 

Depending on the rendezvous sequence, there may also be some switches that have 
to be modified as a function of the maneuver number. 

Excluding the maneuver call variables, all the input variables can be set 
prior to the flight. 

The output parameters for the initial maneuver in the sequence are more 
complete than for the succeeding maneuvers. 



Output Parameters for the Initial Maneuver 


Av 


Av 


LOS i 


A ^LYi 


- 1c 


■N 


Maneuver magnitude 

Maneuver in line of sight coordinates 

Maneuver in local vertical coordinates 

Target vector used in Powered Flight Guidance 
Routine (See Ref, 5) 

Unit normal to plane used in same routine 


Other parameters such as delta altitude, phasing angle, elevation angle and perigee 
altitude can be computed as required. 

Output Parameters for the Other Maneuvers in the Sequence 

t Time of the maneuver 

Av Maneuver magnitude 


Illustration of Inputs 

Table 2 contains a set of inputs for the Orbiter targeting program based on the 
five maneuver Skylab rendezvous configuration. The following switches and para- 
meters are not used as inputs to the Orbiter program: 

s astro* S exit' S opt’ S outp ‘ S soln ’ 1 


The inputs in Table 2 are set prior to the mission so they will not have to be in- 
serted by the astronaut. The astronaut will have to modify the following quantities 
upon resetting the maneuver number as well as inserting the time of the next ma- 
neuver. 


i = 2: 

s . =0, 

term^ 

s, =32 

term . 

4 

i = 3: 

s = 5, 



man^ 

At 3 = ~ At NSR-TPI 

1-4: 

s , =0 



term^ 
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| INPUT VARIABLES FOR SKYLAB RENDEZVOUS CONFIGURATION 

4 

I Maneuver 


-a 




4. DESCRIPTIO N Ol EQUATIONS 

The only eqmui - -ntained in this document which are not trivial are those 
involved in computing me traverse between two specified position vectors. The 
rcquii ed equations can be derived from the equation of the conic expressed in the 


r = r F / rj = 0“ / 1 + e £ cos ( 6 + 0 p ) J 


where 


°c = [ » I 2 S I 0( % 2 -1) 2 ) ,/2 

" i r 2 i 

6p = cos [( fij - 1 )/ e l (perigee angle) 
v c = ^/^ 1/2 

*1 and *1 ar<; the normalized (with respect to v q ) radial and horizontal components 
of velocity. 

The above equation can be expressed 


P s / r I " &\ = c 2^°'l sin 0 /^j ' Cj) 


where 


Cj = cos B - l/r 

C 2 - 1 “ cos 0 

P s = semilatus rectum 

For a maneuver that is constrained to be in a horizontal direction, Eq. (1) 
can be solved for ^ as a function of the specified & . 

[ orj sin 0 + (ftj 2 sin 2 0-4c 1 c 2 ) 1 / 2 ]/2c 1 

As there has to be both a positive and negative ^ solution to this equation (one 
trajectory in each rotational direction), the sign choice is resolved in favor of plus 

V 

Tor a maneuver that is applied along the velocity vector, the flight path angle 
Vj is to be held fixed. Using Eq. ( 1 ) 
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tan y ^ - a. ^ /5j = (Cj + c 2 )/ sin 0 

Therefore, 

/?j = [ c 2 /(sin 0 tan y ^ - c ^ )] 

* 

By interchanging the I and F subscripts, Eq. (1) can be expressed 

P s = r F c 2 /(o; F sin 0/0 p - cos 0 + r F /rj.) 

Combining with Eq. (l) using the apogee/perigee constraint &pf Pp = ^ results in 

tan y q = = (1 - 1/ r)/ tan (6/2) 

Inserting into Eq. (1) gives the required horizontal component of velocity for apogee/ 
perigee designation maneuvers. 

Pj = [ r Cg/(r - cos 6) ] 

The derivation of the equation 

6 = cos 1 ( r p cos (e p )/r p ] - 


where 


f e L if g l - n 

e, - 77 if e T > 7T 

_ Lj 1 — / 


for computing the desired central angle 6 between two positions (r p , r p ) which 
satisfies the TPI constraints (e p , Ah) is discussed in Ref. 12. This equation is 
used in the Desired Position Routine. 







l 


3 


5 . DETAILED FLOW DIAGRAMS (!*% 

Figures 6 and 7 contain the detailed flow diagrams' of the main CVbiter rendez- 
vous targeting program and the general maneuver routine, respectively. The follow- 
ing six routines are called by these two programs. 

Iteration Routine 

This routine contains a Newton Raphson iterative driver based on numerically 

computed partials. The routine computes a new estimate of the dependent variable 

x and returns the old values of the error e and x. If the iteration counter c 

exceeds 15, a convergence switch s is set equal to one. 

b conv N 

* 

Coelliptic Mane uver Routine 

This routine computes a coelliptic velocity vector v^ based on a target 
vehicle’s state vector and a delta altitude. 


Phase Match 



) 


This routine phase matches the target state vector to the primary state vec- 
tor. The controlling switch (s ) equals tv o if the leading vehicle leads the 

p n 3l s 6 

other vehicle by more than one revolution: otherwise the switch equals one. If the 
primary vehicle leads to target vehicle, the switch is negative. 

Desired Position Routine 


This routine updates a specified, state vector to the time tp and then offsets 
the updated state vector through either (Ad , Ah) or (e^ , Ah), depending on the 
setting of the switch s, to obtain r^. The routine contains an iterative search to 
solve the (e^. Ah) offset problem, where is defined in Figure 2 and Ah (posi- 
tive when the target orbit is above the primary) is defined as shown below. (This 
represents the TPI geometry used in Apollo and Skylab.) 


TPI point on target orbit 
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Update r.'tiuiine 

This rm . * » dates a state vector based on the update switch 


s update * - 


i Updates through the time Ip * t 
Updates through the time 
Updates through n revolutions 

4 Updates through the angle 8 

5 Updates to where the orbit intersects 
the line defined by r^ 


s 


update 


Search Routine 

This rou » • . , j s the following computations depending on the setting of 

the search swi: m ; 

search 


= n P’inds the time of the n*h apsidal crossing 
(> 0 ) and updates the state vector to that time 

- - 1 Finds the time of the next perigee crossing 
and updates the state vector to that time 


s 


search 


1 . 


~2 


-3 




Updates the state vector through the time 

t - t and computes the elevation angle 
F 

Finds the time associated with a specified 
elevation angle and updates the state vector 
to that time 


The detailed flow charts for these routines are shown in Figures 8 to 13. 
The iterative algorithm used to determine the time associated with the elevation 
angle is described in Ref. 8. 


Each input and output variable in the routine and subroutine call statements 
can be followed by a symbol in brackets. This symbol identifies the notation for 
the corresponding variable in the desired description and flow diagrams of the 
called routine. When identical notation is used, the bracketed symbol is omitted. 
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UNIVERSAL 

CONSTANTS 


INPUT VARIABLES 


£p* £p» £'j ' 1 £■]’ ’ £]) ’ ' ’ 


Av, Ah, A0, n r> 


c» S S S 

‘direct* ene* "man* " o 


— unit 


direct 


— 

< ^4^direct ; 

= 1 

- i = unit (n x r p ) 

jT^° 


i 


+ Av i 


Vpp - [2^1 (r p -f- Ah) / r p 2 (2 + Ah/r p )J 1 / 2 £ 


Call Update Routine 

Input: £ t [£]-!t[ v ]’ £p[£d]' S F 
Output: £rp p [r p j, £ Tp [£ p j 


s ~ 5 

pert * update 


Call Coelliptic Maneuver Routine 

Input: * ^^TF ^ ^ 

Output: Vpp 

Figure 7a. General Maneuver Routine - Detailed Flow Diagram 


(Figure 7c] 
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update^, ^ 



r 

I - 1 r 

0 - 0 - sign 

<r p * r T ) • (r p * Y p ) 

cos 


Call U p date Routine 

Input: r [r], v [v] . t. t , At, n . 


r ' pert* update 


Output: r TF [r F ], S TF [ilp] > t p 




Call Desired Position Routir 


Input: r Tp [r], Y TF [yJ, tpftj. t p . e L> Ah. a£> , s pert , s Ur [s 
Output: r TF [r F ], r TA [r D ] 


_ Cal l Coelliptic Maneuver Routine 

Input: r TF [rI. Y^fv], (r TF - r TA > [Ah ] 



Output: i XA [Y N ] 


Call Update Rout i n e 


Input: r TA [r], v^ A [v], t p [t], At, 


S ,, s 


pert' update 


Output: r TA [r F ], t p 


Call Precision Required Velocity Det ermination Routine {Ref. 1) 

Input: Pp[r 0 ], Sp^], Ita^]. t[t Q ]. t^tj. m. " r [ s^. s #|lg 

Output: Inrll’n] » r«„[r*, 3 . Y_[y-J. r^. i N< s faiV s proj , Ar pr . 


Set Alsrm 

j 

EXIT 

Code a A 
4 
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Figure 7b. General Maneuver Routine - Detailed Flow Diagram 



























Figure 7d. General Maneuver Routine - Detailed Flow Diagram 
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INPUT VA5UA1U.PS 


A 

u 

-■if 

% 

I 

I 



Figure 8. Iteration Routine - Detailed Flow Diagram 



UNIVERSAL 

CONSTANTS INPUT VARIABLES 



Figure 9, Coelliptic Maneuver Routine - Detailed Flow Diagram 
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Figure 10. Phase Match Routine - Detailed Flow Diagram 
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PROGRAM CONSTANTS 

ra 


INPUT VARIABLES 
X. Y. L t p , e L> 

Ah ' ' S pert’ 8 



Call Update Routing 

Input: r, v/t, t p . S pert> S update = 1 

Output: r[r p ], v[v p ] 


C all Update Routine 

Input: r, v, tpft], Afl [ $] # s pert' s update ^ 

Output: , Vp 


£ f = L > £f = - ; - 


c = At = 0, t = t T 


e L>* 


e T = e T - ft 

x~j Lj 


Call Update Routine 

Input: r, y. t, At, s^, S update = 2 

Output: Vp 


Set Alarm Code a. 


c i “ ( r F " *^) cos(ep) / r 




- 

e = cos* (c.) - e L - cos* (r • r p / r r p ) sign ] 

(rp x r) • (r x v) 


~D ~ — F 


Call Iteration Routine 


Input: c, e, At [x] , e^, ^^q[ x q] 
Output: c , At[x], e^, At^fx^], 


OUTPUT VARI ABLES f 

| x F ' V 



Set Alarm 


EXIT 

Code ag 




Figure 11. Desired Position Routine - Detailed Flow Diagram 
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UNIVERSAL 

CONSTANTS 


PROGRAM 

CONSTANTS 


INPUT VARIABLES 



\ 
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0 = sign [(£ S X r) • (£X v) J [cos (r • / r r s ? ' 77 J 


c = c + 1 

0 " fi s ‘ 8 


0=0-2 T7H®C 



At = t „ - t c 


t-J 

1 

M 

< 

1 

ii . 
< 

F S 


|-r - -r - | 


OUTPUT VARIABLES 


I F > v F . At, t F 


Figure 12b. Update Routine - Detailed Flow Diagram 
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Submittal 41: RENDEZVOUS TEPJTEMI PHASE AUTOM ATI C BRAKIN G SEQUENCE, AND TAFXiETII 


INTRODUCTION 


The purpose of the Rendezvous Terminal Phase Braking Program is to 
provide the means of automatically bringing the primary vehicle (Orbiter) within 
desired station-keeping boundaries relative to the target vehicle (or satellite). 

To accomplish this task, the program of necessity contains navigation, targeting 
and guidance functions. 

The program is initiated subsequent to the last midcourse maneuver of 
the rendezvous targeting sequence. Line-of-sigbt corrections, braking correc- 
tions, and filtering of rendezvous measurement sensor data to improve vehicle 
and target state estimates are performed in a sequential manner. At program 
initiation, the relative range is on the order of three to five miles. 

When the primar}' vehicle has achieved a position (and velocity) relative 
to the target which places it within the desired station-keeping boundaries so that 
the station-keeping function can be initiated and maintained, the program is 
terminated. 








N CM EN C I v A T UR E 


Number of biases to be estimated in Unified 


Navigation Filter program 


Measurement code identifying i th measurement 


Thrust of the engine selected for the maneuver; 
used in the Powered Flight Guidance Routines 


Previous range gate passed; subscript used in 
braking (range) gate loop 


Unit vector in direction of relative position vec- 
tor, JL 


Unit vector which defines center of station-keeping 
boundary, relative to target vehicle 


Constant used to determine the range at which 
each range gate search starts when approaching 
that particular range gate 


Constant used to determine how often the line- 


of-sight targeting loop is entered; integer 
number of terminal phase program cycles 


Constant value of range rate added to the mini- 
mum range rate at a given range to insure 
primary vehicle intercept of target vehicle 


Constant used to determine how often the range-rate 
correction targeting loop is entered 


Current estimated primary vehicle mass 


Transformation matrix from reference coordinate 


frame to body axes coordinate frame 
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imt k 1 m r-^v n .'V i aid - u^a l isaSiiaL ■ - *t£k M*i ^ si ^»,ri.- « .tn -1 , ai^Utf - 


? 

¥ 





M ^ y Q2 Transformation matrix from reference coordinate 
frame (in which vehicle states are expressed) to 
LOS coordinate frame axes 


M 


R-M 


Transformation matrix from reference coordinate 
frame to measurement coordinate frame 


M R-SM 

M NB-B 

M NB-M 
M 

J l SM-NB 


Transformation matrix from reference coordinate 
frame to stable member coordinate frame 

Transformation matrix from navigation base frame 
to body axes 

Transformation from navigation base to measurement 
coordinate frame 

Transformation matrix from stable member coordinate 
frame to navigation base 


n 

q i 


Number of discrete braking gates in the range J 
range rate correction schedule 

measured relative parameter at t 

* m 


q PN 


Process noise acceleration 


Local vertical relative position vector (target vehicle 
local vertical) 


-P 


-T 


— ( 1 A > 


B 


Primary vehicle position vector 
Target vehicle position vector 

Aimpoint vector used in Lambert targeting calcula- 
tions 

Switch which controls braking gate targeting cycle 


s 


eng 


Engine select switch 


S freq 


Switch which controls measurement processing 
(navigation) cycle 
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Switch which indicates guidance mode to be used 
*GM * , u • 

in Powered Flight Guidance Routine; "l - two axis 

thrusting; 11 3" modified Delta- v mode; '4 n - 

modified Lambert mode 

s.^ ...Signifies first entry into Unified Navigation Filter program 

s Switch used to select type of targeting scheme 

used in the Terminal Phase Braking'Sequencing 
Program- 


s LOS 


s search 


s 


mk 


s 


Av 


S £v LOS 


Switch which controls line-of-sight targeting cycle 

Indicates target search routine is needed in the 
Unified Navigation Filter program 

Switch which controls navigation cycle 

Switch which indicates if a velocity correction is 
to be made or not 

Switch which indicates LOS correction is to be made 
Current time 


ig 


Maneuver ignition time 


t Measurement time 

Time. associated with primary and target vehicle 
state vectors 


ly Time of bias estimate 

— l Relative velocity vector in target vehicle local 

vertical frame 
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i 

'3 



4 


-P 


Primary vehicle velocity vector 


Target vehicle velocity vector 


W 


Initial filter weighting matrix 


r* . 

1 


Measurement variance used in filter to process 
i th rneasurement data 


Elevation angle of line-of- sight in measurement 
frame 


6 1 


B 


Delta time to ignition for a range-rate correction 
maneuver 


61 


LOS 


Delta time to ignition for a line-of-sight correction 


6t . 


m 


Time between successive measurements within 
the measurement loop 


At 


m 


Basic sequencing cycle time 


a x b 


Velocity change expressed in the body coordinate 
frame 


Av 


LIM 


Magnitude of velocity change below which no maneuver will 
be applied 


a -los 


y 


Velocit}' change expressed in line -of-sight 
coordinate frame 

Value of station-keeping boundary cone angle 
Current estimate of bias 
Gravitational constant of the earth 


v 


Relative velocity vector 
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- UM 


Upp- 1 bound on station -keeping velocity 


Pc ■ * _ ■'. und on station- keeping velocity 


Ar. ■’-. "/‘e, 1 ; ity lower limit below which no line- 
of . '-'. . v<-tion is mode; value to which line- 
of y ;. ! = i* Car velocity is driven if a line-of- 
sltv v.*7 : ■ o - ;on is made 


Angular velocity vector of the line- of- sight 
between the primary and target vehicle 

Magnitude of Wj Qg 

Mab-i- f relative position vector, R 

Range rate between the primary and target 
vehicles 

Relative position vector 

Range of the ith braking gate 

Lov- . • u d on station- keeping position 


max . 
l 


Range rate desired at ith braking gate 

and maximum between braking gates i and i + 1 


P - 
min . 


Minimum range rate desired between 
braking gates i and i + 1 


j?off (LV) Offset aimpoint relative to target point expressed in 
target local vertical frame 

P u Uppe,, , .id on station- keeping position 

6 Azimuth angle of line-of-sight in measure- 

ment Frame 

[ ] Vector expressed in measurement coordinate frame 

( )< 

Prime indicates previous values of a variable, 
e. g. prior measurement parameters, prior 
measurement time, etc. 
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2 . 


FUNCTIONAL FLOW DIAGRAM 


The functional flow diagram for the Rendezvous Terminal Phase 
Braking Program is shown in Figure 1. The program is initialed after the last 
rendezvous midcourse correction maneuver of the rendezvous targeting sequence. 
The relative range between the primary and target vehicle at this point is on the 
order of three to five miles and closing. 

The program sequencing begins with the updating of the estimated 
primary and target vehicle relative state parameters with the appropriate sensor 
data. 

These relative parameters are then used in the Terminal Phase Target- 
ing Program where the necessary calculations are perfo'rmed to see if a line-of- 
sight and/or a braking correction is required to maintain the desired character 
istics of the rendezvous trajectory. The line-of -sight corrections (if performed) 
maintain the intercept by nulling out line-of-sight rates which exceed a desired 
rate. At selected ranges between the primary and target vehicles, braking 
corrections are performed to reduce the closing rate to that specified in the ter- 
minal range/range rate profile, if the closing rate exceeds the desired value. 
During the program sequencing a continuous check is made to insure that the 
closing rate is sufficiently high so that the primary vehicle will intercept the 
target. 

If either a line-of-sight correction and/or range-rate correction is 
necessary, the velocity correction is applied using the appropriate guidance 
mode. 

The program sequencing is then repeated. The program is terminated 
when the desired relative position and velocity conditions are achieved so that the 
station-keeping mode can be initiated and maintained. 
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ENTER 



Exit to 


Station-Keeping Mode 


Figure 1. Rendezvous Terminal Phase Braking Program , 
Functional Flow Diagram 
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3, INPUT AND OUTPUT VARIAB LES 

The Terminal Phase Braking Program consists of three basic functions 
navigation, targeting and guidance. The following is a description of the input 
and output variables for the basic sequencing program, the navigation program 
and the targeting program. The Powered Flight Guidance Program is described 
in Ref. 3. 


Terminal Phase Braking Sequencing Program 


Input Variables 


— P ( t s > ' 
^P (t s } 


Estimated primary vehicle state vector at time t ( 


r T (t ), 

£<* > 


Estimated target vehicle state vector at time t ( 


Number of discrete range gate corrections 


BO ' * 


Range values of the n braking gates 


p 

BO' # ’ 


Range rates desired at the n braking gates 


Switch which controls navigation cycle 

Switch used to select type of targeting scheme 
used in the Terminal Phase Braking Program 

Output Variables 


Ip' Xp 


Primary vehicle state vector for use in station- 
keeping phase 


r T , v ^ 


Target vehicle state vector for use in station- 
keeping phase 


Time tag of above state vectors (can be different 
for active and passive vehicles) 













,-: tc' ,v. \ : 7 . , i 



1 

♦* 


3 


I 


1 



3. 2 Terminal Phase Targeting Routine 


-P 1 -P 


It 


— T 


P 


P 


Input Variables 

Primary vehicle state vector 

Target vehicle state vector 

Relative range between primary and target 
vehicle 

Range rate between primary and target 
vehicle 


i 

”P 


-LOS 


Unit vector in direction of relative range vector 

Angular velocity vector of the line-of-sight be- 
tween the primary and target vehicles 


w 


LOS 


Magnitude of ^ 


M 


R-B 


Matrix transformation between the reference 
coordinate frame and body coordinates 


M 


R-LOS 


Matrix transformation between the reference 
coordinate frame and the line -of-sight coordinate 
frame 


LAM 


Switch used to select type of targeting scheme 
Current time 


— N 


Output Variables 

Unit normal to the trajectory plane (in the 
direction of the angular momentum at ignition 


-lc 
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Offset target position 



Time of upcoming maneuver 




1 . 
ig 

AVg Velocity change of upcoming maneuver in body 

coordinates 


^-LOS Velocity change of upcoming maneuver in line- 

of-sight coordinates 

Av LV Velocity correction in local vertical coordinates 



s 

proj 


S GM 


Engine select switch 

Switch which indicates if velocity correction is to 
be performed during this sequencing of the Ter- 
minal Phase Braking Program 

Switch set when the target vector must be 
projected into the plane defined by i^ T 

Switch which indicates guidance mode to be used 
in the Powered Flight Guidance Sequencing Program 
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4 . 


DESCRIPT ION. Or, ■ W ' AT IONS 
4.1 Terminal PL 0 M.g Seqfuencing Program 

The Terminal Phase Braking Sequencing Program (Figure 4), which is 
the main sequencing program for the terminal phase, is initiated after the last 
midcourse correction in the rendezvous targeting sequence. 

The range/ rar ^ i'Uc terminal braking schedule used in the program 
is determined prior to the initiation of the program and consists of discrete range 
gates and their associate } 1 • • : rd range rates. A minimum range rate is also 
specified throughout the - * • ■-•nul phase to insure primary vehicle intercept of 
the tai gel vehicle. An example of such a braking schedule is shown in Figure 2. 

The sequencing begins with the processing of rendezvous sensor data 
to obtain estimates of range, range rate, line-of -sight rates, etc. These 
estimates are derived from processing the sensor data in the Relative 

State Updating Routine (which is also used throughout the rendezvous sequence. 

Ref. 2) 

These relativ •* ’ r * oter estimates are then used in the Terminal 

Phase Targeting Routine * .> determine if a maneuver (either a braking maneuver, 
line -of-s ight correction or a combination of both) is to be performed. The associ- 
ated maneuver time and guidance parameters are also computed. 

If a maneuver is to be performed, the Powered Flight Guidance 
Sequencing Program (similar to the Servicer Routine in Apollo) is entered with 
the appropriate inputs to accomplish the maneuver.’ 

This basic sequencing is repeated until the primary vehicle is within 
desired station-keeping boundaries relative to the target vehicle (Figure 3). ' 

4. 2 Termin a l Pha .. Tar geting R outine 

The Terminal phase Targeting Routine (Figure 6) computes the neces- 
sary maneuvers to maintain the p r i m a r y vehicle on an intercept with the target 
vehicle while keeping the range/range rate profile within the desired boundaries. 

Two modes of operation are available. The first mode is referred to 
as automatic-line- of-sight control braking and the second automatic Lambert 
braking. 
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Cone angle of station-keeping zone 

p u i Upper and lower values of station- keeping 

boundaries 

W LV > ' Relative offset vector in target vehicle local 
vertical, used to target Lambert braking cor 
rections; primary vehicle will intercept this 
point in the station-keeping zone 


Figure 3. Station-Keeping Boundaries -Station-Keeping Abov, 





M l „ * ,■ ,... ,, ;, 



When s is set to zero, the automatic line-of-sight control braking 

mode is used. If the line-of-sight rate as determined from processing the sensor 
data is above a set limit (typically 0. 1 mr/sec), the line-of-s ight correction neces- 
sary to dri\ e the line-of-sight rate to some level is computed and the appropriate 
ignition time, engine selection and guidance mode switches arc set. Since these 
line-of-sight corrections are made frequently, the maneuver magnitudes are small 
(several feet/second or less) and hence the small RCS thrusters are used to effect 
the maneuver. The maneuver is accomplished by using two-axis thrusting normal 
to the line -of-sight. 

The line-of-sight correction check is typically made every two cycles 
of the main program. (Line-of-sight cycling is determined by k£) 

The range/range rate checks, to insure that the desired terminal profile . 
is being followed, are made after the line-of-sight checks. If the range rate at 
certain pre-selected ranges exceeds the desired range rate a braking maneuver 
is performed to reduce the closing rate. Continuous checks are made to insure 
that the closing rate is above the minimum value to maintain intercept. If it is 
not, then the closing rate is increased. 

The ignition time which is set 6 seconds from the present time allows 
the necessary burn preparations to be made before ignition since these correc- 
tions typically involve significant maneuver sizes. 

The second mode of operation, the automatic Lambert braking, targets 
for an intercept point (either the target vehicle or a point offset from the target 
vehicle indicated by R 0 ^> Figure 3) at each pre-selected braking gate. 

Line of sight rate is implicitly corrected to maintain the intercept trajectory when 
using this mode of operation. 

When the range between the vehicles reaches (1 + k^) times the pre- 
selected range gate, the time of arrival at the range gate is computed. The cal- 
culation assumes the present range-rate remains constant until the range gate 
is reached. The primary and target vehicle state vectors are then advanced to 
this ignition time. 

The time of arrival at the intercept point is redefined by the equation 


(Range at ignition ) 

(Desired range rate at 
this range gate ) 
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Thus t gQ is then used to calculate a new target vector for use in the Lambert 
routine to determine the necessary velocity correction. 

By redefining the intercept point in this manner, the Lambert solution 
forces a reduction in range rate to the desired range rate, insuring intercept in 
a length of time equivalent to the time it would take to travel the present range 
at the constant desired range rate. The line-of-sight rate is automatically cor- 
rected in the Lambert solution to assure intercept. 

The new target vector, time -of-arrival, ignition time and guidance 
mode switches are then used in the Powered Flight Guidance Routines (Ref. 3) 
to effect the maneuver. 

Between braking gates, line of sight corrections are made when necessary 
(as in the fiist mode of operation) to insure arrival at subsequent braking gates and 
to insure intei cept, based on the latest navigated state estimates. (These additional 
line of sight corrections are not normally needed until the last braking gates has 

been passed since the Lambert targeted corrections at each gate are adequate to 
maintain rendezvous intercept. ) 
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5 . 



DI. -! ».£D FLOW DIAGRAMS 



■ ’ on tains detailed flow diagrams of the Terminal Phase 
Braking Sequent.-. ’ J rograi)i, and the Terminal Phase Targeting Routine. 

Each input and output variable in the routine and subroutine call stat- 
ments can be followed by a symbol in brackets. This symbol identifies the nota- 
tion for the corresponding variable in the detailed description and flow diagrams 

e called r, ./hen identical notation is used, the bracketed symbol is 

omitted. 




UN I VERSA L 
CONSTANTS 


PROGRAM 

CONSTANTS 


INPUT VARIABLES 



Figure 4a. 


Terminal Phase Braking Sequencing 
Detailed Flow Diagram 


Program, 
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Call Powered Flight t»uidanco 
Sequencing Program 

Input: t c ,r p (t) r(t) . v p ( t ) |^v ( t ) j 


t . , ITS, s f 2 

ig eng 


S perl J S G1 


/ either Av ^ 
< or t . ft 


r . , s . , \ XT 

— 1c proj - N 


\ depending on 
Output: r p ( t g ) , v p ( t g ) 


Call Precision State and Filter Weighting 
Matrix Extrapolation Routine 

Input: r T (t g ) [r Q ] . v“( l g ‘ ) f v Q 


‘.M* ‘.[‘Fj't-pert]* 

0 tk]. 0 [ Sq ]. 0 [q mag ], 0 [s cont ], 0 t 6 ], 

0W.0[l con 1.0[^ on ].0[T].0[x-], 

0[T’] 

Output: r T (t s ) [£ p ], V T (t s ) [Vp], t g [t p ] 


It * I P | 

— T - Ip 



Figure 4c. Terminal Phase Braking Sequencing Program, 
Detailed Flow Diagram 


S41-20 


























































Figure 5e, 


Terminal Phase Targeting Routine, 
Detailed Flow Diagram 
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6 . 


SUPPLEMENTARY INFORMATION 


The Rendezvous Terminal Phase Braking Targeting and Sequencing Pro- 
grams utilize an inertial state vector formulation of both the target and primary 
vehicle. This formulation is the same as that employed throughout the rendezvous 
phase and therefore the navigation filter used to process the relative measurements 
is the same in both phases. (see Ref. 1). 

The .Targeting Program contains two options; (1) 11 Lambert" maneuvers 
at the discrete braking gates with line of sight, corrections performed as needed 
and (2) corrections down the line of sight at the discrete braking gates with line 
of sight corrections as needed. 

Powered flight guidance studies have shown that for reduced thrust cap- 
ability in the terminal phase powered flight guidance is required. 

Thus, in this situation, option 1 above would be preferred. 

Final studies are presently being done to answer the following questions. 

1. What is the relative range between the target and 
primary vehicle below which it is necessary to 
switch to a relative state formulation of the problem 
as is done in the station-keeping phase? 

2. What modification of the Powered Flight Guidance 
Routine (if any) is necessary to improve per- 
formance in the reduced thrust situation and provide 
the best performance in the nominal thrust case ? 

3. What modification is necessary to a standard 
(nominal) range, range rate braking schedule for 
reduced thrust cases to provide adequate time 

. between range gates for thrusting and navigation 
functions? 
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1. INTRODUCTION 

The purpose of the station-keeping guidance system is to automatically keep 
one orbiting vehicle within a prescribed zone fixed with respect to another oibiting 

vehicle. The active vehicle, i. e. the one performing the station-keeping maneuvers, 

is referred to as the shuttle. The other passive orbiting vehicle is denoted as the 
workshop. The passive vehicle is assumed to be in a low -eccentricity near-earth 
orbit. 

The primary navigation sensor considered is a gimballed tracking radar lo 
cated on board the shuttle. It provides data on relative range and range rate be- 
tween the two vehicles. Also measured are the shaft and trunnion axes gimbal an- 
gles. An inertial measurement unit (IMU) is assumed to be provided on board the 
orbiter. The IMU is used at all times to provide an attitude reference for the ve- 
hicle. The IMU accelerometers are used periodically to monitor the velocity-cor- 
rection burns applied to the shuttle during the station-keeping mode. 

The guidance system presented here is capable of station-keeping the shuttle 
in any arbitrary position with respect to the workshop. This objective is accomp 
lished by periodically applying velocity-correction pulses to the shuttle. These 
velocity corrections are computed by the guidance routine with the objective of 
minimizing the average expenditure of propellant (by the shuttle) pel orbit. 

2, FUNCTIONAL FLOW DIAGRAM 

A functional flow diagram for the station-keeping guidance routine is shown 
in Figure 1. The overall structure of the routine is simple and straight-forward. 
There are two basic subroutines: one is used for computing the normal velocity 
corrections ( s rnode = D and the small midcourse corrections (s mode = 2); the 
other is used for computing boundary-avoidance velocity corrections. The guidance- 
routine call times and mode selection are accomplished by the Station-Keeping Ex-, 
ecutive Routine (Ref. 7). 

Both subroutines use relative position and velocity (shuttle w.r.t. workshop) 
from the Station-Keeping Navigation Routine (Ref. 6) as a basis for computing the 
required velocity corrections. Local-vertical coordinates are used in the normal 
and midcourse -correction modes, workshop fixed coordinates are used in the 
boundary-avoidance mode, In-plane and out-of-plane velocity corrections are com- 
puted separately in the normal and midcourse correction modes. 
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NOMENCLATURE 


Notational Conventions 


* Upper-case letters represent matrices 

® Lower-case and Greek letters reserved for scalars 
and vectors 

^ Vector quantities are underlined, e, g. x 

• Vectors are assumed to be column vectoi'S 

unless explicitly noted 


Dummy 2x2 matrix used in velocity-correction 
computations 

a Elements of A 

B Dummy 2x2 matrix used in velocity-correction computations 

b Elements of B 

C Dummy 2x2 matrix used in velocity-correction computations 

c Elements of C 

d Dummy variable used in velocity-correction com- 

putations 

hg c Height of desired station-keeping limit cycle 

i Unit vector along r (local-vertical coordinates) 

— RL — * 


Symbols 

A 
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.ector alonp; r (workshop-fixed coordinates) 


-NW 


— XWL’ 
-YWL’ 

— ZWL 


i'n ' or normal to station-keeping cone boundary 
( v ‘ ^’lop-fixed coordinates) 


Workshop- fixed frame unit vectors (local-vertical 
r rdinates) 

Midcourse-correction fractions 

Tr. ^formation matrix from local-vertical to workshop- 
fixed coordinates 

Transformation matrix from local- vertical to stable - 
member coordinates 


Dummy test variables used in boundary avoidance 
calculations 

Shi . \e position w.r. t. workshop (stable-member 
coordinates) 


Shuttle position w.r.t. workshop (local-vertical 
coordinates) 

Shuttle position w.r. t. workshop (workshop coordinates) 

Magnitude of component of £ in workshop frame 
X-Z plane 

Lower limit on along workshop Y-axis for which 

boundary-avoidance velocity corrections may be required 

Desired target position for orbiter w, r, t, workshop 
at terminal time t— 
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Switch used to select out-of-plane guidance mode 
Switch set at unity if velocity correction is required 
Switch to select current mode of operation of routine 
Switch set at unity when new tp is required 
Current time 

Terminal time for current guidance limit cycle 

Shuttle velocity w. r. t, workshop (local -vertical 
coordinates) 

Velocity-correction level used for boundary avoidance 

Lowest part on desired limit cycle w, r. t, workshop 
(along vertical axis) 

Maximum desired out-of-plane distance for orbiter 
Minimum desired out-of-plane distance for orbiter 
Dummy variable 

Workshop angular velocity (stable-member coordinates) 
Dummy variable equal to (tp - t) 

Dummy time interval (t n - t^) 

Lower limit on computed velocity-correction magnitude 
Required velocity correction (local-vertical coordinates) 
Station- keeping cone half angle 
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Special Notation 

{ )< A-priori estimated value prior to measurement 

incorporation 

( ) Ensemble average of ( ) 

| ( ) | Magnitude of ( ) 

nr 

( ) Transpose of ( ) 

unit ( ) Unit vector for ( ) 
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4. DESCRIPTION •'•• ■ ■ 'ATI ON’S 
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4. 1 General Infor* 


The station-k ; ■ _ idance routine is capable of maintaining an active ve- 

hicle (shuttle) in a small zone which may be arbitrarily located with respect to a 
passive orbiting vehicle (workshop). The passive vehicle is assumed to be in a iow.- 
eccentricity orbit around *he earth. The station-keeping is accomplished by the 
periodic application ^ umciA velocity-correction pulses. The size and location of 
the station-keeping zone are specified as program constants and input variables 

(e.g. h „ , y . , r v , z . , r ). 

& £c Tmn v - x min — DL 

The guidance iVuiine has three primary modes: (1) normal station-keeping, 

(2) midcourse correction, and (3) boundary avoidance. In the normal mode the 
velocity corrections required to hold the shuttle in the specified zone w. r.t. the 
workshop are computed. In typical situations these corrections are relatively small 
(e.g. 2-5 ft/sec or less). The magnitude and frequency of these corrections is 

dependent on the size and location of the station-keeping zone. The midcourse- 
correction mode uses essentially the same relations as the normal mode. The 
basic idea here is tha< '-s applying small velocity corrections in between the normal 
velocity-correction •:*: • : the total velocity-correction propellant expenditure may 

be reduced. In the he-. . a ry-avoidance mode, special tests are made to see if the 
shuttle is outside of the station-keeping zone and heading away from it. Appropriate 
velocity-correction pulses are applied to the orbiter to return it to the desired zone. 

The times at which each mode of the station-keeping guidance routine is called 
are determined by the Station-Keeping Executive Routine (Ref. 7), 



Three coordinate systems are used in. the station-keeping guidance routine: 

(1) stable-member, (2) local-vertical, and (3) workshop-fixed coordinates. All 
three systems are orthogonal right-handed systems. The relationships between 
these frames are shown in Figure 2. The stable-member system is fixed w. r. t. 
the inertial measure ■ -rtf unit (IMU). The local-vertical system rotates with the 
workshop, as shown in Figure 2, with its. X-axis along the local vertical and its 
Z-axis along the workshop-orbit angular-momentum vector. The normal-mode 
and midcourse-correction computations in the guidance routine are done pri- 
marily in this local-vertical frame. The workshop-fixed frame is fixed w. r.t. the 
desired station-keeping zone. The boundary-avoidance mode computations in the 
guidance routine are performed in this frame. 

4. 2 Normal Guidance Mode 

The analytical development of the basic guidance concept has been extensively 
documented in Refs. 1 and 3 for AAP missions, and in Ref. 2 for SSV missions. 
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P-Frame; inertial reference 
L-Frame: local vertical at workshop 
W-Frame: fixed to workshop 
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Figure 3. Geometry for Guidance Limit Cycles 



Extensive pe A * ■'* ‘nee data is given in these references. It is most convenient to 
consider the in-plane and out-of-plane guidance equations separately. This ap- 
proach will be i ’"v ed here. 

The in- pi: '• problem will be considered first. The basic idea is to put the 
shuttle on a trf : -ry that will terminate at a specified position w. r, t, the work- 

shop (£ dl > at a fixed terminal time { t p ). A typical limit-cycle trajectory is 
shown in Figure 3 for the case where the station-keeping zone is above and in front 
of the workshop. 

The terminal time (tp ) is based on the desired limit-cycle trajectory 
height ( hj c ) and desired minimum altitude of shuttle w. r, t. the workshop 
The basic relaf . Is (Ref. 2): 



where t is the current time, and O)^ is the workshop’s angular velocity. 

The required correction (6v p) that must be made to the current shuttle 

velocity ( v ^ ) in order for the vehicle to arrive at the position ( r^p) a * the 

terminal time t_ is computed in a straight-forward manner (Refs. 1 and 2). 
r 

The basic relatic is 

6 -L = $ RV 1 (t ’ t F^-DL' *RR (t * V* -id'-L 

where r. and v T represent the position and velocity of the shuttle w. r.t. the 
— L» — L» 

workshop, expressed in local-vertical coordinates. 

The matrices and are submatrices of the matrix $ , which is 

KV RrC 

used to extrapolate the shuttle state w. r.t. the workshop forward in time, using 
local-vertical coordinates. The relation is 




%r££ L *R V ( T | 

4 VR * T ^ I 



where t , and t are arbitrary times (t > t . ). The detailed relations for 
n - i n n n-i 

a , $ , $ and . are given in Refs. (1) and (2) as a function of workshop 

RR RV VR VV 

angular velocity (w^) and the time interval from to t n (referred to as t). 



In the detailed flow diagram for the guidance routine (Figure 4) the required 
elements of and ^RR are represented by the matrices A, B, C, and 

the dummy variable d . 

Two out-of-plane guidance modes are provided (Ref, 2). If the desired sta- 
tion-keeping zone is centered in the workshop orbital plane, then Eqs. (2) and (3) 
can be used to compute the required velocity correction ( <5v ). The basic 

equation in this case is simply: 

® V ZL * w W r L,2 cot6 V L, 2 < 4) 

where r^ ^ an ^ ^ are the out-of-plane components of shuttle position and 
velocity w. r. t, the workshop. The quantity o> ^ is the workshop's angular ve- 
locity w. r. t. the earth. The dummy variable 0 is given by: 

es “ W (t F' t) (5) 

where tp is the desired arrival time at the terminal or target point. 

If, on the other hand, it is desired that the station-keeping zone be displaced 
from the workshop orbital plane, then the required velocity correction (Ref. 2) is 
given by 


6v ZL " w W</ Z max ' (z min ) - v L.2 (6) 

The parameters z max a nd z^.^ specify the desired maximum and minimum 
displacements of the shuttle w. r.t, the workshop in the out-of-plane direction. A 
velocity correction is applied only if I r T „ I is less than z . and the relative 
velocity is such as decreases I 2 I further (i. e. v^ ^ r^ ^ is nega- 

tive). 

The boundary-avoidance guidance scheme assumes an inverted truncated 

cone as the desired station-keeping zone. The apex of the cone is at the workshop, 

and the- cone's axis * s assurne d fixed w, r. t. the workshop. The lower 

boundary of the zone is specified by the parameter r , which is its minimum 

min 

distance from the workshop. The size of the station-keeping zone is specified by 
the cone half angle 0 7 . 
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Two boundary-avoidance tests are made. First, if the shuttle is too close to 

the workshop ( r w , < r min > and its velocity is taking it towards the shuttle, then 

a correction is required. The shuttle in this case is given a preselected velocity 

(v ,) away from the workshop. This is accomplished by a velocity correction 
DN 

( 6 v T ) of 


fi - L = (V DN " V W, 1 } - YWL ( ‘ ’ 

where i r T is a unit vector along the workshop-frame Y-axis (station-keeping 
zone con™ is), and L and v w> L are the components of relative position and 

velocity along this axis. 

Next, a test is made to see if the shuttle is inside the desired zone. The lest 
quantity (q) is computed from: 

'> = 1 RW,l‘ COS °Z 181 

A second test is now made to see if the shuttle's velocity is directed away from the 
zone's center line, i. e. the angle between £ w and the cone's axis is increasing. 
The test quantity q is given by: 


q =[v w>1 - <v w 


r* ^RW, 1 ^ 


If both q and q are negative, then the shuttle's component of velocity (w. r. t. 
the workshop) in the direction normal to the station-keeping cone boundary ( j_ 
is given a prespecified value of v p>N' * n i° war< ^ s the cone axis. The 

quired velocity correction to accomplish this is (Ref. 1). 

W = ^ V DN ’ -W * — NW ^ -NW ^ 

where the required velocity correction 6v w is in workshop-fixed coordinates as 
is the relative velocity ( -W*’ 
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DETAILED FLOW DIAGRAMS 


A detailed flow diagram is shown for the Station-Keeping Guidance Routine 
in Figure 4. To operate this routine, navigation information is required from the 
Station-Keeping Navigation Routine, The mode selection and routine-call 

times for both the Station-Keeping Guidance and Navigation Routines are controlled 
by the Station-Keeping Executive Routine, 
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(Figure 4e) 


Figure 4b. Station- Keeping Guidance 
Detailed Flow Diagram 
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Figure 4e. Station-Keeping Guidance, 
Detailed Flow Diagram 
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Su bmi t tal 28: Deorbit Targetin g 
1. INTRODUCTION 

The large entry crossrange capability of the shuttle permits deorbit to a 
specified landing site to be accomplished with a single maneuver. Since the re- 
quired velocity change is smallest when no plane change is made, the equations 
presented here are designed to target the Powered Flight Guidance Routines (Ref- 
erence 3) for an in-plane maneuver. The ignition time for this maneuver is se- 
lected to satisfy entry interface and landing site constraints with minimum fuel 
expenditure. 

If the shuttle had no crossrange capability, then an in-plane deorbit maneu- 
ver to a specified landing site could only occur when that landing site, which ro- 
tates with the earth, intersects the orbital plane of the vehicle. Assuming the 
landing site latitude is less than the orbital inclination angle, and neglecting the 
effects of precession, the landing site will intersect the orbital plane twice every 
twenty-four hours. However, the time difference between these two intersections 
is in general not twelve hours. In the case when the landing site latitude is equal 
to the orbital inclination there will be only one intersection every twenty-four 
hours. 

Since the shuttle has a high crossrange capability, deorbit does not require 
intersection of the landing site vector and the orbital plane. It is possible when- 
ever the angle between the landing site vector and the orbital plane is less than 
approximately 20 deg. In general, there will be two sets of opportunities 
every twenty-four hours. Within each set, there may be several deorbit oppor- 
tunities occurring on consecutive orbits with varying crossrange requirements. 

When the latitude of the landing site approaches the inclination of the orbit, these 
two sets merge to become one. It should be noted, in addition, that if the landing 
site latitude is greater than the orbital inclination, the landing site may still fall 
within the crossrange capability of the vehicle. With these facts in mind, this 
routine has been designed to continue stepping through successive solutions, allow- 
ing the crew to select a particular deorbit opportunity based upon entry crossrange, 
time-to-ignition, required velocity change, landing site lighting conditions, urgency 
of the return, etc.. 

The desired entry range and flight path angle will be considered inputs to this 
routine, since available data relating to footprint size and shape, entry heating at 
various ranges, and optimal entry flight path angle are only preliminary. In future 
revisions, consideration should be given to computing the optimum values of these 
quantities for the particular situation. 




o 


FUNCTION A LF L OW DIAGRAM 


A functional flow diagram presenting the basic approach to the deorbit target- 
ing problem can be found in Figure 3. In addition to the state vector, the primary 
inputs to the routine are the landing site location (latitude and longitude), the entry 
downrange distance, the entry angle (at 400, 000 ft) and the earliest desired time of 
landing. Since the high crossrange capability may make deorbit possible on two or 
more consecutive orbits, after each solution the crew lias the option to recycle 
the program to determine the next possible deorbit opportunity. To give the crew 
the flexibility to evaluate solutions in the future without stepping through all earlier 
opportunities, the earliest desired time-of-landing is included as an input. How- 
ever, the vehicle is assumed to be in coasting flight until the deorbit maneuver, 
and therefore the effects of any maneuvers prior to deorbit are not accounted for. 

After the vehicle stale vector is extrapolated forward to the earliest desired 
time-of-landi ng, the solution process is initiated. This consists of tin ee major 
steps. During the first step the vehicle state is further advanced until the landing 
site, which rotates with the earth, lies sufficiently near the orbital plane so that it 
is within the crossrange (or out-of- plane) capability of the entry phase. During the 
next step an iterative process is used to select the ignition time for this deorbit 
opportunity which requires the smallest velocity change, thus minimizing the fuel 
expenditure. Since the first two stops involve several conic ctpproximations to 
minimize the computer time used, the third step fine tunes the solution by gener- 
ating a precision trajectory which satisfies the constraint on the desired entry angle 
while accounting for gravitational perturbations and the non-impulsive nature of the 
deorbit maneuver. After completion of this step the results are displayed to the 
crew. They may then elect to accept the solution, recycle the routine to solve for 
the next deorbit opportunity, or exit. If they accept the solution, a few minor 
computations are required to initialize the Powered Flight Guidance Routines for a 
modified Lambert aimpoint maneuver. 

To aid the reader in understanding the functional flow diagram, each of the 
three major steps in the solution process is discussed in more detail below. 

2. 1 Determination of the Next Deorbit Opportunity ($tep j) 


To determine the next possible deorbit opportunity, it is necessary to calcu- 
late the inertial location of the landing site (which rotates with the earth) at the 
time-of -landing. Then the angle between the orbital plane and the landing site can 
be used to estimate the crossrange required during entry. To accomplish this, an 
estimate of the time-of-flight difference A t between (1) the interval from de- 
orbit through entry to landing, and (2) the time spent in orbit over the same total 


S28-2 





r^r-js &,. ♦jss*.-: ' ; - ■ : ;' &*J 1 


r. >_£*■:>»«*: f.lfiK^. W-*)-' £ ,WR» ':-•> i 




central angle is used. Analysis has shown that, a constant is probably adequate to 
represent this difference since more precise calculations in the following step will 
compensate for any error. 

Upon completion of the initialization process, the state vector is extra- 
polated forward to the earliest desired time-of-landing. Then the inertial location 
of the landing site at the present state vector time, biased by the time difference 
A t de i is computed. This landing site vector is projected into the orbital plane, 
allowing the in-plane central angle 0j p between the vehicle position and the pro- 
jection of the landing site to be determined. 



The conic routines can now be used to determine the time-of-flight Atj p 
required to coast in orbit through the central angle © ip . If the state was then 
propagated through this central angle, its position vector would be aligned with the 
previously determined projection of the landing site vector. Unfortunately, the 
landing site will move slightly due to earth rotation while the vehicle transfers 
through the central angle. Therefore, the inertial location of the landing site must 
be recomputed, accounting for the time difference At pp explained previously. 
Thus, an iterative process is required to precisely determine the location of the 
landing site at the expected time-of-landing. During the first pass through the de- 
orbit targeting routine, the previously described steps are repeated once to insure 
convergence. However, on subsequent passes no iteration is required, since the 
initial guess achieved by extrapolating the state vector one orbit beyond the previ- 
ous solution guarantees a small value for the time-of-flight correction Atj p . 

Assuming the deorbit maneuver is in-plane, the angle between the orbital 
plane and the landing site location at the estimated time-of-landing can be used to 


measure the crossrange required during the entry phase. If the crossrange is 
within the capability of the vehicle, the solution process continues on to the next 
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step. If not, the vehicle state is extrapolated forward one revolution to the next 
potential deorbit opportunity and the process of estimating the crossrange is re- 
peated. 

It should be noted that the process used to determine the crossrange require- 
ment is only approximate, and therefore a small increment is added to the tolerance 
used in the cross range check to allow for this. A small number of cases which 
pass this check will actually lie outside the vehicle crossrange capability, how- 
ever, a more precise check later will screen these out. 

2, 2 Ignition Time Selection (Step 2) 

During this step in the solution process, an ignition time is selected which 
minimizes the impulsive velocity change required. For these computations the 
projection of the landing 9ite into the orbital plane is assumed to be the real landing 
site. Then, based upon the desired entry downrange distance, a target position at 
entry interface which also lies in the orbital plane can be defined. This target 
position is set 400, 000 ft above the Fischer ellipsoid. 






Using this entry interface target, and the desired entry flight path angle, a 
search is made on the central angle 6 D traversed between the deorbit maneuver and 
entry interface to locate the position and time of the minimum Av maneuver. 
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Then the time-of-flight required for the deorbit and entry phases can be accurately 
determined. Using this time-of-flight, an accurate calculation of the inertial loca- 
tion of the landing site at the time-of-landing can be made, and the entry interface 
target can also be updated. To preserve the central angle of the deorbit phase, the 
impulsive maneuver time is adjusted. Then the ignition time is biased from the 
impulsive time by half the expected length of the maneuver and the state vector is 
extrapolated to this time. 

Since the location of the landing site at the time-of-landing is now known 
accurately, the angle between the orbital plane and the landing site is recomputed 
to precisely measure the entry crossrange required. Then a precision check is 
made, and any solution exceeding the crossrange capability is rejected, thus 
returning the routine to step one to search for the next opportunity. 

2. 3 Precision Solution (Step 3) 

During this step a precision integrated trajectory from deorbit to entry inter- 
face is generated which accounts for both the finite length of the thrusting maneuver 
and the effects of gravitational perturbations. Since the time-of-flight from de- 
orbit to entry interface is known, the Precision Required Velocity Determination 
Routine can be used to generate this trajectory. However, the effects of conic ap- 
proximations in the previous steps and the finite length of the maneuver can cause 
significant error in the reentry angle. Therefore* the resulting entry angle is 
checked and if it is in error, a slight modification is made in the time-of-flight 
from the deorbit maneuver to entry interface to adjust the entry angle. Then the 
precision trajectory is recomputed. After satisfying the flight path angle constraint, 
pertinent data relating to the maneuver can be displayed to the crew or transferred 
to the Mission Planning Module. 
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ENTER 


Input state vector, target location, entry range, entry 
angle, earliest time-of-landing, acceptable crossrange 


Calculate the orbital period 


© Set t 3 equal to the earliest time-of-landing 


INITIALIZATION 

~l ! 


Increment by the 
orbital perioa. 


Calculate inertial land- 
ing site vector at time 


© Extrapolate state vector ahead to t 

<5 

• Use At n p (the approximate time-of-flight 
difference 11 between the sum of the deorbit and entry 
times-of-flight and the time spent in orbit over the 
same central angle) to calculate the inertial landing 
site vector at the estimated time -of- lan ding t 3 + At ^ 


1 

© Calculate in-plane central angle 8 jp betwe.en ve- 
hicle position vector and the projection of the 
landing site vector into the orbital plane 

• Find the conic propagation time Atjp correspond- 
ing to this central angle 




First 

Pass 

INo 


t 3 + Atjp 


STEP 1 


© Recalculate inertial landing site vector at time 


• Estimate entry crossrange by calculating the 
angle between the landing site vector and the 
orbital plane 


'Cross range! 
acceptable 
V ? ^ 


(rough check) 


Figure 3a. Functional Flow Diagram 
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• Calculate the entry interface target position based on the 
projection of the landing site vector into the orbital plane, 
the entry range angle 0p, the radius of the Fischer El- 
lipsoid and the orbital plane 

• Search on 0 D (the central angle traversed between the 
dcorbil maneuver and entry interface) to find the minimumj 
fuel deorbit maneuver 

• Estimate entry time-of-flight based on entry interface 
conditions 

• Update the estimated lime-of -landing based on deorbit 
and entry times -of-flight 

• Recalculate the inertial landing site vector at the time-of- 
landing and its projection into the orbital plane 

• Recalculate the entry interface target position 

• Update the time of the impulsive deorbit maneuver 
to preserve the central angle of the deorbit phase 

• Bias ignition time by half the expected length of the 
maneuver 

• Extrapolate the state vector to the ignition time 

• Recalculate the entry crossrange 


y 


STEP 2 


Use Precision Required Velocity Determination Routine to 
compute a precision traje ctory from deorbit to entry interface] 1 

I 



Adjust entry inter- 
face time-of-arrival 
to give desired y ^ 


STEP 3 


FINAL 

CALCULATIONS 


A 


X 

Initialize Powered 

F° 

j Flight Guidance 


;rew x 

Yes y/ option to 

Continue to nexp 
^opportun- 
sityj/ 

fNo_ 

EXIT 


EXIT 


Figure 3b. Functional Flow Diagram 
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NOMENCLATURE 


a 


a i 
a 2 


3 T 


d 


d ACR 


d DR 


Semimajor axis 

Alarm c ~ih ' ■- taiUre in Av minimization loop 

Alarm c^vhv nature in Precision Required Ve- 
locity DcW n\ motion Routine 

Semimajor axis of Fischer Ellipsoid 

Estimated magnitude of the thrust acceleration 

Semiminor axis of Fischer Ellipsoid 

Number oi columns of navigation filter weighting 
matrix (s -it- 0 In this routine since the matrix 
is not required) 

Maximum acceptable crossrange distance of Orbiter 
Estimated entry crossrange distance 
Entry downrange distance 



L ACS 


f 


QMS 


El 


Magnitude of the engine thrust 

Magnitude of the attitude control system 
translational thrust 

Magnitude of the orbital maneuvering 
system engine thrust 

Entry interface altitude (400,000 ft) 


S28-8 





'L 


— El 


-El 


i 


El, z 


ih 


- LSP 


In 


m 


max 


rev 


P D 

P 


r 

p PFr 


Unit vector formed by the cross product of the 
angular momentum and the landing site vectors 


Unit vector in the direction of r 


El 


First estimate of _igj 


Z-component of the unit vector (z-axis 

assumed North) 


Unit vector in the direction of the angular momentum 


Unit vector in the direction of the landing site 
projection into the orbital plane 


Unit normal to the trajectory plane (in the direction 
of the angular momentum at ignition) 


Sensitivity coefficient used to compute adjustment 
to time-of-ar rival at entry interface 


Estimated vehicle mass 


Iteration counter 


Iteration limit 


Integral number of complete revolutions to be made in 
the transfer (set to zero in this routine) 


Semilatus rectum of deorbit trajectory 


Secant squared' of the desired entry flight path angle 


Secant squared of the offset entry angle used by the 
Powered Flight Guidance Routine 


Precision position vector 
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Entry i '.v'*'^ ■ ■ V ; xlion from Precision Required 
Velocity f;^Tev]Tjv:,-tian Routine 

Position of the impulsive deorbit maneuver 




— El 


- ig 


r 


LS 


— PFT 


Entry intern ~e position 
Position vector at ignition 

Estimate' 'Iso site position at the time of 
landing 

Powered flight offset target vector 


S eng 
s fail 


S FP 
s pert 


S proj 


Engine select switch 

Switch sc J ; nr 1 :''ate non-convergence of Pre- 
cision Re. - >■(•; ri alocity Determination Routine 

Switch set equal to one after the first pass through step one 

Switch indicating which perturbations are to be 
included in the Precision State and Filter Weighting 
Matrix Extrapolation Routine (See Reference 5) 

Switch set when the target vector must be projected 
into the plane defined by _i ^ 





Precision -v: - ector time 

Time of impulsive deorbit maneuver 

Time-of-arrival at entry interface 

Estimated time at which in-orbit position vector is 
coincident with the landing site projection into the 
orbital plane 
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1 LTL 


Desired latest time-of-landing 


Vq velocity vector 

v'A Entry interface velocity from Precision Required 

. 4 J ^termination Routine 

v^ Pre-impulse velocity 


v Entry interface velocity 


-ig 


Ignition velocity vector 


Vpprp associated with the powered flight 

J u et vector 

Vpp Post-impulse radial component of velocity 


— req Required velocity 

v ! r Required velocity on the coasting trajectory 

Post-impulse horizontal component of 

•• ;v * 

* vy 

6t^ Adjustment to At ^ 

Ad^cR Increment added to the acceptable crossrange, 
:ed in rough crossrange check 

IV--. 

Ap Diffwonce between the predicted and desired p 
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Out-of-plane target miss due to the projection 
of the target vector 

Transfer time (t^-t^) 

Transfer time { t ^ “ t ^ ) 

Transfer time ( t ^ - t 2 ) 

Estimated duration of the powered maneuver 

Time-of-flight difference between (1) the interval 
from deorbit through entry to landing, and (2) the 
time spent in orbit over the same total central 
angle 

Time-of-flight required to transfer through the 
central angle 

Required velocity change 

previous value of j Ay j 

Increment in in-plane angle 0^ 

Initial increment in in-plane angle 0 ^ 

Convergence criterion on A 

Convergence criterion on angle A 0 

Post-impulse flight path angle 

Desired entry flight path angle measured 
from the horizontal 

Landing site longitude 

In-plane angle between precision state vector 
and entry interface 
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0 h ''n-plane angle between precision state vector and 

Hnn?-bit position 

0^ :ln **v«ne angle 'over which search is made to 

\* • . minimum deorbit Av 


0 


K 


In-plane central angle traversed during entry 


e 


ip 


Angle between precision state vector and the 
projection of the landing site into the orbital 

r -.V.ne 



e 


p 




LS 


T 


Previous value of & ^ 

Gravitational parameter of the earth (product of 
the earth’s mass and universal gravitation con- 
stant) 

Landing site latitude 
Orbital period 
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4 . DESCRIPTION OF EQUATIONS 

To minimize the size of the Deorbit Targeting Routine, extensive use is 
made of other routines. Therefore, this routine consists primarily of simple equa- 
tions, logical operations, and calls to other routines. Since most of the complicated 
equations requiring detailed explanation are contained in the description of the 
other routines, this section will be limited to a list of items not covered in the text 
describing the functional flow diagram. These items will be listed in their order 
of occurence, and are intended to supplement the detailed flow diagram in sub- 
section 5. 

4 . l Selection of Perturbing Acceleration during 

Precision State Extrapolation 

During the first step in the solution process, which may require long term 

state vector extrapolation, it is desirable to maximize accuracy by including all 

significant perturbing accelerations in the extrapolation process. Therefore, the 

switch s . , which controls the selection of perturbing accelerations in the 
pert 

Precision State Extrapolation Routine, is set to 2. During the later portion of the 
routine, referred to as step three, the switch is reset to 1 , thus limiting the dis- 
turbing acceleration to the J 2 term, the second harmonic of the earth's gravita- 
tional potential function. Since extrapolation during step three is limited to the 
interval from the deorbit maneuver to entry interface, the effects of smaller 
perturbing accelerations are not significant. In addition, extrapolation over this 
interval lies within an iterative loop, and thus may be repeated several times. 

The simplified model can therefore significantly reduce the running time of this 
step. 

4.2 Selection of 0 j p Quadrant 

During the discussion of the functional flow diagram, it was mentioned that 
successive solutions to the deorbit problem (when successive solutions exist) are 
about one revolution apart. To find succeeding solutions to the problem, the state 
vector is extrapolated forward one revolution and then the in-plane central angle 
0 ip between the state vector and the projection of the landing site into the orbital 
plane is computed. Analysis has shown that for some selections of orbital inclina- 
tion and landing site, the correction to the assumption of one revolution may be as 
large as 29°. A lower limit on $ ip of -30° was chosen, thus allowing a small 
margin from the empirically determined limit of -29° . The upper limit on 0j p 
is +330° . Large positive values for 0^ p only occur in situations where no solu- 
tion existed on the previous revolution. 

To determine Gj p , the following equation is used, 

0 IP = cos' 1 [unit (r 0 ) . l’ LSp ] sign [(£ 0 Xi' LS p) • _L h ] 
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where 

Tq = vehicle position vector 

i* * v-'ht vector in the direction of the 
— Lo Jr 

t ; ,‘c u" ing site projection 
i ^ ” ufut angular momentum vector 

This places 0 TP between -180° and +180° and therefore an additional test, shown 

1* Q 0 

in Figure 4b, is mrde 1o f orce B ^ between -30 and +330 . 

In order to make the first entry into step one compatible with subsequent 
entries, the state vector is initially extrapolated forward beyond the earliest desired 
time-of-landing t^' ■ ■ one-twelfth of the orbital period, to the time t^, where 

*3 • *ETL + T ' n 

One -twelfth of the period is nearly equivalent to a central angle of 30° for typical 
(near circular) orbits, and hence makes the first entry into step one compatible 
with later entries, 

4. 3 Effect of Approximate Entry and Deorbit Times- 

of- Flight on ^ntry Crossrange Calculation ’ 

During the (\ st step in the solution process, an estimate of the time of 
landing is necessary to compute the inertial location of the landing site and the 
associated entry crossrange. Since the parameters of the deorbit trajectory have not 
been computed, the deorbit and entry times-of-flight are not known. To estimate 
the landing time, a constant At DE is used to approximately represent the differ- 
ence between the sum of the deorbit and_entry times-of-flight and the time spent 
in orbit over the same total central angle. Preliminary analysis has shown that 
if an average value is selected for this time difference, the maximum error will 
be about 6 minutes. This analysis, described in Reference 7, did not include 
variations in entry time-of-flight for the particular entry range, but further 
analysis is expected to show this effect is small. 

During the first step in the solution, this error will affect the calculation 
of the inertial landing site vector and subsequent entry crossrange computation. 

This effect on the crossrange estimate will be largest for deorbit from a polar 
orbit, and result in a maximum error of less than 90 n.mi. To insure that potenti- 
ally acceptable solutions are not rejected due to errors in the initial crossrange 
estimate, the rough check on crossrange during the first step uses a test criterion 
90 n.mi, larger than the acceptable crossrange Input to the routine. In step two, 
after the time-of-landing has been refined, the crossrange is recomputed and a 
precision check is made. Thus a few cases which pass the first test will be re- 
jected later. 
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Velocity Change Minimization Method 


4. 4 

Step two of the routine includes an iterative search to determine the location 
of the impulsive maneuver which minimizes the velocity change Av, As shown in 
Figure 4d,this iteration uses 0^ , the central angle traversed between the impul- 
sive maneuver and entry interface, as the independent variable. A very simple 
halving step iterator is used to search for the minimum. Although this does not 
converge quickly, it is safe and reliable. The more efficient technique of using a 
slope iteration was not selected because analysis has shown that inflection points 
exist in the relationship of Av and 0^. These inflection points would greatly 
complicate any iteration designed to determine the minimum by driving the slope 
to zero. 



4,5 Required Velocity Equations 

The equations used in the previously described iterative loop to determine 
the required velocity can be found in Reference 2. These equations, shown in 
Figure 4d of the detailed flow diagram, use the initial vehicle position r the 
entry interface position r^, and the desired entry angle Y^j as follows. First 
the tangent of the initial (post-impulse) flight path angle Y^ is computed by 

tan Y ^ = (1 - r D /r EI ) cot ( 0 D /2 ) - r D /r EI tan(V EI > 

where 0^ is the central angle between r ^ and r^j and also the independent 
variable in the search. The semilatus rectum p ^ of the deorbit trajectory can 
then be determined from 

p = 2r D (r D /r EI ~ U 

P ° (r D/ r EI )2 p y - (1 + tanYj 2 ) 



The parameter p y , the secant squared of the desired entry angle, is computed 
once during initialization of the routine. 


The horizontal and radial components of the required velocity are then 
obtained from 


HD 


^P D / r D 


RD 


v HD tan Y l 


The required velocity is then formed and differenced with the premaneuver velocity 
to obtain the impulsive Av. 


v 

— req 


v 

-RD 


unit <r D ) + v HD 


unit 


[*- D x -D* x -d ] 


Av 


V 

— req 


-D 
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4. 6 Entry Time-of-Flight Computation (TBD) 





In Figure 4e of the detailed flow diagram, the time-of-flight from 

entry interface to landing is shown as a function of entry velocity, flight path angle, 
and range. Functionalization of this time -of- flight will be included later when entry 
guidance analysis is complete. 


4.7 In-Pl ane Effect of Approximate Deorbit 



As discussed in subsection 4. 3, the first estimate of the inertial location of 
the landing site is dependent upon an estimate of the time-of-landing. A constant 
time difference At , used to estimate the landing time, may be in error by as 
much as 6 minutes. This led to a significant error in the crossrange estimate for 
a high inclination orbit. For orbits of lower inclination, where the movement of 
the landing site can be nearly parallel to the orbital plane, this same error can 
affect the definition of the entry interface location used in the Av minimization 
iteration. 

The entry interface location, computed early in step two, is based upon the 
projection of the landing site vector into the orbital plane and the desired entry 
range. After the minimization process is complete, the deorbit and entry times- 
of-flight can be accurately calculated. As shown in Figure 4e, another calculation 
of the inertial landing site position is made, thus removing the error due to the 
^ ^ DE a PP I ' ox ^ rna f ion. To maintain the desired entry range input to the routine, 
the. entry interface position is recalculated. This new position will be, at most, 

I. 5 (equivalent to 6 minutes of earth rotation) from the entry interface used in the 
A v minimization. To maintain the geometry of the deorbit phase, the time of the de- 
orbit maneuver is adjusted accordingly so that the central angle from deorbit to 
entry interface is preserved. This adjustment in deorbit time <5t Q1 is computed 
from the following equation 



where j_ f EI is a unit vector in the direction of the entry interface position used 
during minimization, is the new value, is a unit angular momentum vec- 

tor, and t/ 2 it is the inverse of the mean orbital rate. The cross product of the 
unit vectors is nearly equivalent to the angle between them, and the dot product 
gives the proper sign. The mean orbital rate is used to calculate the deorbit time 
adjustment from the angular adjustment. Following this adjustment to the impulsive 
deorbit time, the ignition time for the maneuver is biased from the impulsive time 
by one-half the expected length of the maneuver, thus centering the finite thrust 
maneuver about the impulsive maneuver. 
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*| % 8 Compcnsati Ohiateness and Finite 

Maneuver I . ••• ih 

Step three c " *\ solution process contains calculations which account for 
the finite length of * ^nruufuig maneuver on the required velocity change, and 
compensate for the i PiV.-i^ of the J 2 gravitational perturbation on the deorbit 
trajectory. The Precision Required Velocity Determination Routine is used to 
accomplish these objectives, and the reader should refer to Reference 1 for a 
description of the technique. That routine, however, is designed to maintain the 
terminal (entry interfr ) time-of-arrival, and this can cause changes in the entry 
angle. Preliminary analysis, described in Reference 7, has shown that the nomi- 
nal entry flight path angle error resulting from the oblateness and finite maneuver 
length is about 0. 2°, ’ : can be as large as 0. 6° in extreme cases. Therefore, 

to preserve the desir.*; entry angle, the time-of-arrival at entry interface is ad- 
justed slightly. Delaying the time-of-arrival tends to loft the trajectory and thus 
increase the entry angle. An earlier time-of-arrival will depress the trajectory 
and result in a shallower flight path angle* 

To determine the time-of-arrival adjustment, the approximate sensitivity 
of changes in time-of-flight to changes in entry angle is used. Analysis has shown 
that this sensitivity varies by a factor of about 13, depending on the characteristics 
of the pre-maneuver ■ Tcjectory. However, the sensitivity divided by the deorbit 
time-of-flight varies 1 ' m factor of less than 3. This variation is sufficiently small 
such that a constant can be used as the sensitivity coefficient for all cases. 

To reduce the computations required to constrain entry angle, both here 
and in the Powered Flight Guidance Routines*, the secant squared of the entry angle 
p is used rather than the actual angle. In particular, no inverse* trigonometric 
function evaluations are required. 




-I 
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The sequence of calculations designed to reduce the entry angle error are 
shown in Figures 4f and 4g. First the error Ap^ in the secant squared of the entry 
flight path angle is computed from the following equation; 


A P 


1 - [unit ( r" ) ♦ unit ( 


1 


*2 >r 


where r ” and v '2 are the terminal position and velocity determined by the 
Precision Required Velocity Determination Routine and py is the desired value. 
If the error is too large, the entry interface time-of-arrival t 2 is adjusted as 
follows; 


zi 


1 2 = * 2 ' k y At 12 Ap - 


— — 

The Powered Flight Guidance Routines, described in Reference 3, 
use the same basic technique described here to maintain entry 


angle in the event of off-nominal thrusting conditions. 
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where k ^ is the sensitivity coefficient described earlier and At i2 is the time-of- 
flight from deorbit to entry interface. After adjusting the time -of-arrival, .'the 
Precision Required Velocity Determination Routine is recalled with the adjusted 
time-of-arrival and the results are checked. 

4.9 Offset Entry Angle 

In the process of computing a required velocity, the Precision Required 
Velocity Determination Routine computes an offset target for use during the 
powered flight. For the deorbit maneuver, the powered flight guidance also 
requires an offset entry angle. This offset entry angle, actually the secant squared 
of the angle, is computed from the following equation 


P PFY 


X 


-[ 


1 - unit ( r 


-PFT 


) * unit ( v 


- PFT 


) 


where r ppj is the offset target for the powered flight guidance and XpFT * s 
the associated velocity. 
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Su binittal 25; Powered Flight Guidanc e 
1. INTRODUCTION 

The objective of the Powered Flight Guidance Routines is to issue the proper 
steering and engine cutoff commands such that the desired terminal conditions of the 
maneuver are satisfied. The basic powered flight guidance law used in the orbiter is 
a velocity- to- be- gained concept with cross-product steering. 

The two principle modes of the Powered Flight Guidance Routines are. 

1. Delta-V Maneuver Guidance Mode 

2. Real-Time Required Velocity Up- 
dating Guidance Mode. 

The Delta-V Maneuver Guidance Mode is essentially equivalent to the 
External Delta-V Maneuver Guidance Mode used in APOLLO. The input desired 
velocity change is modified to compensate for the estimated central angle to be 
traversed during the maneuver. Then the object of the powered phase is simply to 
steer the vehicle to achieve this velocity change. 

The Real-Time Required Velocity Updating Mode is a generalized version 
of the Lambert Aim Point Maneuver Mode used in APOLLO. The object of these 
maneuvers is to place the vehicle on a coasting trajectory which will intercept a 
specified target at a specified time. Two new concepts which greatly improve the 
accuracy of these maneuvers are introduced. First, guidance during the maneuver 
is based on a state vector navigated from ignition in a spherical (Keplerian) gravity 
field. Second, the required velocity is not determined using the present vehicle 
position but rather an offset position which accounts for the finite length of the ma- 
neuver. Since this is primarily an equations document, these new concepts are 
treated only briefly in the text. A detailed description and derivation can be found 
in Reference 5. 

Because the calculation of required velocity can be a lengthy process, the 
ability to update the required velocity every major cycle is dependent upon the speed 
of the computer. The APOLLO Guidance Computer required portions of several 
major cycles to complete the solution. The guidance equations described here will 
assume that the orbiter computer will also need portions of several major cycles to 
complete the solution for required velocity. A faster computer would not alter the 
basic concepts presented here, but would simplify the mechanization somewhat. 

The Real-Time Required Velocity Updating Mode may select a specific 
required velocity routine to accomplish one of the following maneuvers. 


A 
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1_ Lambert Aim Point Maneuver 

2. Deorbit Maneuver 

3. Other maneuvers such as a maneuver to an 
orbit with certain specified constraints (TBD). 

The required velocity routines will be subjects of separate documents. 
Since this report is mainly concerned with the documentation of guidance equations, 
logic or computations concerned with monitoring or controlling system operation 
will not be presented. 












2. FUNCTION J. ; - -LOW DIAGRAM 

Powered Flight guidance involves both the prethrust and thrusting phases 
O' f v ^ N.‘' : ; . *'/&’ pier rust computations, shown in Figure 1, are a single step 
F ; pr rto/^-^^-v^!- ninutes prior to the maneuver to prepare the vehicle 

RTNhrdT- paired to process targeting parameters to determine the 
desired vehicle attitude at ignition. In addition, the state vector is advanced to a 
specified time prior to ignition. At this time, an integral number of major cycles 
prior to ignition, the thrusting phase computations, including Powered Flight 
Navigation, are initiated^ Of course, the attitude maneuver necessary to align the 
vehicle to the desired attitude at ignition should be completed before entering the 
t n r u stir r ph :: ; > * «•* . \ X v ' iS . 

The seqeet’, ?. A lunctions performed during the main branch of the powered 
flight phase is illustrated in Figure 2. The guidance computer program known as the 
Servicer Routine, which controls the various subroutines to create a powered flight 
sequence, is not included in this document. The Servicer Routine will call the main 
branch every guidance cycle until engine shutdown has occurred. 

Each guidance cycle begins with the reading of the accelerometers and is 
followed by the updating of the state vector in the Powered Flight Navigation Routine. 
Then the vrJocU -\. v $ + - - ■ 1# ned i s updated in the Cross-Product Steering Routine. 

If steering R, -.v :: vj#. : v * a latter also computes the time-to-go and the steering 
command beginning V' T R*ed time after ignition. 

The targeting calculations used to predict and compensate for gravitational 
perturbations establish an offset target which assumes that the vehicle is under the 
influence of only a spherical gravity field after the expected ignition time. There- 
fore, in the Real-Time Required Velocity Updating Mode, it is necessary to main- 
tain an additional state vector navigated in a spherical gravity field. This dual navi- 
gation should begin at the ignition time assumed in the targeting program if it differs 
from the actual. 

In the Rer v- Required Velocity Mode, another branch of the Powered 

Flight Guidance , u \‘ . nvolving the calculation of required velocity is operated 
independent of the main guidance branch. This separate branch, called the Velocity- 
to-be- Gained Routine, is initiated and controlled by the Servicer Routine and may 
require portion^ of several major guidance cycles to complete its solution. Of 
course, simple velocity-to-be-gained updates computed by decrementing the previ- 
ous value by e V" * velocity change continue in the Cross-Product Steering 

Routine every majt : r le. Normally, the Velocity-to-be-Gained Routine operates 
on a lower priority than the main guidance loop so that the new velocity-to-be-gained 
vector is not used by the Cross-Product Steering Routine until the next guidance cycle. 
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Prethrust Routine 


* Compute Thrust Direction 
Desired at Ignition 

• Advance State Vector to a Specified 
Time Before Ignition 


t 

Attitude Maneuver 


• Orient Vehicle to the Thrust 
Direction Desired at Ignition 


i 

Powered Flight Guidance Routines 
(Fig. 2) 


Figure 1. Powered Flight Program 
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The characteristic of the transfer in the Real-Time Required Velocity 
Updating Mode in relation to the singularity cone of the Lambert problem is deter- 
mined by the targeting program before the powered phase is initiated. This in- 
formation is passed on to this guidance program through the s prQ - switch and is 
used by the Conic Required Velocity Determination Routine to define the transfer 
plane. (See Ref, 3 for a detailed explanation of the singularity cone and Ref. 6 
for the targeting procedure). 

If the s . switch has been set, the transfer will take place in the plane 
proj 

defined by the unit vector J_ N in the direction of the angular momentum vector at 
ignition. If this switch has not been set, there are two possibilities. Under normal 
circumstances the transfer will take place in the plane defined by the vehicle and 
target position vectors. However, unexpected degradation in engine performance 
during flight may prolong the powered maneuver to such an extent that the input 
position vector to the Conic Required Velocity Determination Routine is inside the 
singularity cone. The procedure to cope with this situation is presented below. 

If the s . switch has not been set by the targeting program, the s rone 
proj 

switch, which is an output of the Conic Required Velocity Determination Routine, 
is checked at each guidance cycle. If it is found that this switch has been set, 
indicating that the input position vector is inside the singularity cone, the Servicer 
Routine is directed to bypass the Velocity-to-be-gained Routine for the remainder 
of the powered maneuver. In other words, the remaining powered maneuver will 
be completed simply by decrementing the previous value of the velocity-to-be- 
gained by the sensed velocity change as is done in the Delta- V Mode. 

When the time-to-go becomes less than some predetermined value, 
active steering is suspended and an engine cut-off command is set to be issued at 
the proper time. 




1 
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NOMENCLATURE 


Estimated magnitude of thrust acceleration 


Matrix to rotate the target vector to compensate 
for earth rotation due to change in time of flight 
during deorbit maneuver 


~x'~y 


-TD 


‘ steer 


n min 


Dimension of navigation filter weighting matrix (d = 0 
in this routine since the matrix is not used) 


Thrust 


Magnitude of orbital maneuvering system engine 


thrust 


Magnitude of attitude control system engine trans- 


lational thrust 


Gravity vector in the oblate gravity field 


Gravity vector in the spherical gravity field 


Unit vector in the direction of the angular momentum 
vector normal to the transfer plane 


Unit vectors of local vertical coordinates 


Unit vector of desired thrust direction 


Iteration counter in acceleration computation 


Sensitivity used in computing the desired change 
in flight time to control entry angle during deorbit 


Steering gain 


Intermediate variable in t computation 

go 


Current estimated vehicle mass 


Number of guidance cycles used in thrust 
acceleration magnitude filter 
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3 


rev 


N 


— s 


J^V 


A r 


3 accel 


cone 


eng 


guess 


pert 


proj 


Tnteger number of 360° revolutions used in 
Conic Required Velocity Determination Routine 

Normalized semi-latus rectum of conic 
transfer orbit 

Parameter defining the desired terminal flight 
path angle 

Parameter defining the projected terminal 
flight path angle 

Position vector navigated in the oblate gravity 
field 

Position vector on the coasting trajectory 

Position vector navigated in the spherical 
gravity field 

Offset target vector at tg 
Initial position offset 

Switch indicating whether the acceleration is computed 
from sensed Av's 

( = 0 prethrust estimate^ 

= 1 sensed Av ! s 


) 


Switch in the Conic Required Velocity Determina- 
tion Routine to indicate if the transfer is near 180 c 
(see Ref. 3 for details) 


Engine select switch 

Switch to indicate whether estimate of independent 
variable T will be input to the Conic Required Ve- 
locity Determination Routine (see Ref. 3 for details) 

Switch indicating the perturbing accelerations to be 
included in Precision State and Filter Weighting 
Matrix Extrapolation Routine (see Ref. 2 for details) 

Switch indicating whether the initial and target position 
vectors are to be projected into the plane defined by 
J. ^ (see Ref. 6 for details) 
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^soln 


SwiVh. * vj',. v ting which of two possible solutions 
is i i <’•; It' 1 !, he multi-revolution case (see Ref. 3 
for v’« r - ' ’ 

V* - * 


steer 


Steering enable switch 


/ - 0 inhibit \ 
\ = 1 enable / 


tgo 


Switch indicating whether initial t computation 
har > -e 

V / = 0 t not yet computed \ 

( go 1 

^ = 1 t computed J 

Current state vector time (during thrusting phase, 
this is the time at which the accelerometers are read) 


At 

At 

At 


cut-off 


Guidance cycle time step 
Dumnv V'^nsfer time set to 0 

Value "t- , used to define time to issue engine cut- 

off coi vvn d and terminate active steering 


At ’ Value of t which distinguishes between long or 

enable go 

short maneuver 


At 


At 


to 


t 1 


Time interval before t. to start thrusting phase 
computations 

Time interval prior to t. when initial t^ 0 predic- 
tion is made 


At 


At 


t 2 


tail-off 


Time ' Vr . il after t. when steering is permitted 
* 

Time interval representing the duration of a burn 
at maximum thrust equivalent to the tail-off impulse 
after the engine-off signal is issued 


At 


tail-off, 

OMS 


At, tr of orbital maneuvering system engine 
tail 
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A 1 tail-off, 


At . rr of attitude control system engine for 
tail-off 

translational m an e u v e r 


Change in time of arrival required to satisfy ter- 
minal flight path angle in a deorbit maneuver 


Time of arrival at r(t 2 ) 

Time-to-go before engine cut-off 
Nominal engine ignition time 

Velocity vector navigated in the oblate gravity field 
Velocity vector navigated in the spherical gravity 


v 

— req 


— req 


A^LV 


Velocity-to-be-gained vector 
Magnitude of 

Required velocity vector at the offset initial posi- 
tion (defines the coasting trajectory) 

Required velocity at current position (no initial 
position offset) 

Measured velocity increment vector due to 
thrust in one guidance cycle 

Magnitude of Av 

kth value of sensed Av saved for acceleration 
computation 

Desired velocity change vector input to Delta- V 
Guidance Mode 

Minimum sensed Av which will allow acceleration 
filter computations to be made 


Av k for n min cycleS 
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Av^,- In- plane components of Av^y 

■ v mpensated in- plane components of Av 



i umponents of Av ^y 
Exhaust velocity 

Exhaust velocity of -the orbital maneuvering 
system engine 

"xhaust velocity of the attitude control system 
engine for translational maneuver 

Reciprocal of normalized semi- major axis of 
conic transfer orbit 

Tolerance criterion establishing a cone 
around the negative target position di- 
rection inside of which the Conic Required 
Velocity Determination Routine will define 
the transfer plane by j _ ^ 

Projected terminal flight path angle with respect 
to local horizontal (negative downward) 

Converged value of iteration variable used in Conic 
Required Velocity Determination Routine 

Previous value of T 
Estimated value of t 
Time rate of change of T 
Ratio of | r { t 2 ) J to |r'(t)| 

Estimated central angle traversed during thrusting maneuver 
in Delta- V Guidance Mode 

Earth’s gravitational constant 
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T 


Time associated with current required velocity 
Tp Previous value of T 

^ Angular velocity command 

W earth Magnitude of the earth's angular velocity 
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Prethrust 

Maneuver to ignition attitude 

Initialize powered flight navigation 
Start powered flight routines (including v 

6 

routine in the Real-Time Required Velocity 
Mode) 

Initial t prediction 
Ignition 

Start dual state vector navigation if 
in the Real-Time Required Velocity Mode 

Short maneuver engine cutoff (t < ^ t ena b^ e ^ 
Enable steering 


Terminate steering and issue cutoff command 

( t < At , », ) 

v go cut-off 

Engine cutoff 


Terminate powered flight routines 


Figure 5, Sequence of Events 


















Call Conic Required Velocity 
Determination Routine (Ref. 3) 


Input: r(t. g )[ r Q ] , r(t 2 )[ r j] , 

(t 0 - t. )[ At] , n 
2 lg 1 J rev 

s s r 

soln* guess* guess ’ 
C 0 € cone ^ ' S proj’ — N 
Output: v req (t ig )[ v 0 ], 

r guess ^ r 0^ 


V (t. ) - V (t. ) - v(t. ) 

-g lg -req' lg - lg 


■ - 


| V (t. ) I V (t. ) 

*-£ 1£ ‘ “g lg 


r(t. ) x v(t. J Av 


t * I (t ig )] a T 


Av = Av 
— c 1 — xz 


unit (Av xz ) < 


+ unit (Av x i ) sin 
— xz y 

(t. ) - Av _ + Av.. i 


I ’(t ig ) - r(t ig ) + Ar(t. g ) 


s =1 

guess 


Call Conic Required Velocity 
Determination Routine (Ref. 3) 

Input: 

I'VIZO 1 ’ i (t 2 )[ ^l ] * 


(vV^'i'Vv 


s s r 

soln* guess* guess 


€ 0G ^ € cone ^ * S proj' — N 

Output: 

v f (t. ) [ v n ] , r [ r n ] 

— req lg 1 — 0 J L 0 


v (t. ) = v* (t. ) - v(t. ) 
-g ig -req lg - ig 


Figure 6b. Prethrust Phase, 

Detailed Flow Diagram 
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UNIVERSAL 

CONSTANT 



PROGRAM CONSTANTS 


At, At . : , n . , Av 

I skip min min 



INPUT VARIABLES 


t, r(t), v(t). t. , v (t), 
- iR ~g 

a T ‘ At tail-off' v exh 



Call Powered Flight Navigation Routine (Ref. 1) 


Input: r(t), v(t), At'[ At] ( Av [ Av 

Output: g(t)[ jg(t + At)] 


sensed^ 'i*^ ,s pert 


Enter every major 
cycle after initial 
cycle 


Read Accelerometers 
Output: Av 


First 

Pass 


Exit 

Return to ® at start 
of next major guidance 
cycle 


s = 1 
pert 


Figure 7a, Powered Flight Routines, 
Detailed Flow Diagram 
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Call Cross-Product Steering Routine (Fig. 9) 


Input: v (t - At). Av, t.^, t, r(t), v(t), a T 


^tail-off* v exh ' ^tgo 


Output: v (t) 

-g 



Call Powered Flight Navigation Routine (Ref. 1) 


Input: r (t - At)[ r(t)] , v (t - At)[ v(t)] , At, 

*1 [Av sensed ].fi s (t- At)[ £ (t)]. S pert 
Output: r (t) [ r(t + At)] f v (t)[ v(t + At) ] , 

£g(t)[ £(t + At)] 


l OUTPUT VARIA BLES (Req. - Vel. Up- 

dating Mode only) 


t. £ r W < Z s (t) ' £ s (t) 1 -g (t) ' 


a T ’ s accel 


c ' 7 Xi c 

r (t), v(t) 


EXIT 


(Return to @ at start of next 
major guidance cycle) 


Figure 7c. Powered Flight Routines, 
Detailed Flow Diagram 
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PROGRAM CONSTANTS 


INPUT VARIABLES 


u, At u , At t2 , k steer , 
^enable ’ ^cut-off 


v (t - At), Av , t. , l 
-g - ig 

r(t), v(t). a T . At tail _ o 


V u . S. 

exh tgo 


v (t) = v (t - At) - Av 
-g -g 


* <i i g - ",1 


s 4 = 1 


s tgc 

= 1 

> 

J 

■ 

Compensate v 
if required. 

r (t) for ullage 


i 

t = - 

go 

1 - fi (o 1 
a rj. 

! 

V 


t > At i i 
.go enable^ 


Compute t for short 
maneuver ^°(TBD) 


s =1 

steer 


s , - 0 

steer 


Figure 8a. Cross-Product Steering Routine, 
Detailed Flow Diagram 
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Figure 9a. Velocity-to-be-Gained Routine, 
Detailed Flow Diagram 
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Figure 9c. Velocity-to-be-Gained Routine, 
Detailed Flow Diagram 
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1. INTROD p/ i 

The Entry-Guidance Routine presented here is designed to take the orbiter ve- 
hicle from ert : -T ? ^ (h » 400,000 ft) through the critical heating phase of 

entry down to tt r- n' \ri of the approach phase (h 100, 000 ft). The basic ideas are 
outlined in Re' X" v J : Y). Simulation results demonstrating the feasibility of the 

concept are given in Ref, (2). 

There are three basic guidance modes: 


( 1 .) 


(2. ) 


programmed-maneuver mode in which the 
vehicle is oriented with a zero roll angle (wings up), 
and an angle--of~attack corresponding to maximum 

v T - iV,. t the point of pullup, 

*- 

j\ constant heating -rate mode during which the stag- 
nation-point heating rate is held constant at a pre- 
selected value, chosen essentially to minimize heat 
loads on the vehicle without violating maximum 
temperature constraints. 


(3, ) A reference trajectory mode during which the vehicle 
follows a prestored stored trajectory designed to get 
* t i; ■ le to the terminal point with a minimum re~ 

' ’ : - e ' YPS weight, and without violating operational 
constraints on the vehicle. 

Thermal control is provided by varying the magnitude of the roll angle so as 
to follow a density-vs. -speed profile. Density information is derived from IMU 
measurements of the aerodynamic specific force acting on the vehicle. A-priori 
knowledge of the vehicled mass, effective aerodynamic area, and drag coefficient 
(Cp) are required in the process. 


Range control is provided by changing the angle-of- attack of the vehicle. 
Upper and lower limits on angle-of-attack are required in order not to violate 
operational con< v v^ ; s on the vehicle. Lateral trajectory control is obtained by 
reversing the d rection of the roll angle. 
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2 * FUNCTIONAL FLOW DIAGRAM 

The basic information flow in the Entry Guidance Routine is shown in Figure 1. 
This is based on the -guidance concept of Ref. (2. ). 

After the routine is entered, a series of targeting computations are made. 

This involves the computation of quantities such as the current vehicle heading (\b), 
the desired great-circle heading to the target point ($ D ), range to the target point 
(0 ), cross-track distance to the target point ^CT*’ and down-range distance to 
the target point (0 ^ ). 

The particular guidance mode to be entered is next determined. There are 
three possible guidance modes: 



(1. ) Initial programmed maneuver 

(2.) Constant stagnation-point heating rate guidance 

(3. ) Stored reference-trajectory guidance 


The constant heating -rate mode is entered when the vehicle's vertical velocity is 
greater (more positive) than a preselected value. The reference-trajectory mode 
is entered when the magnitude of the vehicle's relative velocity is less than a pre- 
selected value. 


In the programmed maneuver mode the vehicle is oriented with a zero roll 
angle (wings up) and an angle -of- attack corresponding to the maximum aerodynamic 
lift coefficient. This orientation is maintained until the heating-rate mode is entered. 

In the constant heating-rate mode, the density altitude (h^ ) required to attain 
the desired stagnation point heating rate (q^) is first computed. An angle-of-attack 
command (c*q ) is next computed, based on the desired range-to-go for the heating- 
rate mode. Finally, roll-angle magnitude commands ( 0 ^-) are computed to control 
the vertical-plane motion of the vehicle so as to follow a density- altitude vs. speed 
profile. No roll reversals take place in this mode. 

In the reference-trajectory mode, the required reference-trajectory quantities 
are first obtained from the stored table at the current speed. These include angle- 
of-attack (» D )> altitude (h D ), range-to-go (r GD )> and the ratio of cross-track to 
down-range distance-to-go ( 77 ). The angle of attack command is then computed 

as a perturbation from the reference value (or^) based on the difference between 
the stored and measured values of range-to-go. Roll-angle magnitude commands 
are computed in the same manner as for the constant heating-rate mode, except 
that the desired density altitude (h-^ ) is from the stored table. Roll-angle direction 
is based on a comparison between the current estimate of t) and the reference-tra- 
jectory value. 
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NOMENCLATURE 
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Notational Conventions 

Upper-case letters represent matrices 

Lower-case and Greek letters reserve for scalars and vectors 
Vector quantities are underlined, e, g. x 

Vectors are assumed to be column vectors unless explicitly noted 

Symbols 

a Effective aerodynamic area for vehicle 


C 0 ,C 1 Coefficients used to compute a q from c 


'D 


D. 


D 


MIN 


~D 


MAX 


— N 


h 

h 


D 


Aerodynamic drag coefficient for vehicle 

Coefficient in desired-density relation for constant heating- 
rate mode 

Coefficient used in relation for desired c^ 

Desired value of c^ for constant heating-rate mode 
Lowest permissible value of c^ for constant heating-rate mode 
Highest permissible value of c^ for constant heating-rate mode 
Drag force per unit mass 

Aerodynamic force per unit mass on vehicle 
Stored array of reference trajectory altitudes 
Vehicle altitude above Fischer ellipsoid 
Density altitude 
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F prog 

>V:Y 



hr 


h 0 ,h 

p p 


iRE 


— pole 

— EE 


-SE 
— GCE 


— RDE 
-VR 


-iXA 


-YA 


-HOR 

k 


Verticf*.., ,:.ocity at which programmed mode is terminated. 

v£V^' : v^<-- iCtiS ; c; altitude and derivative w* r. t„ time 

• * * ■ . j - " 

- „ . , _ 

* * . ^ “*• a. ■. * -- 

. -^w-v to vertical velocity and derivatives w. r.t. time 

Scale height for exponential density-variation model 
Filter ■' al velocity 

. Altit^rk* i .which p would occur and its derivative w. r. t. time 
Unit vector along r ^ (earth-fixed coordinates) 

Unit vector along North pole (earth-fixed coordinates) 

Unit vector directed towards local East (earth-fixed coordinates) 

Unit ve^or directed towards local South (earth-fixed coordinates) 

: .’\t . 

f , : vunt .-* >ng desired great-circle heading direction (earth- 

fixed " > 'V 1 nates) 

Unit vector along desired terminal position (earth-fixed coordinates) 

Unit vector along v^ (stable-member coordinates) 

Unit vector along vehicle longitudinal axis (stable member 
coordinates) 

.Unit v£rh:;V ?ong vehicle lateral axis (stable-member coordinates) 

i : 

Unit vector normal to plane of vehicle's position and relative 
velocity vectors 

U -variable used for reference -trajectory lookup 
V 
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Sensitivity factor in angle -of- attack relation 

Value used for k in thermal control portion of ref. traj. mode 
a 

Value used for in final position control portion of 
ref. traj, mode 

Coefficients used in roll -command relation 

Fraction of T? at which roll angle should be reversed 

Desired latitude at the end of entry 

Desired longitude at the end of entry 

Transformation matrix from stable member to earth- 
fixed coordinates 

Mass of vehicle 

Index for computation- cycle time 
Mach number 

Vehicle position (stable -member coordinates) 

Target-point position vector (earth-fixed coordinates) 

Earth radius (nominal) 

Vehicle position (earth-fixed coordinates) 

Range to go to target point 

Desired value of r^ at current v^ (from stored trajectory) 

Desired range to be covered in the constant heating- 
rate mode 
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r r 

REF 

U APP 

r r 

°MIN 

r G 


— GN 


s 


0 


s 


q 


S REF 


s 


v 



- t 


V 


a 



Xn 


'H 


Nominal range to be covered in the reference trajectory mode 
Nominal range to be covered in approach phase 
Lower limit for r^ in angle-of-attack computation 

G q 

Range-to-go used in reference-tbajectory mode guidance 
computations 

Stored array of reference-trajectory range-to-go 
Dummy variable used in roll-reversal logic 
Switch used to start constant heating rate mode 
Switch used to start reference-trajectory mode 
Vertical component of specific force on vehicle 
Desired vertical component of specific force 
Current time 

Vehicle velocity (absolute in stable-member coordinates) 
Speed factor used in angle-of-attack command relation 

Relative velocity at which constant heating-rate mode 
is terminated 

Reference-trajectory array of vehicle speed w.r.t. 
air mass (22 elements) 

Horizontal component of vehicle ! s velocity (absolute) 
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-R 
— RE 
-RITE 


Vehicle velocity w. r.t. air mass (stable-member coordinates) 
Vehicle velocity w. r. t. air mass (earth-fixed coordinates) 
Horizontal component of 



v RLO 

"c 

MAX 
a CMAX 

a CMIN 

a D 

"max. 

q 

a MIN • 

q 

— N 

0 T 

Ail) 

A *0 
A v 

At 


Lower limit on v 
a 

Angle-of-attack command 

Angle -of- attack corresponding to maximum c^ 

Maximum permissible value of a q 
Minimum permissible value of a ^ 

Desired value of a at current (from stored trajectory) 

Maximum permissible angle-of-attack in constant heating- 
rate mode 

Smallest permissible angle of attack in constant heating- 
rate mode 

Stored array of reference-trajectory angle of attack 
Roll angle command 
Computed vehicle roll angle 

Different between current and desired heading of vehicle 
w, r. t. air mass . . 

Value of A<) on first pass 

Accelerometer-measured velocity change from previous to 
. present computation-cycle time 

Time interval from previous to present computation- cycle time 
Local heading (w. r.t. South) 
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^ D 
ClJ^ 


0 MAX' 

>V 

0 MIN 

0 MAX’ 
0 

MIN 


0 

y OLD 

e 

®CT 
6 DR 

i n 


”d 


J2n 

Po 

p 


p 

£ 


L * . ,’ed heading (w. r* t. South) 

- rotor in vertical velocity filter 

i, ' * 

9 - *■ 

-* r - « f, r" 

Levels used in ^ computations 

Levels used in 0^ computations 
V 

Dun any variable used in roll-angle computations 
Previous value of 0 q 

Great-circle angle from the current position to the desired 
target point 

Cr* ♦ > • track component of 0 
Di:.»r- range component of 0 
Ratio of 0£irp ®DR 

Desired value of rj at current v-^ (from stored trajectory) 
Stored array of reference-trajectory r\ 

Sea-level value of earth 1 s density 

Estimated density from specific force measurements 
Desired density for constant heating-rate mode 
Earth’s gravitational constant 

Dummy variable used in reference trajectory lookup 
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Spe cial Notation 

( )» A-priori estimated value prior to measure- 

ment incorporation 


( ) 

Ensemble average of ( ) 

l< )l 

Magnitude of ( ) 

( ) T 

Transpose of ( ) 

unit ( ) 

Unit vector for ( ) 

sign ( ) 

Algebraic sign associated with ( ) 


+ 1 or - 1 , with s ign (0) ^ +1 



aMadigii »**» «*«<* 





r 



v 

A v 
At 



Vehicle position (stable-member coordinates) 

Vehicle velocity w. r.t. air mass (stable 
member coordinates) 

Vehicle velocity (stable-member coords. ) 

IMU-measurement velocity change 

Time interval over which Av is taken 
Output Variables 
Angle-of-attack command 

Roll angle command 
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Submittal 59 


CLOSED FORM ENTRY GUIDANCE LOGIC FOR 
THE HIGH CROSS-RANGE ORBITER 


1.0 SUMMARY 


Entry guidance logic has been developed for the space shuttle which 
controls the entry trajectory by roll modulation while using a preselected 
angle of attack profile, which is a function of velocity. Range predictions 
are based upon an analytic solution to the equation of motion for equilibrium 
glide and constant load factor profiles. Inplane range errors are nulled by 
changing the magnitude of the roll angle and cross-range errors are nulled 
by roll reversals. 

The basic guidance concept consists of three phases: a constant 

temperature phase, an equilibrium glide phase, and a constant load I actor 
phase. The constant temperature phase is entered first and is designed 
to control the trajectory to a constant temperature profile until an 
inertial velocity of 25 000 fps is reached. At this point in the trajectory, 
the initial descent rate has been controlled and near equilibrium flight 
conditions exist. At this point, the equilibrium glide phase is entered and 
entry range predictions are intiated. These range predictions are based on 
an equilibrium glide trajectory until a load factor of 1.5g is reached, 
followed by a constant load factor trajectory of 1.5g until transition. 

The roll angle during the equilibrium glide phase is selected to 
null the inplane range errors. When the resultant equilibrium glide 
trajectory intersects the constant g trajectory required to reach the 
target, control is transferred from the equilibrium glide phase to the 
constant g phase. At Mach 6, the entry guidance is terminated and control 
is transferred to the transition guidance* 
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2.0 INTRODU CTION 


Analysis of entry trajectory shaping studies of the high cross-range 
orbiter has resulted in an understanding of the relationship between tra- 
jectory shaping and entry constraints and objectives (such as temperature 
limits, minimum. TPS weight requirements, and load factor constraints). 

This analysis indicated that all know orbiter constraints and objectives 
could be met through proper entry targeting, and therefore, direct con- 
trol of the trajectory to minimize constraint parameters is not necessary. 
This analysis also indicated that ranging could be accomplished early in 
the entry with negligible effect on the trajectory shape. In fact, indi- 
cations are that delaying ranging until after the major aerodynamic heating 
has been passed could cause an impact on other constraints, such as load 
factor, later in the entry. 

The analysis further indicated that several simple control modes can 
be used to satisfactorily control the orbiter trajectory. Analysis of 
these modes indicated that a combination equilibrium glide and constant g 
mode will not only produce a satisfactory trajectory but can also be used 
as a basis for closed-form guidance logic. This document presents an 
analytical guidance technique based on this concept. Roll angle is used 
to control inplane ranging and roll reversals are used to control cross 
range. The angle of attack profiles are predefined functions of velocity. 
Section 4 discusses the guidance concept and subsequent sections present 
a description of the guidance logic. Equation derivations, guidance flow 
charts, and a detailed description of the guidance logic are presented in 
the appendixes . 


3.0 SYMBOLS 


A 


vector pointing at target 


A T 

An 

Arz 


A projected into EFT 


components of A^ 
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ATM* 


c * 

FNM 


C * 
rnin 


eg . 


D * 
ref 


H * 
up 



yation matrix from EFF to EFT frame 

conversi^.i from feet to n. mi. 

conversion from radians to n. mi. 

cos sine 0 
drag 

dt„ b \ h V? .eg level for constant g range prediction 

Ji - *■ ’ ■ 

drag required to reach target 

drag reference 

gravity acceleration of earth 
drag limit 
altitude 
altitvvv. rate 

afiiiohi'^-.rlc density altitude constant 
11 ft* sector orientation flag ro preentry 


* 

These f 



v 


ab*-;ls appear in the guidance flow charts in appendix D 



1 

1 



1 


I 


i 


•I 


1 



V 






m 


■m 


■ ’&SH4 r **insej 








4 


IFT* 

IG* 

ISTP* 

K10* 

K2R0L* 

L 

l/d 

L 

DyCommand 

L 

V ef 

LATSW* 

LMN* 

LOD* 

LODI* 

m 

N 


flar for first pass through range prediction 

flag to transfer to constant g phase 

flag to determine sequence in range prediction 

constant in D „ equation in constant heat rate 
rei 

roll direction indicator 

lift force magnitude in the vertical plane 
lift to drag ratio of the orbiter 

commanded L/D in the vertical plane 

reference L/D in the vertical plane 

flag to inhibit roll reversals through l80° 

L/D command for 5° deviation from lift vector up for 
cross-range control 

vehicle L/D 

desired inplane L/D 

vehicle mass 

total load factor 



R* 


stagnation point heat rate 
commanded Q 

radius vector 




J 


*These symbols appear in the guidance flow charts in appendix D 
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R* 


R * 
CG 


REQ* 

RPT* 

RT* 

RTD* 

RTG 


R * 
ref 


Select* 

TPS 

UR 

UT 

V 

V 


V * 
CG 


V * 
E 


V 


El 


V 


V * 

Q 


v * 

V Q2 


altitude rate 

predicted constant g range 

predicted equilibrium glide. range 
predicted transition range 
total range to target 
conversion from radians to degrees 
total range to transition point 
reference R 

flag to determine guidance phase 
thermal protection system 
unit position vector 
unit target vector 
velocity 

time derivative of velocity 

inertial velocity to enter constant g phase 

relative velocity 

inertial velocity at entry interface 
inertial velocity 

relative velocity to start transition 
inertial velocity to start transition 


These symbols appear in the guidance flow charts in appendix D 
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local satellite velocity 


local satellite velocity 


V * 
XX 


velocity to start range prediction 
velocity to start transition a modulation 


earth fixed frame (EFF) 


earth fixed topocentric frame (EFT) 


earth rotation rate 


vehicle weight 

lateral deadband switch point 

angle of attack 

angle of attack command 

flight-path angle 
time derivative of y 

inertial flight-path angle at entry interface 


central angle to target 


These symbols appear in the guidance flow charts in appendix D 
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0 


0 * 


\|i* 

^ * 


05GSW* 


Pi 

density 

density at sea level 

density at sea level 

roll angle 

roll angle command 

relative azimuth 
relative azimuth to target 

flag to begin guidance 



*These symbols appear in the guidance flow charts in appendix D 
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4.0 GUIDANCE CONCEPT 


* 



1 

l 

•J 

■l 

i 

d 

-j 


The entry guidance must keep pe&k acceleration levels, maximum 
temperatures, and heat loads within limits while maintaining ranging 
capability. The guidance must operate over a wide. range of initial con- 
ditions and vehicle lift to drag ratios with a minimum of changes to the 
guidance software. The guidance must also be insensitive to navigation 
system errors. One means of accomplishing this is to develop a set of 
analytic trajectory prediction equations based on a flight profile that 
satisfies the objectives previously mentioned. Trajectory shaping 
studies showed that two control modes can be combined to satisfy the tra- 
jectory limits and objectives, and would also be amenable to analytic solu- 
tions of trajectory parameters for constant and near optimum angle of 
attack profiles. These modes are equilibrium glide and constant g. This 
document presents the guidance logic for both a constant and a near 
optimum variable angle of attack profile. A detailed description of the 
guidance logic can be found in section 5-0, however, a brief overview of 
the guidance concept follows. 

From 0.05g to an inertial velocity of 25 000 fps, the guidance con- 
trols the trajectory to a constant temperature profile. This profile 
controls the initial descent rate and stabilizes the trajectory prior to 
initiating ranging at an inertial velocity of 25 000 fps. Between an 
inertial velocity of 25 000 fps and a load factor of 1.5g, the entry 
trajectory is controlled to an equilibrium glide flight mode. During 
this phase the roll angle for equilibrium flight is analytically computed 
to satisfy the entry ranging requirements. The resultant equilibrium 
glide trajectory is maintained from the point in the trajectory where 
the equilibrium glide drag level is greater than the constant heat rate 
drag level (point 1 in fig. l) to the point in the trajectory where the 
constant drag level required to reach the target is equal to the drag 
level resulting from the equilibrium glide trajectory (point 2 in fig. l). 
From this point until transition, the guidance commands the roll angle 
required to maintain the constant g level required to reach the target. 

At Mach 6, the guidance transfers to the transition guidance mode. 
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Figure 1 Guidance concept for the high cross-range orblter. 


5.0 GUIDANCE LOGIC DESCRIPTION 


The basic guidance logic must perform three primary functions, 
these are trajectory parameter prediction, targeting, and attitude 
command generation. The guidance first performs trajectory and range 
predictions and then the controller converts these data into attitude 
commands which are provided to the autopilot for execution. An analytic 
reference trajectory is recomputed each computer cycle to correct for 
range errors. Based on this recomputed reference trajectory, a refer- 
ence lift to drag ratio (L/D), drag level, and altitude rate are 
analytically computed and provided to the controller. 

The total guidance logic can be divided into four major phases as 
depicted in figure 2. These phases are preentry, constant heat rate, 
equilibrium glide ranging, and constant g ranging. Several service 
routines are used during each phase such as targeting, lateral logic, 
roll command, and controller. The major phases are described in 
sections 5.1 through 5*^» and the service routines are described in 
section 5-5. A complete derivation of the range prediction equations 
and math flow is presented in the appendixes. 
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, 5.1 Preentry Phase 

The primary activity of the preentry phase is the computation of 
the attitude ho < ,-j;< ... ~n aa prior to the atmospheric entry and the be- 
ginning of the computation of the entry targeting data. This computa- 
tion defines a total range to target (RT) and the current heading to 
target The equations used for targeting are discussed in 

section 5- 5-2. Until 0.05g> the spacecraft will be in a three axis 
attitude hold i mx .. At 0.05g, rate damping will be initiated and the 
guidance will transfer to the constant heat rate phase. 


5.2 Constant Heat Rate Phase 

During the constant heat rate phase a stable trajectory is estab- 
lished at an acceptable temperature prior to the initiation of ranging. 

A roll command is generated which will control the spacecraft along a 
desired constant temperature profile through pullout (y = 0). This 
phase is required to stabilize the trajectory prior to entering the 
equilibrium glide phase. The equilibrium glide ranging phase is entered 
after pullout at an inertial velocity of 25 000 fps. Appendix A pre- 
sents the derivation of the guidance equations for the constant heat 
rate phase . 


5.3 Equilibrium Glide Ranging Phase 

At an inertial velocity of 25 000 fps, the guidance enters the 
equilibrium glide ranging phase. During this phase entry range pre- 
dictions and reference trajectory parameters are computed which are re- 
quired by the trajectory controller to eliminate range errors. However, 
trajectory control is not transferred to the equilibrium glide mode 
until the drag command from the reference equilibrium glide profile is 
greater than the drag command from the constant heat rate phase. This 
control mode transfer prevents a discontinuity in the total drag refer- 
ence trajectory, thus eliminating an undesirable transient in the tra- 
jectory. 

Closed form solutions of the equations of motion are used to pre- 
dict the entry range and the reference trajectory parameters. These 
parameters are based upon an equilibrium glide flight at constant bank 
angle. If the equilibrium glide profile is flown at low speeds, higher 
than desired load factors may result; therefore, the trajectory profile 
is based upon a constant load factor starting when the load factor 
reaches 1.5g. The range prediction is accomplished by analytically 
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predicting the inertial velocity at which the equilibrium glide tra- 
jectory will reach a total load factor of 1.5g (V GG ), and then 

analytically predicting the range from the current orbiter velocity 
to V n _ based upon an equilibrium glide trajectory, and analytically 

predicting the range from V GG to transition by assuming a constant 

1.5g trajectory. The equilibrium glide roll angle is selected to make 
the resultant range prediction equal to the current range to the target. 

Once the desired equilibrium glide roll angle has been determined, 
a reference trajectory is analytically computed and a reference vertical 
L/D, a drag reference, and an altitude rate reference is computed and 
sent to the controller. The controller then computes a vertical L/D 
command based upon the difference between the reference drag and alti- 
tude rate commands and the actual trajectory drag and altitude rate. 

This vertical L/D command is converted into a roll command in the 
ROLL COMMAND service routine (section 

A new equilibrium glide roll angle is computed each pass through 
the guidance logic until the constant g ranging phase is entered. In 
addition to the equilibrium glide and constant g reference trajectory, 
a constant g reference profile is analytically computed based on the 
constant g level required to reach the target from the current spacecraft 
velocity. This constant g level is compared to the g reference level 
from the equilibrium glide trajectory. When the equilibrium glide drag 
reference is greater than the constant g reference profile required to 
reach the target, the equilibrium glide phase is terminated and control 
is transferred to the constant g ranging phase. Appendixes B and C 
present the derivation of the equations used in the equilibrium glide 
ranging phase. 


5.U Constant g Ranging Phase 

The constant g phase predicts the constant g level required to 
reach the target and then computes the reference parameters required by 
the controller to fly the desired constant g profile. The range pre- 
diction' is based on an analytic solution of the equations of motion which 
predicts the range flown from the current velocity to transition (assumed 
to start at Mach 6). The L/D reference, the desired drag reference, and 
the altitude rate reference is computed and sent to the controller. The 
constant g phase is terminated at the velocity for beginning transition. 
Appendix D presents the derivation of the equations used in this phase. 
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5.5 Service Routines 

Four service routines are used by the guidance system: controller, 

targeting, lateral logic, and roll command. 

5-5-1 Con troller .- The controller generates an L/D command in the 
vertical plane based upon the reference L/D, the reference drag level, 
and the reference altitude rate computed in the guidance phases pre- 
viously described. The basic controller equation is defined as follows. 


~ command = £ ref + Cl(D - D J + C2(R - R ,) (l) 

Dy Dy ref ref 

The constants Cl and C2 vary depending on the particular guidance 
phase . 

5.5.2 Targeting .- The targeting program computes the total range 
to target, the spacecraft heading to target, and the initial roll direc- 
tion. These computations are made in the earth relative coordinate 
system. The total range is computed as the great circle range between 
the present vehicle position and the target position. As shown in 
appendix E, the current heading to target ^ is computed based upon 

the current position and the target position. Knowing the heading to 
target, the initial roll direction is chosen to reduce the angle be- 
tween the present heading and the heading to the target. 



5.5-3 Later al logic .- The lateral logic consists of a lateral 
deadband about the spacecraft heading. When the magnitude of the 
difference between the spacecraft heading and the heading to the target 
exceeds the lateral deadband and the roll direction is such that this 
difference will increase, the guidance commands a roll reversal. The 
azimuth deadband method of cross-range control was chosen because a 
cross-range deadband technique will cause a high L/D vehicle to spiral 
above Mach 1. Direct control of azimuth eliminates the spiral. For 
vehicles with a low roll response, it may be necessary to prevent a roll 
through negative lift at high g levels. This capability has been in- 
cluded in the guidance logic as presented in appendixes E and F. 

5.5.^ Roll and alpha command . - The roll and alpha command subroutine 
generates angle of attack and roll commands for the autopilot. This 
subroutine also converts the vertical L/D command from the controller into 
a. roll command. The direction of the roll command is determined by the 
lateral logic. 
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6.0 CONCLUSIONS 


An entry guidance logic for preselected angle of attack trajectories 
has been developed and initial studies using this guidance demonstrate 
excellent performance. This guidance logic combines control of load 
factor and temperature with ranging by means of an analytically computed 
reference trajectory. Analysis of this guidance concept has indicated 
the following: 

a. Closed loop ranging can be provided by an analytical guidance 
logic while implicitly controlling temperatures and load factor. 

b. The guidance system affords at appropriate times close control 
of all critical constraints (i.e., temperatures, load factor, and heat 
load) . 

c. The closed form range predictions afford fast computational 
capability which is desirable for an onboard guidance system. 

d. Preliminary navigation error analysis indicates that this system 
is sensitive to navigation system errors. 
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APPENDIX A - CONSTANT HEAT HATE PHASE 

The constant heat rate phase computes a reference trajectory which 
is used until the ranging solution from the equilibrium glide and constant 
k phases is valid. The purpose of the constant heat rate phase is o 
stabilize the trajectory at a constant temperature during the 
entry into the atmosphere prior to the initiation of ranging which begins 
at an inertial velocity of 25 000 fps. This reference trajectory con- 
sists of a vertical L/D reference, a drag level reference, and an i- 
tude rate reference. These reference trajectory parameters are used by 
the controller during the constant heat rate phase. 

Stagnation point heat rate for a 1-foot radius sphere is defined as 


Q = 17 60C W(F 


Specific aerodynamic drag is along the negative velocity vector with 
the magnitude computed as follows : 

pV E C D S l a: 


Solving (Al) for p 


(17 600) 2 (V e /26 OOO) 3 * 15 


Substituting (A3) into (A2) gives 


D = D 


26 0 00 6 ‘ 3 p o C D SQ 2 
ref " (2 x 17 600 2 )mV E U ‘ 3 


Equation ( AU) provides an expression for constant heat rate in te 
a reference drag force. The reference drag is used in the controller. 
The altitude rate reference term used by the controller can be e 
as follows: 


S59-U 




2k 


Assume p = 


-H/HS 


_ 3^ 3H 

aH at 




- HsP " “ HS P 


»VS 8 


V¥ 

Then taking the derivative assuming that 1 b a constant gives 


1+mDV 2mD 


Since V ■ -D 


Since R ■ - %S 
P 


2fil + £ D 

V E D 


£.. . 2V + D 


p v E D 


"■- HS (^ + 5 


* ref ^ref 

R . - -Hs[ 

ref V V„ D . 

\ E ref 
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Equation (A4) gave 


26 000 4,J p C SQ 2 3 

~ . Kiov **-3 

(2 x 17 600 2 )mV E 4 ' 3 


where 


26 000 6 ’ 3 p o C D SQ 2 
jqO = ~~Jj2 x 17 gOO 2 !^” 


(A10) 


\ _ ref 3V 

'ref " " 3V 3t 


D , » -U.3K10V " 5 * 3 v 
ref h 


' >1.3K10V E - 5 ' 3 D ref 


' , - 3D ref 2 


(All) 


Substituting (Al+) and (All) into (A8) gives 


(‘ 

= -HSl- 
■ -HSl- 


2D ref D ref 
, V E D ref 


,2D ref , U - 3D ref 
v V E V E > 


R - -6.3HS-— 1 
ref V E 


The nominal L/D required to fly the desired profile, ^ 

*1 


(A12) 


reference 


is derived in the following manner: 


V 2 cos y 

Vy “ + L - g cos y 

K 
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or 


V Y = 


cos y 
_ 


therefore , 




w- 


cos y 


g cos y 


g cos y " 




+ VyJ/d 


Assume cos y = 1» R S = v sa ^ 2 


(A13) 


Since 


i. . . AA * xx 

D v ®V D 

h = V sin Y ~ V Y 
h = Vy + Vy 
Vy ■ h - Vy 


for constant heat rate 


therefore 


h = -6.3HS- 

5 = -6.3Hs(%-^) 

v-f * -6.3HS^- p)- vy 


(AlU) 


(A15) 


(Al6) 
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.(-6.3h4)(- £)-6.3h4 


Combining (Al6) and (A17) gives 


Vy = -6 


•38s(| 


+ A 


Therefore 




However, since D - U. 37 - for constant heat rate 




3^z* + 


i- 


- 39.69HS^5- 


Therefore 


V ef D ref 


1 - r~~J - 39.69HS-£f£ 
k V sat 2 / V 


(A17) 


(Al8) 


(A19) 


(A20) 


(A21) 


(A22) 


Evaluating the second term in the m- equation for the constant 

ref L 

angle of attack case produces a maximum change in - of 0.0097 units 


and for the variable angle of attack case 0.0216 units. Since this term 
is negligible 
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(A23) 


Equations (a4j, viu2), and (A23) provide the D reference, R ref- 
erence, end L/D reference that are required "by the controller to maintain 
a constant heat rate trajectory. Figure A-l shows a time history of the 
commanded and actual heat rate during the constant heat rate phase. 


S59-19 







29 



Figure A1 .-Comparison of commanded and actual heat rate. 



APPENDIX B - EQUILIBRIUM GLIDE PHASE FOR 
A CONSTANT ANGLE OF ATTACK 


The purpose of the equilibrium glide phase is to predict the range 
capability of the spacecraft and to compute a reference trajectory which 
will terminate at the target point. This Is accomplished by predicting 
analytically the range flown from the current orbiter velocity to the 
velocity at which a load factor of 1.5g is reached. Then the resultant 
range for a constant 1.5g trajectory is predicted in the remainder of 
the entry. The initial range prediction assumes an equilibrium glide 
trajectory with a roll angle selected to correct for range errors. Once 
the equilibrium roll angle has been predicted that will satisfy the 
range requirements, a reference drag trajectory is commanded that will 
correspond to the desired equilibrium glide trajectory. 

The basic equilibrium glide equation is 

V* cos y 

Vy ■ p + L - g cos Y (Bl) 


For equilibrium glide, y B 0; therefore 

cos y 

0 ■ r + L - g cos Y 


Assuming coa y a 1, equation (Bl) reduces to 



. Since Rg = V^^ 2 , ‘ thie equation reduces further to 



(B2) 


(B3) 
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Solving for D gives 


D 




(B4) 


Since L/D in the vertical plane a L/D * cos (f>, equation (Bl) 13600063 


L 

- cos 4 



(B5) 


Using equation (B5), it is possible to predict the range that will 
be flown during the equilibrium glide phase by means of the following 
equations . 


Assume that the equilibrium glide trajectory will be based on a 
constant roll angle, 4 , and will be flown to the inertial velocity at 
which the predicted trajectory reaches 1.5g (V CG ) . V CQ can be predicted 

by solving for Vj in equation (B5). 




D V 


2k 


cos 4 


(B 6) 


Where D is the drag along the velocity vector equivalent to l,5g 

'■'o 


D a — 

Cg ^/l + (L/D) 2 


(B7) 


Equation (B 6) is valid for all equilibrium glide roll angles that result 
in trajectories that reach 1.5g« However, trajectories based on small 
equilibrium glide roll angles do not obtain 1.5g. For this class of 
trajectories, the guidance can determine this by checking for a negative 
square root in equation (B 6). When this occurs, the guidance must assume 
that the constant g phase is eliminated and the equilibrium glide trajec- 
tory required to reach the target is flown all the way to transition at 
Mach 6. 
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The range fir v . ;ig the equilibrium glide phase can be predicted 
by the following c u v :ious: 


SR _ 9R 9T 
SV “ ST SV 


V 

D 


ref 


(B8) 


Using equation (B r .) 


r- cos 4> V . 2 
D sat 


frt 2 ) 


(B9) 


(L/D) cob <(> V ■ 2 f CG V_ 

R = _ _ . sat . / v2~T~v T *V 

5 J I sat 


(BIO) 


Integrating equal 1 tf-10) 


W»i CO. ♦ /V * - V » 

R “ bpz/ ■ 


(Bll) 


The range from V „ to transition can be analytically predicted by the 
equations 


_SR _ SR ST 
SV = ST SV 


V V I 
D " ** D 

eg 


(B12) 


V, 


R = 


eg J 


TRAN 


V dV = 


V 2 - V 2 
CG TRAN 


2D 


(B13) 


eg 


CG 
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% 

4 


3 


Therefore equations (Bll) and (B13) represent the total predicted range 
for the entry from the current orbit er velocity to transition. 

R = R + R (BlU) 

M EQ CG • 

A comparison between Rp and the actual range to the target 

(assumed to be the transition point) will produce a range error which 
can be nulled by changing 4> > the equilibrium glide roll angle. 

Figure B-l presents the range correction capability as a function of . e 
equilibrium glide roll angle. This figure shows that for an equilibrium 
glide roll angle below 1+3.5°, the equilibrium glide trajectory will not 
intersect 1.5g. Thus for targets that require these roll angles, an. 
equilibrium glide trajectory will be flown throughout entry. This figure 
also shows that for large equilibrium glide roll angles (to the right of 
the line marked V CQ greater than in fig. B-l) , the desired equilib- 

rium glide roll angle will intersect 1.5g prior to the current velocity. 

So for these cases, the guidance will immediately transfer into the 
constant g ranging phase whenever \ r CQ is computed to be greater than 

v r 

Once the equilibrium glide roll angle has been determined, the con- 
troller reference parameters must be computed in order to fly the desired 
equilibrium glide trajectory. The controller requires a L/D reference, 
a drag reference, and an altitude rate reference. The drag reference 
term is simply equation (B5). 


A. 


ref | cos * 


V I 


The L/D reference term is simply 


L ^ L 
- ref = - cos 

where ~ cos 4> is the inplane ^ required to reach the target. 

The altitude rate reference can be derived as follows: 

From equation (A8) . 

ref \ V E D ref. 


(B15) 
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V * * 


p B — JL— 
ref L 

- cos $ 


v ! 

Vt 7 


Taking the derivative of D 


2 * D ref T I 


ref L , „ o 
D C08 * V BOt 


(Bl6) 


Combining equation (A8) with equations (B5) and (Bl6) gives 


L 

D C0S 



(B17) 


Equations (B5), (B15) and (B17) are sufficient to establish an equilib- 
rium glide trajectory. 
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APPENDIX C - CONSTANT g PHASE 


The purpose of the constant g phase is to predict the constant g 
level required to reach the target and to generate a ^ re f» ^ re f* an( * 

a ^ for the controller, 
ref 


Equation (B 13) presents the equation that analytically predicts 

the range that will be flown if a constant g profile (D ) is flown 

eg 

between and transition. This equation is as follows: 

V 2 - V 2 
CG TRAN 

R CG * 2D 

eg 


The range to the transition point is obtained from the targeting logic 
and Is equal to the total range to target minus the desired range to 
the target at transition 


^TG ” R : 


PT 


The constant g level to reach the target becomes 


D ■ 

o 


V 2 - V 2 
V CG V TRAN 


2R, 


(Cl) 


TG 


The constant g trajectory is controlled by means of the drag controller 
where 

D = D 
ref o 


R . * -HS 
ref 


2 D ref + D ref 
V E D ref> 


For const ant g D 


ref = 0, therefore 


D 


R « -2HS-— 
ref V E 


(C2) 
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As was the case for constant heat rate, a L/D reference term can he 
derived from the equation of motion 


VY 


Y = ^SSSJL + L - g cos y 

Jrv 


or 




D„ D 


V 2 
v sat 


and from equation (A15) 


For constant g 


Vy ■ H - Vy 


H ■ -2HS y ■ 


H “ -2HS 


D DV\ 
V 2/ 


(£- 

\V 


-2HS 


V2 


vy * -2HS - Vy 


V - ^ - (- £HE v)(- v) 


2HSD 2 
V 2 


Vy “ -toiS “J 


Therefore , 


i m J l . J !L,\ . 

e v »\ wy 


HSD 
+ V2 


so 


v ef 


D ref 


(>-& 


ref 


27 - UHSr vT 


(All*) 


(C3 ) 


(C4) 


(C5) 


:(C6) 


T7) 


( C8) 
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1 

3 

\ 

| Relative velocity was assumed for the constant g phase because of the 

» requirement to switch from inertial velocity to relative velocity late 

f in the entry when velocity is approximately equal to V /2. 
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A’vFJflVlX. D - MATH FLOW LOGIC FOR THE CONSTANT 
. ANGLE OF ATTACK GUIDANCE 


All necessar v .- .auctions have been developed In the main text and 
in appendixes A tiao..-i. C. The final step Is to connect these equations 
with decision logic to convert trajectory data from the navigation into 
a commanded roll end a commanded angle of attack for the autopilot. 

This appendix pr-vinu*? the guidance flow logic and all necessary equa- 
tions and constai ' . . ’’ v _ the constant angle of attack case. Table D-I 
presents the constants and initial variable values for the guidance. 

Flow charts 1 through 13 presents the math flow logic for the guidance. 



.3 

9 

1 


i 

4 






TABLE E-I 
i 

ALMN 
C„ 


CONSTANTS AND INITIAL VARIABLE VALUES 


D 

: fnm 

: rnm 


1.5 

0.336 

0.000164578836 

3437-7468 

32.2 
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g 

n.d. 

n. mi. /ft 
n. rax. /rad 
ft/ sec 2 


G2 

2.5 

g 

HS 

28 500. 

fps 

H 

1 

n.d. 

up 



IFT 

0 

n.d. 

IG 

0 

n.d. 

ISTP 

1 

n.d. 

LATSW 

0 

n.d. 

LOD 

1.497 

n.d. 

*c 

80 . 

Btu/ft 2 / 

R 

21 04 1 776 . 

ft 

RPT 

182.4 

n. mi. 

RTD 

57.29577951 

rad/ deg 

SELECT 

1 

n.d. 

VQ 

6443. 

fps 

VS 

• 25 766.1976 

fps 

VSW 

25 000. 

fps 

Vxx 

6443. 

fps 

WE 

.72921149 E-4 

rad/ sec 

WT 

i 4 o 000. 

lb 

y 

20. /RTD 

rad 

IT 

P 

O 

3.14199265 

0.076474 

n.d. 

lb/ft 2 

05GSW 

0 

n* d. 








Flow chartDl Overall flow logic. 
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Flow chart (2,- T agisting - Corwluckd. 
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F low chart d 3 . - Se lector. 


S59-37 







'*?? 'f . 


S59-3S 





















































o|r~ 


72 



Flow chart d 7,- Controller. 
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Submit tal 53 -Approach Guida nce 
1. INTRODUCTION 

The Approach Guidance Routine presented here is designed to take the orbiter 
vehicle from the end of the entry phase (altitude * 100, 000 ft) down to the start of the 
terminal guidance phase 5. 9 n. mi. from the runway at an altitude of 6900 ft. and a 
velocity of 490 ft/sec. It is based on the ideas of Refs. (1) and (2). 

The guidance routine consists of six modes: Acquisition, Energy Dissipa- 
tion, Turn-in, Initial Approach, Heading Alignment, and Final Approach. The 
horizontal geometry is illustrated in Figure 1 in which the circled numbers refer 
to the various modes. The Acquisition Mode begins at 100.000 ft altitude, contains 
an angle- of- attack transition maneuver, and ends when the vehicle is within about 
15 n. mi. of the runway. Energy dissipation involves flight in the vicinity of the runway 
around a cylinder of radius 13. 5 n. mi. During this mode the vehicle descends from about 
50,000 ft altitude to 26,000 ft. This helical flight usually comprises less than one 
half of a revolution around the cylinder. The next three modes, i.e. Turn-in, Initial 
Approach, and Heading Alignment, constitute a two-turn maneuver to place the ve- 
hicle on the appropriate final approach path. The Final Approach Mode establishes 
the- proper interfaces with the Terminal Guidance Routine for the final maneuvers 
required to land on the runway* 



i 





Figure 1* Horizontal Geometry, Approach Guidance 
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2. functional flow diagram • 

The basic flow of the Approach Guidance Routine is shown in Figure 2* 

After the routine is entered and initialized, targeting computations are made 
to obtain the current values of position and velocity, and direction parameters of the 
vehicle relative to the desired touchdown point. Next, the mode is selected based 
on the current trajectory conditions, and quantities unique to the specific mode 
are computed. 

The angle-of- attack command is used for vertical control and is computed 

during the first part of Mode 1 so as to accomplish a constant (-0. 3 deg/sec) angle- 

of- attack transition maneuver. During the remainder of Mode 1 and for Modes 2 

lb . . . 

and 3, an angle-of-attack which will yield a constant (210 g ^ dynamic pressur 
is commanded. Finally, during the last three modes, the angle-of-attack which 
will cause the vehicle to fly at a constant flight path angle (-11 deg) is commanded. 

The roll-angle command is used for horizontal control and is computed for 
the various modes as shown in the following table* 

Table 1 

Geometric Criteria for Roll- Angle Command 

Mode 

1 
2 

3 

4 

5 

6 

Finally, the rudder flare or speed brake is deployed during the last three 
modes in order to achieve a speed of 480 f/ s at the end of approach guidance. 


Geometric Criterion 

Tangent to Energy Dissipation Cylinder (EDC) 
Fly on EDC 

Turn toward center of EDC 

Tangent to Heading Alignment Cylinder (HAC) 

Fly on HAC 

Align into vertical runway plane 






Initialization 



Figure 2. Functional Flow Diagram, 
Approach Guidance 
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NOMENCLATURE 
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Notational Conventions 


Upper-case letters represent matrices 

Lower-case and Greek letters reserved for scalars and vectors 
Vector quantities are underlined, e. g. x 

Vectors are assumed to be column vectors unless explicitly noted 


Symbols 

a Effective aerodynamic area of vehicle 


c l 


cl 


RT 


turn • 
g 
h 
h 
h 

hi 




Acceleration (RW coordinate) 


Coefficient of lift 


d(c£ )/d(a) 


Rudder-flare command 


Distance between touchdown point and ED center 


Turning factor 


Gravity 

Altitude 


dh/dt 

dh/dt 


Altitude at beginning of Initial Approach 


Reference altitude 


d(h R )/ 


dt 
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CoT-'rc-i gain 


’ r; -: ol gain 


j_,elt or right HA selector 


Lr A lit specific force 


\ ehicie mass 


Mach no* 


M Run-way to NED Coordinate Transformation Matrix 

R-NED 


pressure 


Dynamic pressure at last guidance call 


Desired dynamic pressure 


dq/ dt 


dc / ^ 


Radians into degree 


Distance from vehicle to Tacan 


Radius of ED cylinder 


— NED 


Horizontal component of vehicle position vector (NED-Tacan Coord. ) 


Horizontal component of position vector (RN-Tacan coord. ) 


Integration of steady state heading error in Mode 2 


Mode selector 


Clock time 
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At 


Hi 
-2 
H 3 


Hned 


Hr 

Het 


y 

y 


a 


O' 


0 
0 , 


Ao 


D 


A l/> . 


* 


c 


"RW 


HA 


Time interval between guidance updates 
Unit vector (1,0,0) 

Unit vector (0, 1, 0) 

Unit vector (0, 0, 1) 

Vehicle velocity w. r. t. air mass 

Horizontal component of vehicle velocity (NED Coord, ) 

Vehicle velocity vector (R.W. Coord. ) 

Horizontal component of vehicle velocity (R, W, - Tacan Coord. ) 
Cross range to runway 
dy / d t 


Current angle-of-attack 


Angle-of-attack command 


Roll angle command magnitude 


Roll angle command 


Desired angle-of-attack change 


Desired azimuth change 


Azimuth ( 0 ~ 360 


Azimuth command ( 0~360 c ) 


Runway azimuth 


Azimuth of (’Pp^R,) vec * or 
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a 

■a 



v*v 

j 

J 


i 





p 

Pp 

— HAN 

— HAR 
i*HAR 
°R 

— HAT 


Angle to Tacan w. r.t. north (from vehicle) ( 0- 360°) 

Angle to Tacan w. r. t. velocity (from vehicle) ( 0~ 180° ) 
Air density 

Predicted range to go to key point 

Vector from proper HA to vehicle (NED Coord, ) 

Vector from proper HA to vehicle (RW-Tacan Coord, ) 
Vector from other HA to vehicle (RW-Tacan Coord, ) 
Vehicle position vector in (RW Coord. ) 

Distance to touchdown point 

Vector from Tacan to proper HA center (RW-Tacan Coord. ) 


Coordinates: 

(RW Coord, ) 

Runway coordinates, centered at touchdown point 

z in runway landing direction, i. e. directed down-range and forward 
x up 

y x, y, z from right hand orthogonal 

coordinates 

(RW - Tacan Coord. ) 

Runway coordinates, centered at Tacan or center of EDC 
(NED Coord. ) 

Eocal North, East, down coordinates at point of Tacan 


i 


4 



t 


Angle Measurements 

(0,360°) or(0< 0<36O°) 

0 between 0 and 360°, measured clockwise 

(0,180°) or (-180° <0 <180°) 

0 between ( 0 180 ) if measured (clockwise) 

( 0 -180 ) if measured (counterclockwise) 
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Vehicle >n vector (RW Coord.) 

Runway azimuth ( 0 - 360° ) 

Time into v f i-rSveen guidance updates 
- Variable s 

Command i oil . ngle 
Command angle-of-attack 
Command rudder flare 






pet. mi j:d p low diagrams 

This section contains detailed flow diagrams of the Approach Guidance Routine 


UNIVERSAL CONSTANTS INPUT VARIABLES 



X 


0 

>— < 

H 

< 

to 

j 

SJ 

t-> 

S3 

f— i 

1 



i 

I 

I 


I 

I 


§ 

»— i 

H 

W 

O 


Figure 4a, Detailed Flow Diagram, Approach Guidance Routine 
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(Fig. 4d) (Fig. 4e) 


(Fig. 4f) (Fig. 4g) (Fig.4h)* 



Figure 4b. Detailed Flow Diagram, Approach Guidance Routine 
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MODE SELECTION 
(ENERGY DISSIPATION MODE) 






ES 








i Flow Diagram, Approach Guidance Routine 
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Figure 4h. Detailed Flow Diagram, Approach Guidance Routine 
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Figure 4k. Detailed Flow Diagram, Approach Guidance Routine 
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Subndttal 60: Energy Dis sipati on Ra te Guid ance 



1. INTRODUCTION 

The Approach -Guidance Routine presented here is designed to take the orbiter 
vehicle from the end of the Entry Phase (altitude «100,000 ft) to the start of the 
Final handing maneuver (altitude ^7000 ft). A detailed description of the guidance 
concept along with simulation results demonstrating its feasibility is given in Ref. 
(U. 

The Approach-Guidance system is a closed-feedback-loop scheme. The 
vehicle energy is managed by controlling the rate at which energy is dissipated 
during a straight-in approach flight. Energy dissipation rate is controlled by 
flying at a constant value of dynamic pressure and varying the vehicle’s lift to drag 
ratio with the Rudder Flare and or other available drag-increasing devices (e. g. 
body flap). 

The complete approach flight consists of straight , fixed-length segments 
from the vehicle’s initial position to the airport’s main navigational facility (VOR, 
TAG AN) or a suitable artificial checkpoint, then to a point in the final approach 
plane (intersection point) where the final flight path is intercepted, then straight 
towards the runway until the Final Landing Guidance System takes over (Outer 
Marker). Constant-bank turns link the straight flight segments. 

The closed-loop energy management policy automatically compensates for 
any wind component that may affect the energy dissipation rate of the vehicle. 



INTERSECTION 

POINT 
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NOMENCLATURE 


Notational Conventions 

Lower-case and Greek letters reserved for scalars and vectors 
Vector quantities are underlined, e. g, x 
Components of a vector x are denoted 


Symbols 

a TAC 

d INT 


d OM 


d TAC 


e 


e OLD 

e 

e D 
e IN 
e OUT 


Azimuth of VORTAC 

Distance from touchdown point to final approach plane 
intersection point 

Distance from touchdown point to initiation of landing 
guidance system 

Distance from touchdown point to VORTAC 
Current vehicle energy 
Desired vehicle energy at target points 
Value of vehicle energy at previous guidance pass 
Current vehicle energy dissipation rate 
Desired value of vehicle energy dissipation rate 
Value of e D at which flyover mode is entered 
Value of at which flyover mode is left 
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I 


$ 




MIN 


MAX 


-DH 


-ivH 


K P 

V V t 3 




-RT 


— 1 ’— 2 ’ — 3 


— DH 


-OLD 


'FO 


Upper limit of preferred e ^ region 
Lower limit of preferred e ^ region 
Gravitational constant at sea level 
Current target index variable 

Unit vector in direction of difference between horizontal 
components of vehicle and target position vectors 

Unit vector in direction of horizontal component of velocity 
w. r. t. touchdown point 

Longitudinal channel coefficient 

Roll channel coefficient 

Vertical (pitch) channel coefficient 

Horizontal distances from target points at which target 
switching is to occur 

Current value of dynamic pressure 

Desired value of dynamic pressure 

Vehicle position vector w. r.t. the touchdown point 

Target points position vectors w. r.t. the touchdown point 

Vector difference between the horizontal components 
of vehicle and target position vectors 

Vehicle position vector at previous guidance pass 

Switch indicating flyover mode is in effect 


S60-3 














Switch indicating this is the first pass 


Vehicle velocity w, r. t. the touchdown point 


Horizontal projection of velocity vector 


°CMAX 


Maximum permissible angle of attack 


XMIN 


Minimum permissible angle of attack 


Flight path angle 


Command Rudder Flare deployment angle rate 


Maximum Rudder Flare deployment angle rate 


Heading error 


Command vehicle pitch angle 


CMINT 


Minimum pitch angle during transition 


CMC 


Current lower limit for command pitch angle 


0 _ _ Corrected minimum limit for command pitch angle 

■ CiVI 

Maximum command pitch angle 
Command roll angle 

^ Maximum command roll angle magnitude 

Special Notation 

Sign ( ) Algebraic sign associated with ( )* Value is +1 or -1, 
with sign (0) - +1 . 

max ( ) The maximum of all element enclosed in ( ) 
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£ 2 func tional flow diagram 

; | The basic information flow for the Approach Guidance Routine is shown in 

Figure 1. This is based on the guidance concept of Ref. (1). 

The guidance task is made up of three independent channels: 

1. The pitch angle command is proportional to the 

difference between the desired and measured dynamic 



pressure. This command is limited so^as to limit 
the pitch- down during the transition maneuver and 
the maximum and minimum angles of attack. 

2. The desired energy dissipation rate is computed as 
the ratio of energy-to-be-dissipated to the distance- 
to-go. The Rudder Flare deployment angle rate is 
proportional to the difference between the desired 
and the actual energy dissipation rates. A discrete 
rate controller is superposed in order to drive the 
desired dissipation rate to a "preferred" value 
range. 

3, The roll angle command is proportional to the dif- 
ference between the horizontal velocity direction 
and the horizontal direction from the vehicle posi- 
tion to the target point. 

The three channels guide the vehicle sequentially to three target points. 


Switchover from one target to the next is made at a predetermined horizontal distance 
from the current target. 

The longitudinal channel (No. 2) is inhibited when the vehicle is initially too 
close to the first target (desired energy dissipation rate exceeds the maximum avail 
able). The guidance is then in "flyover" mode. Target switching is also inhibited 
in this mode. 


1, Elias, A., "New Approach Guidance Concept for Shuttle", 
23A STS Memo No. 58-72, 4 December 1972, MIT/CSDL. 
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Figure 1. Functional Flow Diagram for Approach Guidance 
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3. INPUT AND O UTPUT VARIABLES 

Input Variables 

Vehicle position w. r*t # touchdown point (touchdown point coordinates: 
x-up, z-dowr * ’^. .ge, y-crosstrack) 

~RT Ve M c l e veloc < w. r. t. touchdown point (touchdown point coordinates) 

q Measured value of Dynamic Pressure 

Output Variabl es 

Pitch angle command 

6^ Rudder Flare F flection angle rate command 
Roll angle command 
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4. 


DESCRIPTION OF EQUATIONS 


4 . l Initial Target Calculations 

During the initial pass, the target point horizontal position vectors ire con- 
structed from their distances to the touchdown point and the VORTAC azimuth 
w. r. t. the localizer direction: 

VORTAC : £j = (0, -d TAC sin ® TAC > "^ T AC COS a TAC } 

Intersection : £2 = (0, ~^INT ^ 

Outer Marker : r 3 = (0, 0, -d ) 

Also, the following initialization tasks are performed: 

1. Flyover mode switch off — G j?o = ^ * 

2. Initial pass switch on — s pp = 1 , 

3* Pitch limit = transition pitch limit (it is assumed that 
the guidance system is initiated with the vehicle flying 
on the back side of the L/D curve), 

4. The target index i is set to 1 (VORTAC) 

5. The "old" values of e and r RT are set to zero. 

This makes the back-difference algorithm of sec- 
tion 4, 4 invalid during the first pass, so commands 
are not issued until the second pass. 

4. 2 Preliminary Calculations 

At the beginning of every guidance pass, the flight path angle and horizontal 
vector from vehicle to target are computed: 

y = tan' 1 ^ - v RTg 

-DH = (0 * r i 2 ' r RT 2 ' r i 3 ' r RT 3 } 

4, 3 Vertical Channel 

The command pitch angle is: 

0 C = k P (q D ’ q> 
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The lower limit for the command pitch angle is the largest of: 

1. The current absolute pitch minimum, 0 

2. The pitch angle corresponding to the minimum 
angle o* -mr-k, y+ ff CMIN • 

The upper limit for kind pitch angle is the pitch angle corresponding to 

the maximum angle of attack, y+ a q]vjaX * 

The current absolute pitch minimum is set to ** *<n.*ition 

maneuver, and to an arbitrary low value (e. g. -1 rad. ) after the dynamic pies- 
sure reaches the desired value for the first time (i. e. at the end of the transition 

pitch-down). 

4.4 Longitudinal .-i_ 

The current vehicle energy is computed from the position and velocity: 
e = " -RT ' -RT J / 2g 

The current energy dissipation rate is computed as the back-difference: 

■ - c " e OLD _ 

/ (r RT 2 ' r OLD 2 )Z ‘ (r RT 3 ' r OLD 3 ) 

then the values of an< ^ Z.OLD are u P^ a ^ e( ^* 

The desired value of ’the energy dissipation rate is then calculated. 




e D = (e i “ e) / I -DH | 

and a Rudder Flare rate is commanded proportional to the difference between 
desired and actual dissipation rates: 

^FtC = k E (e D " e) 

this rate is then limited to the -6 RCIj , 6 RCL range. This command is over- 
riden if the value of e^ falls outside of a preferred range: 

if e D >e MIN , 6 rc = -6 RCL 
if e D ^ e MAX’^RC = 6 RCL 






The "flyover" mode is enabled when e^ is lower than , and is disabled when 
reaches * When the "flyover" mode is enabled, the command rudder 

flare rate is zero. 


Lateral Channel 

Two unit vectors are computed: 


-VH I v T 


is in the current vehicle heading direction 


i nTT = is in the direction from the vehicle to the target 

~ DH IT r^rr\ 


The dot product of these vectors is the cosine of the angle difference be- 
tween the desired and the actual headings. In order to resolve the sign indeter- 
mination of the cos" 1 function, the single component of the cross product 


-DH X -VH = x DH 2 Vh 3 " X DH 3 X VH 2 


is computed. This is the sine of the heading difference, and its sign is used to re- 
solve the indetermination. 

This command roll angle is then calculated: 

0 C = k L 

This command is limited to the - ^CMAX ran S e * 

4. 6 Target Switchi ng 

The target index variable, is incremented when the horizontal distance to 
the current target reaches the target’s switching value. When the last target's 
switching distance is reached, (i = 4), the Approach Guidance is terminated. Tar- 
get switching is inhibited during the "flyover" mode (s^q s 1), 

Commands are not issued during the first pass. 
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Figure 2b. 


Approach Guidance Routine, 
Detailed Flow Diagram 
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Su bmittal 2P: Fi » al Vp proach Guidance 
Introductio - 

These eqe’Aoas are submitted as candidates to fulfill the unpowered Final 
Approach Guidance v eqi:;^ - ants for the space-shuttle Orbiter. They include 
Autoland lateral me lonp'T"Hinal guidance equations. The scheme is all inertial; 
navigations aids a: ^ v’ . Iy to update the navigated vehicle state. Pitch rate 

and speed-brake commands are computed and issued to control in-plane approach. 
Lateral position error and its integral plus heading-angle error are used to form 
the vehicle roll command. (There is no decrab or wings level maneuver; the 
assumption is mad; * that the gear is designed to accommodate the stress of ci ab 
bed landings in design winds). 

Functional 

Figure 1 is functional diagram. Figure 2 is a block diagram, (For gen- 
eral information, the autopilots being used in simulation runs are included in 
Figure 2. ) 

Inputs to the Guidance module are from the Final Approach and Guidance 
Navigation module; the inputs are the navigated state in the Earth-fixed landing 
coordinate system. From this are calculated the range to touchdown target, 
altitude, velocity magnitude, flight-path angle lateral position and heading angle. 
Outputs are pitch cate command, speed -brake position command and vehicle roll 
command to the ar pilot. The guidance roll command drives a roll-rate aileron- 
autopilot inner loop with- roll attitude outer loop. Roll rate command is intercon- 
nected to a rate command rudder autopilot with turn coordination and normal 
acceleration inputs. The acceleration and heading-angle signals are instrumental 
in holding the orbiter to the final approach plane in crosswinds. 

Coordinate System 

The autoland guidance uses vehicle position and velocity relative to a run- 
way coordinate system, as shown in Figure 3. Figure 3 also indicates longitud- 
inal sign convention for the equations. The "altitude of the IMU" at touchdown is 
represented in the .equations as touchdown altitude. 

Equations and Flow 

Figure 4 presents the detailed guidance equations. Autoland guidance is 
initiated with the vehicle established on the final approach path near the plane of 
the runway at 3000 to 10000 feet altitude. It is currently entered 8 times per 
second although little performance degradation is evident at half that frequency. 


2 
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On the first call, an initialization and targeting section is entered. Target- 
ing variables are used to define the flare, shallow glide, pull up and steep sections 
of the reference trajectory, A steep reference flight path angle is calculated such 
that the trajectory passes through the initial vehicle position. If, during the steep 
phase, navigation updates cause large vehicle altitude errors, the steep portion 
of the reference trajectory is retargeted to pass through the new vehicle position, 

A linear desired velocity profile is also computed from the vehicle's current 
velocity to a target value at the beginning of pull up. 

Reference and actual values of h and V are differenced and drive the guidance 
loops shown in Figure 2. Since altitude is approximately equal to the integral of 
Vi , the velocity term in the denominator of the altitude error gain makes that 
loop insensitive to velocity variations. The inner loop controls V which is pro- 
portional to h and provides damping for the outer loop. It is compensated with a 
second order digital filter which effectively cancels two undesirable pole-zero 
pairs arising from the autopilot -vehicle dynamics. This allows stable operation 
of the inner loop at a higher gain level and tighter closed loop control. The 
accuracy of the autoland maneuver is improved by injecting the open loop pitch 
rate commands i and q . The q signal is composed of three parts (q cy , 

1 b u u u u 

Q , q ) all of which are tuned to the specific vehicle being flown. The q 
M cm r*cge ... 

term is a sinusodial pitch rate term added during pull up and again, with a 

different amplitude and period, during flare to lead the vehicle through these 

maneuvers. The q cv term ramps up to a constant value after the pull up maneuver 

and provides the increasing angle -of- attack necessary to maintain lift as the 

vehicle decelerates along the shallow slope. The c l C g e term ramps down to a 

constant value during the flare maneuver which helps the vehicle drive through 

the ground effects and minimizes runway float. Typical plots of these terms are 

sketched in Figure 2. 

Lateral Guidance 

The lateral guidance is all-inertial. A decrab maneuver was not studied; 
the assumption being that the gear is stressed for crabbed landings in design 
winds. The roll gain is halved during flare which levels the wings somewhat 
in steady-state crosswinds. * 

The lateral guidance equations are presented in Figure 4. On the initial 
pass, the roll gain, crossrange integral gain, and the heading gain are stored. 

On a normal pass the crossrange gain, Ky, is calculated as a function of velocity. 
When altitude becomes less than 50 ft, the roll command gain is decreased from 
6 to 3 over a 2-second period. The roll command is the sum of a crossrange. 
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integral of crossrange, and velocity heading angle term. It is limited and issued 
to the autopilot. 

Veloci ty Control 

The speed brake is commanded to a position proportional to the sine of the 
velocity error. Zero error is at 30 degrees brake for bi-directional control. At 
the beginning of pull rp the brake is completely retracted to eliminate pitch rates 
from transients near touchdown. 

Constants /Variables Summary 

Figure 6 summarizes variables and constants. 
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Figure 3. Runway Coordinate System 







































































FIGURE 4-7 
AUTOLAND EQUATIONS 
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Submittal 6: CONIC STATE EXTRAPOLATION 


i. INTROD UCTI ON 

Tin* Conic State Extrapolation Routine provides the capabil- 
ity to conically extrapolate any spacecraft inertial state vector either 
backwards or forwards as a function of time or as a function of 
transfer angle. It is merely the coded form of two versions of the 
analytic solution of the two -body differential equations of motion of 
the spacecraft center of mass. Because of its relatively fast compu- 
tation speed and moderate accuracy, it serves as a preliminary 
navigation tool and as a method of obtaining quick solutions for tar- 
geting and guidance functions. More accurate (but slower) results 
are provided by the Precision State Extrapolation Routine. 
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NOMENCLATURE 


Semi / ir axis of conic 


First ron-c param eter(( • v ^) /sfiP ) 
Secord . l \);y parameter ( / n - 1 ) 

Third conic parameter ( r^ ^ ) 
Power series in £ defined in text 


Eccentrll ' omaly 


True anomaly 

Hyperbo^io analog of eccentric anomaly 


Counter 


Maximum permissible number of iterations 


Unit vector in direction of r^ 


Unit vector in direction of v, 


Semilatu - % ectum of conic 


Normalized semilatus rectum (p / r ) 


Period of conic orbit 


Magnitude of 

Inertial position vector corresponding to initial time 


Magnitude of r (t) 

Inertial position vector corresponding to time t 

Switch used in Secant Iterator to determine vyhether 
secant method or offsetting will be performed 








S{| ) Power series in £ defined in text 

t Final time ( end of time interval through which an 



extrapolation is made) 

Initial time (beginning of time interval through which 
an extrapolation is to be made) 

Difference between specified time interval and that 
calculated by Universal Kepler Equation 

Magnitude of v^ 

Inertial velocity vector corresponding to initial time 
*0 

Inertial velocity vector corresponding to time t 

Universal eccentric anomaly difference (independent 
variable in Kepler iteration scheme) 

Previous value of x 

Value of x to which the Kepler iteration scheme con- 
verged 

Previous value of x 

c 

The n i-th n value of x 
Lower bound on x 
Upper bound on x 

Reciprocal of semi-major axis at initial point 
Normalized semi-major axis reciprocal (ar^) 

Angle from r^ to 




rt SiAi i 


Specified transfer time interval ( t -* t Q ) 


Value of the transfer time interval calculated in the 
Universal Kepler Equation as a function of x and the 
• lie parameters 


Previous value of At 


•X h 7 ‘ i-th" value of the transfer time interval calcula- 
ted in the Universal Kepler Equation as a function of 


the "i-th" value x i of x and the conic parameters 


Maximum time interval which can be used in computer 
due to scaling limitations . 


Increment in x 


A Primary convergence criterion: relative error in 


transfer time interval 


Secondary convergence criterion: minimum permis- 
sible difference of two successive calculated transfer 


time intervals 


Tertiary convergence criterion: minimum permis- 
sible size of increment Ax of the independent variable 


Transfer angle (true anomaly increment) 


Gravitational parameter of the earth 


Product of and square of x 


OGq, Coefficients of power series inversion of Universal 

Kepler Equation 











2 . 


FUNCTIONAL FLOW DIAGRAM 


I 






The Conic State Extrapolation Routine basically consists of 
two parts - one for extrapolating in time and one for extrapolating in 
transfer angle. Several portions of the formulation are, however, 
common to the two parts, and may be arranged as subroutines on a 
com pater. 

2. 1 Conic State Extrapolation as a Function of Time (Kepler 

Routine ) 


This routine involves a single loop iterative procedure, and 
hence is organized in three sections: initialization, iteration, and 
final computations, as shown in Fig, 1. The variable "x" is the in- 
dependent variable in the iteration procedure. For a given initial 
state, the variable ,f x M measures the amount of transfer along the ex- 
trapolated trajectory. The transfer time interval and the extrapolated 
state vector are very conveniently expressed in terms of "x". In the 
iteration procedure, "x" is adjusted until the transfer time interval 
calculated from it agrees with the specified transfer time interval 
(to within a certain tolerance). Then the extrapolated state vector 
is calculated from this particular value of "x", 

2 . 2 Con ic State Extrapolation as a Function of Transfer Angle 

( Theta Routine ) 

This routine makes a direct calculation (i. e. does not have 
an iteration scheme), as shown in Fig. 2. Again, the extrapolated 
state vector is calculated from the parameter ”x". The value of "x" 
however, is obtained from a direct computation in terms of the conic 
parameters and the transfer angle 6 . It is not necessary to iterate 
to determine M x fI , as was the case in the Kepler Routine. 
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ENTER 




Initialization 

(Compute Various Conic Parameters ) 

3 




Compute "x" Corresponding To The Specified 
Transfer Angle 6 


X 

Compute Transfer Time Interval Corresponding 
To The Variable "x" 


Jt 

Final Computations 

(Compute Extrapolated State Vector Corresponding 
To The Variable "x") 


t 

EXIT 


Figure 2 . Theta Routine, Functional Flow Diagram 
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:c input and output variab les 


The Conic Slate Extrapolation Routine has only one universal constant: 

(he gravitational paramoter of the earth. Its principal input variables are the t 
inertial stale vector which is to he extrapolated and the transfer time interval or' 
transfer angle 1 through which the extrapolation is to he made. Several optional 
input variables may he supplied in the transfer time ease* in order to speed the 
computation. The principal output variable of both rases is the extrapolated in- 
ertial state vector. 

3. I Conic State Extrapolation as a Function of Transfer Time 
Interval (Kepler Routine) 

Input Variables 


- 0 * 


A.t 


Inertial state vector which is to be extrapolated 
(corresponds to time t^). 

Transfer time interval through which the extrapola- 
tion is to be made. 


x 


At' 


x 1 


c 


Guess of independent variable corresponding to solu- 
tion in Kepler iteration scheme. (Used to speed con- 
vergence). j^If no guess is available, set x = 0, and 
the routine will generate its own guessj. 

Value of dependent variable (the transfer time inter- 
val) in the Kepler iteration scheme, which was 
calculated in the last iteration of the previous call to 
Kepler, 

Value of the independent variable in the Kepler itera- 
tion scheme, to which the last iteration of the 
previous call to Kepler had converged. 


Output Variables 


(r, v) 


Extrapolated inertial state vector ( corresponds to 
time t ). 
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A t 


Value of the dependent variable (the transfer time 
interval) in the Kepler iteration scheme, which was 
calculated in the last iteration (should agree closely 
with At). 


Value of the independent variable in the Kepler itera- 
tion scheme to which the last iteration converged. 


3.2 


Conic State Extrapolation as a Function of Transfer Angle 


( Theta Routine ) 


<I0' V 


(£, v) 

At 


Input Variables 


Inertial state vector which is to be extrapolated. 


Transfer angle through which the extrapolation is to 
be made. 


Output Parameters 


Extrapolated inertial state vector. 


Transfer Time Interval corresponding to the conic 
extrapolation through the transfer angle 6 . 
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4 . DESCRIPTION O ^ E QUA TI ONS 

4. 1 Conic State Extrr^.d, i> n as a Function of Time (Kepler 

Routine) 

— — 1 » 4 - 
<. ~ 
i, - 

' The universal foY-Yapon of Stumpff-Herrick-Battin in 
terms of the universal eccentric anomaly difference is. used. This 
variable, usually denoted by x, is defined by the relations: 

\| 1 P(E - E q ) for ellipse 
\l]T( tar./ y / ° - tan f Q / 2 ) for parabola 
,i T ^/ for hyperbola 

where a is the semi-major axis, E and H are the eccentric anomaly 
and its hyperbolic analog, p is the semi-latus rectum and f the true 
anomaly. The expressions for the transfer time interval (t - t Q ) = At, 
and the extrapolated position and velocity vectors (r\ v) in terms of 
the initial position and velocity vectors (r^, v^) as functions of x 
are: 


(Dn ~ Y- Kepler Equation) 








where 



and 



S<5) 





Since the transfer time interval At is given, it is desired to 
find the x corresponding to it in the Universal Kepler Equation, and 
then to evaluate the extrapolated state vector (r, v) expression 
using that value of x. Unfortunately, the Universal Kepler Equation 
expresses At as a transcendental function of x rather than con- 
versely, and no power series inversion of the equation is known which 
has good convergence properties for all orbits, so it is necessary 
to solve the equation iteratively for the variable x. 

For this purpose, the secant method (linear inverse inter- 
polation/ extrapolation) is used. It merely finds the increment in 
the independent variable x which is required in order to adjust the 
dependent variable At, to the desired value At based on a 
linear interpolation/ extrapolation of the last two points calculated on 

the At vs x curve. The method uses the formula 
c 


< x n+l -V = ' 


At c (n) - At 


At c (n) - At c (n " l) 







% 


% 



i 


J 

E 

i 

j 
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where At denotes the evaluation of the Universal Kepler 
Equation using the value x.. In order to prevent the scheme from 
taking an increment back into regions in which it is known from past 
iterations that the solution does not lie, it has been found convenient 
to establish upper and lower bounds on the independent variable x 
which are continually reset during the course of the iteration as more 
and more values of x are found to be too large or too small. In ad- 
dition, it has also been found expedient to damp by 10% any incre- 
ment in the independent variable which would (if applied) take the 
value of the independent variable past a bound. 

To start the iteration scheme, some initial guess of the 
independent variable is required as well as a previous point (x_j, 

A t ^ ) on the At c vs x curve. If no previous point is 
available the point (0, 0) may be used as it lies on all At c vs 
x curves. The closer the initial guess is to the value of x corres- 
ponding to the solution, the faster the convergence will be. One 
method of obtaining such a guess is to use a truncation of the 
infinite series obtained by direct inversion of the Kepler Equation 
( expressing x as a power series in At ). It must be pointed out 
that this series diverges even for "moderate 11 transfer time inter- 
vals At ; hence an iterative solution must be used to solve the 
Kepler equation for x in the general case. A third order truncation 
of the inversion of the Universal Kepler Equation is: 




n=0 


where 


*o = °* 
*2 ■ 


1 *<=o_=o )f 

2 r„ 3 ^ 


0 






6r 0 -0 


3 £o ' Xq ^ 

- - (1 - r„ 


Lr 0 \p"‘ 


v } 


with 


= 2 / r 0 - v 0 /" • 
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4.2 Conic State Extr apolation as a Function of Transfer Angle 

( Theta Routine ) 

As with the Kepler Routine, the universal formulation of 
Stumpff-Herrick-Battin in terms of the universal eccentric anomaly 
difference x is used in the Theta Routine. A completely analogous 
iteration scheme could have been formulated with x again as the in- 
dependent variable and the transfer angle 8 as the dependent variable 
using Marscher's universally valid equation: 


* 0 

cot — = 


1 - a Q x 2 S («qX 2 ) 


+ cot y 


2 x C ( x 2 ) 


0 


where 


and 


P . ,IoV si „ 2 




T 0 = angle from ^ to v Q . 


However, in contrast to the Kepler equation, it is possible 
to invert the Marscher equation into a power series which can be 
made to converge as rapidly as desired, by means of which x may be 
calculated as a universal function of the transfer angle 8. Knowing 
x, we can directly calculate the transfer time interval At^ and 
subsequently the extrapolated state vectors using the standard 
formulae. 

The sequence of computations in the inversion of the 
Marscher Equation is as follows: 


Let 

P N = P / r 0* °N = a r 0 


and 


V 
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W&fi* J S 


W 1 = ^r 1 B ■ cot 7 ° K 

* 1 - cos 6 


|Wj '>1. let V x = 1. 


w..., = + dw_ + «„ + WJ ( W. <1 


n N n 


V n + 1 = + fn 2+ “N ( l 1 / W J> + V n < W l >1) 


n n 


( Wj <1] 


l/“ n * <1 V w ! I>/ V „ ( Wj >1). 


* *- s 




“n j = 0 2 i + 1 w n 


where n is an integer > 4, Then 


x /\po 


(Wj > 0) 


2»/\p7 ' * < W !< 


The above equations have been specifically formulated to avoid certain 


numerical difficulties. 
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5 . 


DETAILED FLOW DIAGRAMS 


This section contains detailed flow diagrams of two Conic State Extra- 
polation Routines (Kepler and Theta) and the subroutines used by them. 

Each input and output variable in the routine and subroutine call state- 
ments can be followed by a symbol in brackets. This symbol identifies the nota- 
tion for the corresponding variable in the detailed description and flow diagrams 
of the called routine. When identical notation is used, the bracket symbol is 
omitted. 
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Call Marscher Equation Inversion Subroutine 


Input: 

6, cot V r Q , a N . P N 

Output; 

X, 4 . c v c 2 


Call Universal Keplei 

r Equation Subroutine 

Input; £ > x * 

Output: A i c , 

c l» c 2' r 0 
S<4). C (4) 

1 


, Call Extrapolated State Vector Subroutine 


Input: rQ, Vg, x. C, S(£), C(£) 

Output: r , v 


Yes 



P - 


2tt 


<» N '- 0 ) 3/2 f^ 



Figure 4b. Theta Routine, Detailed Flow Diagram 
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Subroutines Used By The Tran sfer T ime o^ransf erAngle^ 
„ Conic Extra polat i on Routines 


Universal Kepler Equation Subroutine 



Figure 5. Universal Kepler Equation, Detailed Flow Diagram 
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G su P PIGMENTAR Y infor mation 

The analytic expressions for the Universal Kepler Equation 
and the extrapolated position and velocity vectors are well known and 
are given by Battin ( 1964 ). Battin also outlines a Newton iteration 
technique for the solution of the Universal Kepler Equation; this tech- 
nique converges somewhat faster than the secant technique but 
requires the evaluation of the derivative. It may be shown that if the 
derivative evaluation by itself takes more than 44% of the computa- 
tion time used by the other calculations in one pass through the loop, 
then it is more efficient timewise to use the secant method. 

Marscher’s universal equation for cot 0/2 was derived by 
him in his report ( Marscher, 1965), and is the generalization of his 
"Three-Cotangent" equation; 




0 r o < E -V 

cot £ = cot — + cot r 


\fp~a 1 


0 





Marscher has also outlined in the report an iterative method of ex- 
trapolating the state based on his universal equation. The inversion 
of Marscher's universal equation was derived by Robertson (1967a). 

Krause organized the details of the computation in both 

routines, 

A derivation of the coefficients in the inversion of the Uni- 
versal Kepler Equation is given in Robertson (1967 b) and Newman 
( 1967), 
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Submittal 24: Required V e .>city Detetm ination , Corife ^ 


INTRODUCTION 


The Conic ROviJ I Velocity Determination Routine provides 
the capability to solve the following two astrodynamic problems: 


’’The Multiple-. R; ■ ol ution , Lambert Required Velo city Determination 
Problem 11 : up me /he velocity vector required at an initial posi- 

tion to transfer through an inverse square central force field from 
the initial position to a specified target position in a specified 
transfer time interval by making a specified number of complete 
revolutions (/:. ; tiome fraction of another one). Also optionally 
compute the velocity vector at the target position and various 

parameters rf the rrnic transfer orbit. 

V\- : - * 

"The De-orbit Required Velocity Determination Problem”: compute 
the velocity vector required at an initial position to transfer through 
an inverse square central force field from the initial position to a 
specified target radius (which is less than the initial radius) with a 
specified flight-path angle at that radius in a specified transfer time 
interval. Also optionally compute the velocity vector at the target 
position and various parameters of the conic transfer orbit. 

The Conic R- Velocity Determination Routine basically consists of two 

major parts — one i naming the multi-revolution Lambert's problem and one for 
solving the De-orbit problem — which are quite similar. In fact, certain subsections 
of the parts are identical as well as being identical to certain subsections of the 
Conic State Extrapolation Routine (Ref. 7) and these may of course be arranged as 
subroutines on a computer. 

The Conic Lambert and De-orbit Required Velocity Determination Routines 
each involve a single loop iterative procedure, and hence are organized in three 
sections: initialization, iteration, and final computations, as shown in Figure 1. 

The independent variable in the iteration in both routines is the cotangent of the 
flight-path angle a •. initial position measured from local vertical, or equiva- 
lently the cotangen N./.he angle between the initial position vector (extended) and 
the as yet unknown required velocity vector. The dependent variable is the trans- 
fer time interval; it is a function solely of the independent variable and certain 
other quantities which depend explicitly on the input and which are thus constant in 
any one problem. In the iterative procedure, the independent variable (denoted by 
r Q ) is adjusted betw upper and lower bounds by a secant technique until the 
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transfer time interval computed from it agrees with the specified transfer time in- 
terval (to within a certain tolerance). Then the velocity vector at the initial posi- 
tion (i. e. , the required velocity), as well as the velocity vector at the terminal 
position, is calculated from the last adjusted value of the independent variable. 

In the less-than-one-complete revolution case in both routines, the upper and 
lower bounds on the independent variable are explicitly computed since the depen- 
dent and independent variables are monotonically related. However, in the multi- 
revolution case in the Lambert routine, there are two distinct physical ly-meaning- 
ful transfers which solve the problem, and an iterative procedure (entirely separate 
from, and not containing nor contained in the previously described iteration scheme) 
must be used to solve for the value of the independent variable which separates the 
two regions in each of which exactly one solution lies so that upper and lower bounds 
may be established corresponding to the unique solution desired. The multi-revo- 
lution case for the de-orbit problem is not considered in this document. 
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ENTER 



! Initialization 

J 

{ Compute various constant parameters, 

! Establish range over which independent 
.] v triable may vary. Establish first 
i *value of independent variable. 


Iteration 


Compute Transfer Time Interval cor- 
responding to the independent variable. 


'Computed' 
Transfer 
Time = 
Specified 
Transfer 
\Time 


Adjust independent variable 
by Secant technique. 



Final Computations: Compute 
Required Velocity correspond- 
ing to independent variable. 


EXIT 


Figure F Conic Larnbert and De-orbit Required Velocity 

Determination Routines Functional Flow Diagram 
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NOMENCLATURE 


i 

1 






a Semi-major axis of conic 

Cj First conic parameter [ c J = -\fr oTn" ^0 = ^0 ' -0^^ E 

2 

c 2 Second conic parameter [ = 1 = r 0 V 0 / 1 

C or C(£) Power series in £ defined in the text 


E 


H 


l 

— c - 


max 


- N 


i 

~r r 


i 

“r, 


Eccentric anomaly 


True anomaly 


Hyperbolic analog of eccentric anomaly 


Iteration counter 


The negative unit chord vector connecting £q and 
r j • [ l c _ = -unit (r : - r Q )] . 


Maximum allowable number of iterations 


Unit vector in direction of angular momentum 
vector of the transfer and normal to the transfer 
plane. In the Lambert Routine the vector 
always determines the direction of the transfer, 
and will also determine the plane of the transfer 


when either the switch s . = 1, or the switch 

proj 

s . = 0 but the initial position vector r A is 
pro] y —0 

inside one of the cones. In the De-orbit Routine, 


the vector _i^ always determines the plane and 
direction of the transfer. 


Unit vector in direction of 


Unit vector in direction of £ ^ 
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4 


k 


k 


bg 


k 


sm 


m 


m 1 


m 


err 


n 

n 

rev 


N 

P 

Pi 


P2 





Intermediate variable equal to either or k gm 

Constant establishing by what fraction of its permissible 
range (P - P . ) the independent variable P, will 
be biased in the first iteration when no guess ]T^ uess is 
available, in order to establish a second point for the 
secant iteration 

Constant establishing by what fraction of its permissible 
range (P - P ) the independent variable will 
be biased in the first iteration when a guess I^ uess is 
available in order to establish a second point for the 
secant iteration. 

The slope of the line joining two successive points on the 
transfer time interval vs. independent variable curve. 

Previous value of m 

Difference between desired value of the slope m (namely 
zero) and the value calculated on most recent iteration. 

Loop counter in the Marscher Equation Inversion 

Integer number of complete 360° revolutions to be 

made in the desired transfer. [Hence the transfer 

will be between n and n + I revolutions]. 

rev rev J 

Intermediate vector variable normal to transfer plane 
Semi-latus rectum of conic 

Intermediate variable in the Lambert problem equal 
to 1 - cos 0 

Intermediate variable in the Lambert problem equal 
to cos 0 - (tq / r^ ) 

Normalized semi-latus rectum of conic transfer orbit 

(p N = p/v- 

2 

Intermediate variable equal to X/sin Yj 
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«L*i J iiw £.'• iftittiSi « 


Initial or current inertial position vector (corresponds 


to time 1 q). 


Terminal or target inertial position vector (corresponds 


to time t j ). 


Radius at terminal or target position (corresponds 
to time t , ). 


Switch used in Secant Iterator to determine whether 
secant method or offsetting (biasing) will be performed. 

Switch indicating whether the outcome of the cone test 
involving the tolerance criterion c cone was that initial 
position lies outside both of the cones around the 

positive and negative target position vector r ( s cone = 
0), or inside one of these cones ( s cone = 1 )■ [See 
Section 4,7.] 


guess 


Switch indicating whether the routine is to compute its own 
guess of the independent variable Tq to start the iterative 
procedure (s guess = 0), or is to use a guess T guess 
supplied by the user ( S g Uess “ 1) 


Switch indicating whether the initial and target position 

vectors, r Q and r ^ , are to be projected into the plane 

defined by the unit normal _i before the main Lambert 

computations are performed. If S p ro j = no P r °j ect i° n 

will be made unless the initial position r is found to 

lie within one of the cones defined by €' cone j * n which case 

s will be set equal to 1. If s . = 1, the projections 
cone P ro J 

will be carried out immediately, and no cone test will 
be made. 
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1 soln 


Switch indicating which of the two physically possible 
solutions is desired in the multi-revolution case. 

{ Not used in the less- than- 360° transfer case] , In 
particular, s so ^ n = -1 indicates the solution with the 
smaller initial flight path angle V ^ measured from 


local vertical, and s SQ i n 
the larger Vq . 


+1 indicates the one with 


180 


Switch indicating whether the central transfer angle 

is between 0° and 180° (s n = +1), or between 180° 

1 oU 

and 360 (s - -1 ), The determination of which 
180 

one of the above two possibilities is desired is made 
automatically by the routine on the basis of the di- 
rection of the unit normal vector P, . 

— N 

[In the multiple-revolution case, the number of com- 
plete 360° revolutions is neglected; i. e. , s g^ is the 
sign of the sine of the transfer angle, ] 


■m 


S orS(£) Power series in £ defined in the text. 


err 


Difference between specified time interval and that 

calculated by Universal Kepler Equation [ t err = 

At - At ]. 

c J 


m 


-o' 


Inertial velocity required at the initial position r Q 
to transfer to the terminal point in exactly the 
specified time interval At.. 


M 


a 




Inertial velocity at the terminal position r j 


Intermediate scalar variables used in Marscher Equa- 


(n-1, 2. . ) lion Inversi 


ion 


W Intermediate scalar variables used in Marscher Equa- 

n ^ 


(n=l,2. .) lion Invers 


ion 
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Universal eccentric anomaly difference corresponding 
to the transfer from r q to r ^ . 


Normalized universal eccentric anomaly difference 

(x n = x /V^T ] 

Reciprocal of normalized semi-major axis of conic 
transfer orbit = Tq / a). 

Flight path angle at initial position £q measured 
from local vertical, i. e. , angle from r Q to v Q . 

Flight-path angle at terminal or target position 
measured from local vertical (corresponds to 
time t j ) . 

Cotangent of flight-path angle y^ at the initial posi- 
tion £q measured from local vertical; i.e., 
cotangent of the angle between £q and Vq . [ In- 
dependent variable in iterative scheme] . 


Previous value of 


The "i-th" value of Tq 

Cotangent of flight path angle at the terminal 
or target position £ ^ measured from local 
vertical 


guess 


Guess of independent variable corresponding 
to solution (disregarded when S g uess = 0). 


' , Value of corresponding to the physically 

parab 0 

realizable parabolic transfer 


Upper bound on 

Value of r Q corresponding to the minimum energy 
transfer 


S24-8 







r min Lower bound on 


At 




Specified transfer time interval (t^ - t Q ) between Fq 

and £ j 

Value of the transfer time interval calculated in the 
Universal Kepler Equation from the current value 
of and the conic parameters 

Previous value of At 

c 


At ^ The "i-th" value of the transfer time interval cal- 

c 

culated in the Universal Kepler Equation as a 
function of the "i-th M value 1^ ^ of and the 
conic parameters 


* r o Increment in Tq 

AA Increment in A 

e Tolerance criterion establishing small cones around 

cone 

both the positive and negative target position direc- 
tions inside of which the Lambert routine will define 
the plane of the transfer by the unit normal 
rather than the cross product of the initial and target 
position vectors, r Q and r . f e cone ’ s * n 
cone angle ) ] . 


e t Primary convergence criterion: relative error in 

transfer time interval 

€ t' Secondary convergence criterion: minimum permissible 

difference of two successive calculated transfer time in- 
tervals. 


e j, Convergence criterion in iteration to adjust - n 

and r in multiple revolution case: absolute 
max r 

precision to which transfer time interval minimum 
is to be determined 

Tertiary convergence criterion: minimum per- 
missible size of increment Al^ of the independent 
variable 





ri 


i 
. ? 



0 

X 

A 

A’ 


■- T oierance criterion in iteration to adjust TT and 
} ... in multiple revolution case: absolute dif- 
/ierence of two successive values of independent vari- 
able to prevent division by zero 

angle (true anomaly increment) 

Ratio of initial position radius to terminal position 
radius 

Average of the two most recent values of * [A 
is used as the independent variable in the Multi- 
r^-'t/iution Bounds Adjustment Coding Sequence 
IH' -iticn] . 

Previous value of A 


P Gravitational parameter of the earth (product of 

earth’s mass and universal gravitation constant) 

. 2 2 2 
q r he dimensionless variable ax = x /a = x / 

r . [ Equivalent to square of standard eccentric 

or perbolic anomaly difference] . 




> 



I 
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5. D ETAILED FLOW DIAGR AMS 

5. 1 Multiple- Revolution L a mbert Req u ired Velocit y 

Determ i n at iorFfr ou t ine 

This routine utilizes the following subroutines or coding sequences, which 
are diagrammed in Section 5. 3: 

• Lambert Transfer Time Interval Subroutine 

•Marscher Equation Inversion Subroutine 

• Universal Kepler Equation Subroutine 
•Secant Iterator 

• Multi-revolution Bounds Adjustment Coding 

Sequence 

• Secant Minimum Iterator 


UNIVERSAL PROGRAM 
CONSTANTS CONSTANTS 


INPUT VARIABLES 


€ . 1 , 6 r A , r 1 , At, n , s , , 

t t r T A — 0 — 1 rev soln 

^max’ ^sm’ ^bg S guess' ^guess’ cone’ proj 7 -N 


r 0 

\Lq\- r i = HP 

i = 

r 0 

— 0 / r 0 ’ i-r 1 = i 1 / r 1 

N = 

i X i 

“ r 0 ~ r 1 

j sin ej = 

| N | 


sin eke. 


s = i 
cone 


s =0 
cone 


Figure 3a. Multi-Revolution Lambert Routine 
Detailed Flow Diagram 
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5. 2 De-orbit Required Velocity Determination Routine 

This routine utilizes the following subroutines which are diagrammed in 
Section 5. 3: 

• De-orbit Transfer Time Interval Subroutine 

• Marscher Equation Inversion Subroutine 
° Universal Kepler Equation Subroutine 


• Secant Iterator 
UNI VE RS A L P ROG RA M 

CONSTANTS CONSTANTS INPUT VARIABLES 




Figure 4a. De-orbit Routine 

Detailed Flow Diagram 
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5. 3 Subroutin es o r Codin g Sequences used by the Conic 
Requxr ; V' 1 cTvy Determination ITouttnes 


5. 3. 1 Lambert Transfer Time Interval Subroutine 


CC " ryVA^TS INPUT VARIABLES 


. r 

M I! 


Iq. Pj, p 2 . sin 0, cos 0, r n , n 


0' rev 


IT sin & - p, 


2 -P N (1+r o 1 


Call Marscher Equation Inversion Routine 
(Input: sine, cos 0, Tq. r Qt q- n , p N 

Output: x ( c j , c 2 > 


Ca 1 \ T r- • .rsal Kepler Equation Routine 
C y , c 2 - X , 4, r 0 


Output: A t c , S(4), C (4) 


n = CP 
v rev > 


At = At + n 2Tr“\/jr n /Qf KT j /p 

c c rev v » O' N I ' ~ 


G' U N 


OUTPUT VARIABLES 


At c* a N' ^ N J £ * S({), C (4) 


Figure 5. Lambert Transfer Time Interval Subroutine 
Detailed Flow Diagram 
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5* 3. 2 De-orbit Transfer Time Interval Subroutine 


INPUT 

VARIABLES 



Figure 6. De-orbit Transfer Time Interval Subroutine 
Detailed Flow Diagram 
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5. 3. 3 Universal Kepler Eouation Subroutine 




This subroutine is identical to the one used in the Kepler and Theta problems. 


UNIVERSAL INPUT 
CONSTANTS VARIABLES 



§ 


Figure 7. Universal Kepler Equation Subroutine 
Detailed Flow Diagram 
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5. 3. 4 Marscher Eqy n Hon Inversion Subroutine 

This subrou*i;< i identical to the one used in the Theta problem , 



** - : PUT VARIABLES 



(No 

physical! 

realizab! 

solution 

possible' 













































At < At 


soln 


r av - 

A 

max 

^min 

max ( ^min' 
Iparab) 

k bg 

_k bg 

k - 

-k 

s m 

sm 


^rnin ^ 


EXIT 

(to Lambert Rou- 
tine, Figure 3b) 


ERROR EXIT 

(No solution possible to 
this Lambert problem: 
too many revolutions 
for too short a speci- 
fied transfer time in- 
terval) 


Figure 10c. Multi- Revolution Bounds Adjustment Coding Sequence 
Detailed Flow Diagram 
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5, 3. 7 Secant Minimum Iterator 


This subroutine is very similar, though not identical, to the Secant 
Iterator. They can easily be combined into one routine, although they have 
been diagrammed separately here for purposes of clarity, 

INPUT VARIABLES 

^ m - m ’’ m err' A A A- r min - r max , T1 


m 

aA = — ' p - ee -t a a 

1 m - m 1 

AA = sign (in err ) |aA| 


AA = 

sign (m err ) 


k ( r -L:J 


m ax mm 

S - 

0 


A+ aA 


A +AA 


A A = o.9 ( r -A) 


AA = o. 9 ( r -A) 

n * rrn n * x 


OUTPUT VARIABLES 


AA- s 


Figure 11. Secant Minimum Iterator 
Detailed Flow Diagram 
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Subm ilt a! 2 7 : Re quired Velocity De termination, Precision 
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1. INTRODUCTION . 

Calculation of the precision required velocity which satisfies terminal 
position and time-of-flight constraints in a non-Keplerian gravity field is a com- 
putation time consuming process, especially in an on-board computer. There- 
fore, targeting calculations prior to a maneuver are customarily used to predict 
and compensate for the effects of the perturbations from a conic gravity field, 
so that during the maneuver only the much simpler conic related computations 
will have to be performed. 

For Lambert aim point maneuvers (described in Reference 2) an 
adjustment to the terminal (target) position vector will suffice to provide this 
compensation. This adjusted terminal position, referred to as an offset tar- 
get, must compensate for gravity perturbations throughout both the maneuver 
and subsequent coasting flight. Then the required velocity determined by the 
Lambert routine to intercept the offset target in a conic gravity field is 
identical to the velocity required to intercept the true target in the non- 
Keplerian field. 

The traditional technique of predicting the effects of gravitational pertur- 
bations over the trajectory involves approximating the maneuver by an impulsive 
velocity change, and hence assuming a coasting trajectory between the initial 
(ignition) and target positions. However, due to the non-zero length of the maneu 
ver, the actual trajectory will not follow the path predicted by the impulsive 
approximation, but rather a neighboring path. The difference in the perturbing 
acceleration between the two paths accumulates over the entire trajectory, re- 
sulting in a miss at the target. Since the coasting portion of the trajectory is 
generally much longer than the thrusting portion, it is important to accurately 
predict the perturbing effects over this portion of the trajectory. This is accom- 
plished by determining the initial conditions for a coasting trajectory which is 
coincident with the actual trajectory after thrust termination. A detailed deriva- 
tion of this technique can be found in Brand (1971) (Reference 1), and a functional 
description of the procedure follows. 
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2 , FUNCTIONAL FL DIA GRAM 

A functional flow diagram describing the calculations necessary to 
determine the precision rr ,;.’ velocity and offset target is presented in Figure 
1 Since this technique c s for the non-impulsive nature of the maneuver, 
it requires an estimate of.':.;* t ted thrust acceleration. Then the initial 

position can be offset froxL 'iv^ actual position such that a coasting trajectory 
which is coincident with the actual trajectory after, thrust termination can be 
defined. Figure 2 illustrates the concept. 

The calculation ol the coasting trajectory initial position requires an 
estimate of the required velocity change, and therefore two passes are made 
through the Lambert rout*;* ocT re numerically integrating to determine the 
effects of gravitational pei lui’ oauons. The first Lambert solution is used to 
determine the impulsive velocity change required. Based upon this, an estimate 
of the initial position for the coasting trajectory can be calculated. Then the 
second Lambert solution determines the velocity required from the adjusted initial 
position, thus defining the coasting trajectory. 

For transfers angles which are odd multiples of 180 , Lambert's problem 
has a partial physical singularity in that the plane of the transfer becomes inde- 
terminate. A detailed de, Ton of this singularity can be found in Reference 4, 

To prevent possible probu ' in both targeting and guiding a maneuver whose 
transfer angle lies near this singularity, logic has been included in this routine 
to determine whether the transfer angle approaches this singularity at any time 
during the maneuver. If this is the case, the target vector is projected into the 
orbital plane defined by the premaneuver position and velocity, thus preventing 
any plane change. 

If only conic calculations are desired, the routine is exited after the two 
Lambert solutions are completed. If not, subsequent numerical integration de- 
termines the target miss resulting from the effects of gravitational perturbations 
over this path. To compe nsate for these effects, the target vector for the Lambert 
routine is offset from the actual target by the negative of the miss vector. Since 
the adjusted initial position, target offset, and effects of gravitational perturba- 
tions are all interdependent, the process is repeated until changes in the offset 
target position are small enough to indicate convergence. Three passes (two 
iterations) are normally sufficient to establish the offset within a few feet. 


A 

j 
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ENTER 


Estimate thrust acceleration. 
Initialize switches (set s pro; j to zero). 


Initialize iteration counter. 

Set initial value of offset target equal 
to actual target. 


Compute adjusted initial position based upon 
required velocity change and thrust accelera- 
tion (no adjustment on the first pass). 


Use Lambert routine to compute velocity re 
quired to transfer from adjusted initial posi 
tion to offset target. 


if transfer angle is near 180 during the maneuver, 


Numerically integrate the required velocity irom 
the adjusted initial position through the. specified 
time of flight, including gravitational perturbations. 


Figure la. Functional Flow Diagram 
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Figure lb. Functional Flow Diagram 









Figure 2, Coasting Trajectory Illustration 




NOMENCLATURE 


a 


T 


d 


f 

f ACS 

f OMS 


-N 


m 

n 

n 

max 

n rev 


lo 


£o 


II 


-lc 


cone 


Estimated m /.’..d-'t -id- [ the thrust acceleration 

Number of co’lmh:. ♦./ navigation filter weighting matrix (set to 0 in 
this routine since the matrix is not required) 

Thrust 

Magnitude of the attitude control system translational thrust 

Magnitude of ”.-r *.o:v- jial orbital maneuvering system engine 
thrust 

Unit normal to the trajectory plane. (in the direction of the 
angular momentum at ignition) 

Current estimated vehicle mass 

Iteration count 

Iteration limit 

Integral number of complete 360° revolutions to be made in 
the desired transfer 

Initial (ignition) position 

Adjusted initial position used to define coasting trajectory 
Target position .put to the routine) 

Terminal position (output of the routine) 

Offset target position 

Switch set in the Lambert routine to indicate transfer is near 
180° (see Reference 4 for complete description) 
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T'JS/JF ..-W-L *UA L 

• - • v ■■ ; 



Engine select switch 


Switch set to indicate non-convergence 

Switch set to indicate an estimate of independent variable T 
will be input to the Conic Required Velocity Determination 
Routine __ 

Switch set to indicate which perturbing accelerations should be 
included in the offset target calculation ( s p er ^ = 0 indicates 
only conic calculations; see Reference 3 for complete descrip- 
tion of other switch settings) 

Switch set when the target vector must be projected into 
the plane defined by 

Switch indicating which of two physically possible solutions 
is desired in the multi ^revolution transfer (see Reference 4 
for complete description) 

Igpition time 

Target time of arrival' 

Initial (ignition) velocity 

Initial (and required) velocity on the coasting trajectory 
Terminal velocity of a conic trajectory 
Terminal velocity (output of the routine) 

Guess of the independent variable T used in the Conic Required 
Velocity Determination Routine 

Target miss resulting from perturbations 

Out-of-plane target miss due to projection of the target 
vector 
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^ v , 


‘ conv 


Transfer time - t Q ) 

Required velocity change 

Magnitude of the required velocity change 

Convergence criterion: target miss of the numerically integrated 
trajectory 


•a 


f eT 


U> 


Tolerance criterion establishing a cone around the minus 
r direction inside of which the target vector will 
be projected into the plane _i N |e eT = sin (half cone angle) j 

Transfer angle (true anomaly difference) at the start of the 
thrusting maneuver 

Approximate central angle traversed during the thrusting maneuver 

Approximate transfer angle to the target at the termination 
of the thrusting maneuver = 9 - 

Approximate orbital rate 


J 
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UNIVERSAL 

CONSTANTS 


PROGRAM 

CONSTANTS 


INPUT 

VARIABLES 


f OMS' f ACS 


e conv' n max’ 


v V V V El- 

n rev' S Boln' S pert' J eng 


Set f according to s 


= i/m 

■ l l - l o 
- hol/llol 

= unit (r Q x ^q) 


guess 


n = 0 


Elc* El 

-0 = x o 


Eo “ Eo 


Call Conic Required Velocity Determination 
Routine (Reference 4) 

Input: [Eq] ' lie 

M ■ 

At< n rev* s S oln* 

s guess’ r guess 

' C 6T 

r 

[ € cone] ' S proj' — N 

Output: , £i c |£]j 

• ii[ 


T S i 

A guess, cone 

. T 

sin d. 

cos d 


Figure 3a. Detailed Flow Diagram 
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Submittal 51: Boost Abort Guidance 

Equations and flow diagrams are presented in this Section which 
fulfill requirements for abort from boost to an entry path which. achieves 
satisfactory landing at the launch site. Constraints and guidelines are 
presented in Fig. 1. The trajectory and nomenclature are presented in 
Figs. 2 and 3. A general flow diagram is presented in Figs. 4 and 5. 

As shown in Fig. 2, landing is achieved in four phases; An open loop 
phase of powered flight wherein propellant is expended, powered-flight 
constraints are observed and conditions are reached where available AV 
equals AV required to get on the entry trajectory. A closed-loop phase 
of powered flight achieves entry target conditions with very little fuel. 
An unpowered flight phase follows wheje unpowered flight constraints are 
observed and the trajectory approaches the nominal entry trajectory. The 
final unpowered phase consists of holding the entry trajectory, i.e., 
controlling to the trajectory through satisfactory landing. These phases 
and their Guidance Equations are presented in Figs. 6 thru 11. 

Figs. 12 thru 18 are the detailed program flow diagrams. Thereafter 
follows a definition of terms used in this section. 
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BOOST ABORT - CONSTRAINTS/GUIDELINES 
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AVOID LARGE NORMAL DECELERATIONS (2.5Q AND 
DYNAMIC PRESSURES <300 PSF) 




100 200 
RANGE FROM LAU 



Figure 3 


PLANAR REPRESENTATION OF ABORT MANEUVER 
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POWERED PHASE OF SUBORBITAL MANEUVER 
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F-ifrure 7 

PHASE I GUIDANCE EQUATIONS 


o RADIAL t PLATFORM X) 

ATC(1) = K-j CRDG - RG) + K 2 (VDG(1 ) - VG(1)) - G^pj 


o RETURN TO LAUNCH SITE 
ATC(2) = 0 

ATC(3) = AH AH = (AT 2 - ATC (1 > 2 )i 

o DOWNRANGE VDG(3) > VBOZD 

AGC(2) r o 
AGC(3) = AH 

o DOWNRANGE VDG(3) < VBOZD 

AGC(2) = AH . VG(2)/|VG(2)| 

AGCO) - 0 


o TRANSFORMATION , 

ATC(2) = AGC(2)UYGP(2) + AGC(3)UZGP(2) 
ATC(2) - AGC(2)UYGP(3) + AGC(3)UZGP(3) 






F igure 8 
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-Si 



Pl ; - 2 GUIDANCE EQUATIONS 

AGC(1) " tSo? ‘ j§Q ( VDG (D + 2VG(1 ) ) - GEFF 

AGCO) = CVDGC3) - VG(3))/TB0 

AGC(2) = YSIGN(Affy- AGC(l) 2 - AGCO) 2 )® 

\ . 

ATC = AGC (1 ) UXGP + AGC (2 ) UYGP 4 AGC(3) UZGP 
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UNPOWERED PHASE OF SUBORBITAL MANEUVER 
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Figure 10 

PHASE III GUIDANCE EQUATIONS 



u c - e D ’ ^ 

oc c = cx - Q- At 


9n < 9n LIMIT 
9n > 9n LIMIT 


BUT 

^MIN < ^ c MAX 


i 



0 \i c 
^vc 
#vc 

RATE LIMIT 


-15° • VG(2)/|VG(2)| 
-7.5° • VG(2)/|VG(2)| 
0 

,5°/SEC IN PITCH 
J& v IS TBD 




Yv > Yl 

Yi > Yv > Y2 
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D?'?i:r:Tio:j or symbols 


En 1 3 y gui d an cc t ok c c ve r 
Solid rocket r.ct or 


/cc 

;si 

ALPHA 

AbPHAC 

/PC 

EA'tK 

LAICKC 

D7A 

DVZBO 

G 

p ; > r\ t 
U r:_ .*i/i 

GFFF 


Vr^ablcs 

Thrust acceleration co:rrand in guidance cccrdi nater^ jpss 
Horizontal acceleration, FPSS 
Angle-of-nttack, deg 
CoEaanded engle-of -attack, deg 

Tiuast acceleration coasand in platform coordinates, FPSS 
Bank angle, roll about velocity vector, deg 
.Bank angl e connand, deg 
Delta-Y aval labia, FPS 

Delta-V required, velocity to be gained, FPS 

Predicted component of bumout velocity ndnus previous 
value 

Acceleration due to gravity, FPSS 

Relative flight rath angle, deg * ‘ 

Effective gravity acceleration, FPSS 


Grill 

Predicted effective gravi 
powered phase, FPS 5 

tv acceleration at terra, nus 

GLAT 

Latitude of abort landing 

site, degrees 

gi/j:;g 

longitude of abort lor dir. 

g site, degrees 

CMjOT " 

Tine rate of enlarge of no 

rr.sl acceleration, I : ?3S- 


'tOTZ al acsel s-rf-cioa, Fi'£S 


tj 
l A 

Altitude, ft 

. * 

KDOT 

Altitude rate, FPS 
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DEFINITION OF SYMBOLS 


kach 


Variables (continued) 
Velocity heading enrl o, degrees 
Launch azimuth, deg 
Mach number 


FITCH 


PJTCHC 


Q3AP. 

QBIOT 


RISE 


SF30 


SHC5T 


Plat Corn pitch, degrees 
Platform pitch command, degrees 
Intermediate variable 
Intermediate variable 
Intermediate variable 
Dynamic pressure, PSF 

Time rate of change of dynamic pressure, PSFS 
Altitude used in guidance equation. Ft 

Position vector to abort landing site in geographical ' 

coordinates, Ft 

Pc si tier, vector to abort lahding site in ir.ertial (geocentric) 
coordinates, Ft • 

Position vector to abort landing site in platform coordinates, Ft 

Position vector to vehicle in platform coordinates, Ft 

Surface range to burnout, Ft 

Surface range from burnout to abort lending site. Ft 


surface range from 


site, ft 


current vehicle position to abort land! nr 


Tine-to-go until burnout (fuel depletion), Sec 
Tir.e-f rom— turnout, Esc 
Time-f non-liftoff, Sec 

lime to rotate mom current acceleration vector to commanded 

c rys* r» *1 r-V- 
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DEPIMTIO!? OF 52K3CLS 


Variables (Concluded) 


TTL 

Time to thn 

no 

Time require 

UXGP 

Unit vector 

UYGP 

Unit vector 


coordinates 

UZC-P 

Unit vector 


VZBO 

VZBCD 

VZ30U 

W 

WDO? 

WDIR 

VSi? 


Tine required to null lateral (VG(2)) velocity 

Unit vector representing X-guidance axis in platform cocrdin 


^ , - w ■ - — - — » — V 1 V X ,U 

coordinates 

Exhaust gas velocity, FPS 

Relative velocity in guidance coordinates, FPS 
Inertial velocity in platform coordinates, FFS 
Relative velocity in platform coordinates, FPS 
Intermediate variable used in targeting 

A desired value of burnout velocity used in phase I guidance 
Intermediate variable used in targeting 
Weight of vehicle, Lb 

Jime rate of change of vehicle weight, Lb/sec 

Wind vector in platform coordinates, FPS 

Empty weight of orbiter, lbs 

Weight of 0M3 propellant, lbs 

Weight of orbiter (including main propellant) 
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ALIM 

AKAX 

AMIN 

ASZP 

BANKO 

BAIIK1 

bajokp. 

BOOT 

BDOT2 

BSE? ' 

Cl 

02 

C3 

DTI 

DT2 

GLAT1 

GLAT2 

GL0NG1 

glong2 

G'TLIM 

EB01 

E3C2 

HB03 


definition of stisols 

Variable Constants 

Thrust acceleration lirdt, FPSS 

Maximum ar.gle-of -attack licit, radian 

Minimum. angle-of-attack licit, radian 

Desired angle-of-attack for orbiter-tank separation 

Bank. angle commanded during orbi ter' maneuver away from 
tank (45 deg) 

Possible magnitude of bank angle command during phase III 
guidance (15 deg) 

Possible magnitude of bank angle command during Phase III 
guidance 

Bank angle rate limit during orbiter-tank separation 
Bank angle rate limit during entry (phase III) 

Desired bank angle for orbiter-tank separation 

Coefficients used in burnout targeting routine 

Powered guidance (phase I and II) cycle time, 1 • sec 
Unpowered guidance (phase III) cycle time, .2 sec 
Latitude of landing site near launen pad 
Latitude of a dovnrange landing site 
Longitude of landing site near launch pad 
Longitude of a downrange landing site 
Normal acceleration limit ( 1.8 G) 

Burnout, altitude for return to the launch site, 200,000 ft ■ 
Burnout altitude for "downrange landing site, 250,000 ft 
Booster - bum cut altitude, TBD 




1 

1 


i 


M 




A 


r i 




A 


$ 

i 

A 

'-•I 


1 

I 


HEAD1 

HEAD2 

HGAUT 

HI MU 

IP 

K1 

K2 

K3 

K4 

KU 

PDOT 

PDOT2 

PITCH1 

PITCH) 


PQPJJM 

Q3LIM 

PJ)G 

RE 

SHI FT 1 
SHIFT2 
SICFI3 


JL *ix 1 -.Oil k. 4 

Variable Constants 


STVE0L3 

(continued) 



If HEAD EAD1 3°, BAHX C - BAXXI 

If HEAD 2 (1°) HEAD HEAD1, BAh’KC = BAUK2 


Height of gantry 

Distance fron I MU to tail of orbiter 

Inertial to platform transformation matrix 

Altitude gain used in phase I guidance 

Altitude rate gain used in phase I guidance 

Altitude gain used in final portion of phase II guidance 

Altitude rate gain used in final portion of phase II guidance 

Earth gravitational constant 

Pitch rate limit during orbiter tank separation 
Pitch rate limit during orbiter reentry 

Platfom pitch required to- clear gantry for orbiter pad abort' 

Desired angle between X body axis and local horizontal plane 
during entry after powered abort maneuver 

•Attitude rate limit during powered abort 

Dynamic pressure Unit used for guidance purposes cnly 

Desired burr. out altitude, ft 

Radius of the earth, ft 

Shift in latitude, degree 

Shift in longitude, degree 

Shift in latitude, degree 



SHIFT4 

SMIIV'I 

CW %w O 


Shift in longitude, degree 

Acceptable angle-of -attack error for oi’biter tank separation 
Acceptable bark angle error for orbiter tank separation 
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TBOSW 


VE1 

VE2 

VXB01 


VXL'02 


VX303 

UbO 


Y/E 


VEP 

UORB 


DEFINITION OF SH30LS 
Valuable Constants (concluded) 

; - V % - • 

TBO for guidance ecuations (TO sec) 

>v?T. •- ; _ * 

Exhaust JlV 'city for booster 

Exhaust gas velocity for ortiter 

Desired alt* *r . de rate at bumout for return to the launch 
site (0) IFF ' ' . .. 

Desired altitude rate at bumout for dovnrange landing 
site (0) t:':' 

\ ’ v .. 

Desired altitude rate at booster bumout orbiter TBD 
Desired orbiter burnout weight weight, TBD 
Angular rotation rate of earth, deg/sec 

Angular velocity of earth in platform coordinates, aeg/see 
Orbiter li R -off weight, lb 

(for psrcT- r .. . im. weight of orbiter after separation from 
booster p .variable) 
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Flags 



* 

J 

J 


I 



>3ST 

UMT 

i’J.'OMIS 

K0H3IT 

KPASS 

KTKASS 

13R0P1 

KF0P2 

nCRIT 


Set to 1 when booster 1 3 attached 

Set to 1 if initial conditions for booster-orbiter ar8.used 
Set to 1 if nominal mission can not be continued 
Set to 1 if orbit can be achieved . 

Flag used to insure one pass through o. loop 

Flag indicating guidance phase 

Set to 1 when booster propellant is depleted 

Set to 1 when orbiter propellant is depleted 

Set to 1 for return to the launch site following an abort 

Set to 1 if abort is time-critical 
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